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Abstract 

The finite and high temperature equihbrium properties of Yang-Mills theory in Landau 
gauge are studied. Special attention is paid to the fate of confinement and the infrared 
properties at high temperatures. The method implemented are the equations of motion, 
the Dyson- Schwinger equations. A specific approximation scheme is introduced, which 
was previously applied successfully to the vacuum. 

In a first step, the infinite temperature limit is taken. The theory reduces to a 3- 
dimensional Yang-Mills theory coupled to a massive adjoint Higgs field. The equations for 
the propagators of the Higgs, the gluon, and the Faddeev-Popov ghost are obtained. They 
are solved in the infrared analytically and at least in the Yang- Mills sector confinement 
is found. Therefore the high-temperature phase of a 4-dimensional Yang-Mills theory is 
non-trivial and strongly interacting. Solutions for all propagators are obtained numerically 
at all momenta. Thereby also the propagators of a pure 3-dimensional Yang-Mills theory 
are determined. Systematic studies find only quantitative effects of the errors, which are 
induced by the approximations. Good agreement to lattice calculations is found. 

Finite temperatures down to the regime of the phase transition are investigated. It 
is found that the infrared properties are only quantitatively affected, and confinement of 
gluons transverse to the heat bath is estabhshed. The hard modes are nearly inert even 
at temperatures of the order of the phase transition temperature. Therefore the infinite 
temperature limit is a good approximation already at temperatures a few times the critical 
temperature, in agreement with lattice calculations. 

Finally quantities derived from the propagators are studied. The Schwinger functions 
are calculated. It is found that also the gluons longitudinal with respect to the heat bath 
are strongly influenced by higher order or even genuine non-perturbative effects, even in the 
infinite temperature limit. The analytic structure of the gluon propagator is investigated, 
and it is found that at least gluons transverse to the heat bath comply with the Kugo- 
Ojima and Zwanziger-Gribov confinement scenarios. Investigating the thermodynamic 
potential, an approximate Stefan-Boltzmann-like behavior is found. The thermodynamic 
potential, but not necessarily the pressure, is dominated by the hard modes. 

By comparison with calculations below the phase transition and lattice calculations 
it is conjectured that Yang-Mills theory likely undergoes a first order phase transition, 
which changes a strongly interacting system into another. The phases differ mainly by the 
properties of the chromoelectric sector. 
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Note added to the e-print version 

This thesis has been submitted to the faculty of physics of the Darmstadt University 
of Technology on the 19*^^ of October 2004. It has been defended successfully on the 

of December 2004. The supervisor was Prof. Jochen Wambach. This e-print version 
has some formal changes compared to the accepted version, e.g. the german abstract 
has been removed. The accepted offical version can be obtained from the Universitats- 
und Landesbibliothek of the Darmstadt University of Technology online. The URL is 
"http://elib.tu-darmstadt.de/diss/000504/". 

Most parts of the chapters El and IHl have been published in |HII. The unsettled problem 
discussed in subsection 15.5.21 has been resolved since the submission of this thesis. The 
method implemented here has been demonstrated to yield the correct limit for the number 
of Matsubara frequencies N going to infinity, N ^ oo. This will be discussed in some 
more detail in an upcoming publication by A. Maas, J. Wambach, and R. Alkofer, as will 
be most of chapter El 
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Chapter 1 
Introduction 



1.1 Strong Interactions 

Nature is described by only three fundamental forces according to the current knowledge. 
Each of those covers an own realm of physics and each one comes with its own problems. 
These forces are the electroweak and strong interactions, building together the standard 
model of particle physics, and gravitation [HEI- 

Gravitation describes the behavior of macroscopic objects to the largest distances 
known and it determines the gross properties of the universe at the present time. Its 
formulation in the general theory of relativity [2] has been supported by overwhelming 
experimental evidence. It is the weakest of the forces and quite well understood in the 
classical regime. On the other hand, there is still no successful implementation of a quan- 
tum theory of the tensorial gravitation field. Also several observations indicate that the 
current state of the universe is not only determined by the matter fields of the standard 
model. The nature and interactions of approximately 95% of the universe are not known 
today. 

Electromagnetism, electroweak symmetry breaking, weak processes like the /5-decay, 
neutrino-oscillations and various other effects are due to electroweak interactions as for- 
mulated by the Glashow-Salam- Weinberg theory |1J. It is rather well understood and its 
treatment using perturbative methods has been tested with high experimental precision. 
Although being rather tractable in general it is demanding in detail. Fundamental ques- 
tions are still posed by the absence of the Higgs in experiments up to now and the origin 
of the large number of parameters. 

The interaction binding together nucleons to build nuclei, thus forming the core of 
atoms, is known as the strong interaction. It is the strongest of all known forces. It 
describes also the way in which nucleons and all other hadrons are made out of their 
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1.1. Strong Interactions 



constituents, the quarks and gluons. The theory of strong interactions in its quantized 
form is termed quantum chromo dynamics (QCD) P], as quarks and gluons carry a charge 
termed color. In contrast to the electroweak theory, which includes electromagnetism, and 
gravitation, the strong force does not appear on a macroscopic level beyond its bound 
state spectrum. 

QCD offers a rich set of phenomena, which are not yet really understood. Primarily, 
a complete first-principle calculation of the bound-state spectrum of QCD is still lacking. 
Such a calculation must be able to explain why some bound-states like mesons and baryons 
appear in large numbers but why others like glue-balls and hybrids have not yet been 
convincingly found or appear only in rather small numbers like penta-quarks and meson- 
molecules possibly observed recently [Sj. 

At a more fundamental level, three properties of QCD are most striking. Two are 
connected to the fermionic content of the theory: The breaking of chiral symmetry and the 
axial anomaly. The first phenomenon gives rise to the proton mass of nearly 1 GeV, while 
the quarks making it up have masses only of the order of MeV. This can be understood 
as due to the spontaneous breaking of the approximate chiral symmetry of the lightest 
quarks and is also well known in other models. The axial anomaly is connected to a true 
quantum anomaly, and e.g. gives rise to the anomalous large mass of the T^'-meson. 

The third property is confinement, the absence of the colored degrees of freedom from 
the physical spectrum. It is this property which significantly shapes the low-energy reality 
of daily life while simultaneously being least understood of all the genuine non-perturbative 
effects of QCD. At the same time it is one of the properties which has been measured with 
the highest precision available: Free quarks have a unique experimental signature due to 
their fractional electric charge. The absence of such objects in nature has been established 
at a precision of the order of 1 : 10^" P^j. 

The true non-perturbative nature of QCD is the reason for the complications in the 
study of these phenomena. Perturbation theory describes strong effects reasonably well 
only at energies larger than a few GeV, the precise scale depending on the process. At 
smaller energies, which ultimately govern hadron and nuclear physics, the interaction is 
so strong that perturbation theory is not applicable. Therefore different approaches are 
needed. 

Extensive results are available from model calculations. However, only few attempts 
of first-principle approaches have been employed. Out of these, lattice gauge theory ^ 
is by far the most successful and widespread, and several achievements have been made 
using it. Lattice calculations are a numerical tool which discretizes space-time and thus is 
able to calculate the partition function directly, albeit numerically expensively. Especially 
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fermionic contributions pose significant problems. Furthermore lattice calculations are 
restricted at best to two orders of magnitude in momenta and rather small volumes due 
to the limitation by computing time. 

As many interesting phenomena of QCD are likely to be generated in the far infrared 
and include singularities, a continuum formulation is desirable. It is provided e.g. by the 
equations of motions in form of the Dyson-Schwinger equations. It is this method which 
is implemented here and it will be discussed in more detail in chapter IHl Compared to 
lattice calculations, this approach suffers from a larger number of necessary assumptions 
to make it technically treatable. At the same time it allows the direct investigation of 
the continuum manifestation of non-perturbative effects. Both these methods and several 
other approaches are complementary and necessary to untangle the complex structure of 
low-energy QCD. 

I. 2 Thermodynamics of QCD 

The phase structure of QCD is a long-standing problem [Hj. Naively, deconfinement would 
be expected in a sufficiently hot medium in a simple picture like the bag- model [Zj. The 
properties of such a different phase of QCD would grant additional insight into the mech- 
anisms of strong interactions. Such a medium is also expected to have existed at the 
beginning of the universe, and its specific properties are of relevance for cosmology. There- 
fore great experimental efforts have been undertaken in the last decades to generate it. 
Temperatures of around 100-200 MeV have been reached in heavy-ion collisions, but with 
large experimental uncertainties on the highest temperature achieved, especially as the 
thermalization process is not yet clear 0. 

This has also led to great theoretical efforts to solve the many-body problem of QCD. 
However, up to now the results generate more questions than answers. This is especially 
true for the non-equilibrium and dynamical part of the experiments, but even for the 
equilibrated medium it is still the case. 

The general aim of experiment and theory is to map out the phase diagram of QCD in 
the temperature- (baryon) chemical potential plane. Roughly sketched it is shown in figure 

II. H together with the currently available experimental data. The hadronic phase with its 
gas-liquid phase-transition occupies the small temperature and density/chemical potential 
region. A rich set of solid-state-like phases are expected at large chemical potentials, 
such as color-superconductors, see e.g. 0!^. At high temperatures, a so-called quark- 
gluon-plasma is expected. It is not yet settled whether it is reached by a genuine phase 
transition or a cross-over in the real world. A pure gluon system like the one to be studied 
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\Xo few times nuclear |1 
matter density 



Figure 1.1: The left figure from [Oj sketches the current idea of the phase diagram of QCD 
in the temperature (T)-baryon chemical potential (/i)-plane. It is not yet settled whether 
the transition at /i = is a true phase transition or merely a cross-over, /iq corresponds to 
normal nuclear matter density, x denotes chiral. The right panel shows a more quantitative 
picture with fib = yU. This figure from P| contains the available experimental data. The full 
triangle corresponds to normal nuclear matter and the open squares indicate the position 
of a possible tri-critical point from lattice calculations JU]- The circles are from fits of 
particle ratios with a thermal model [SI E] • They have been obtained by the experiments 
at RHIC at high and low energies, SPS at high and low energies, AGS and SIS in order of 
increasing /i. The dashed and dotted curves show freeze-out trajectories for a hadron gas 
at constant density or energy density, respectively jH]. 



here exhibits a first order phase transition P . 

It is expected that chiral symmetry is restored in the high-temperature phase. This 
is supported by lattice calculations, see e.g. [21 . Concerning the other properties of this 
phase, even the nature of the effective degrees of freedoms is currently under debate. The 
initial picture made was a perturbative gas of free quarks and gluons, leading to the name 
deconfinement. As colored currents are gauge dependent, this view can only be true in 
a figurative sense. However, the idea of almost local colorless objects with otherwise the 
quantum numbers of quarks and gluons as weakly interacting quasi-particles remains. 

The main aim of this work is an analysis of the properties of gluons in the high- 
temperature phase as well as the validity of this simple picture. Still addressing merely 
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equilibrium dynamics, the fate of the confining properties at temperatures above the phase 
transition are investigated. In addition, there is recent evidence that at temperatures 
just above the phase transition the matter is in a non-trivial, strongly interacting phase. 
Whether this extends to all temperatures is also a main subject of this work. 

The basic field theoretical concepts used will be compiled in chapter |2l The Dyson- 
Schwinger equations used to study this theory including the introduction of thermody- 
namic features will be laid out in chapter El A first step will be concerned with inves- 
tigations regarding the infinite temperature limit in chapter 0] This will be followed by 
introducing a high temperature expansion in chapter El The thermodynamic potential, 
analytic properties and other aspects of the solutions will be investigated in chapter (HI 
The results will be interpreted and summed up in chapter |3 Appendix El contains con- 
ventions and commonly used symbols. Five further appendices contain technical details 
which would have broken the line of argument inside the main text. 



Chapter 2 

Aspects of QCD as a Gauge Theory 



2.1 Formulation 

QCD describes the interaction of quarks, massive fermions, through gauge bosons, the 
gluons. The gauge group for the physical QCD is SU(3). In this work, the scope is enlarged 
to an arbitrary semi-simple, compact Lie-group. However, only some of the assumptions 
made later have been tested by lattice calculations and just with respect to SU(Ai'c) gauge 
groups at small Nc- Thus although the results will turn out to be independent of the gauge 
group, the assumptions made are potentially not. 

The classical Lagrangian of such a general gauge theory is given by PP 

= -\F^,F'^^'^ + J2{^'}YzD,^'}-mf^'}^j) (2.1) 
/ 

F;, = d,A:-d,A; + g,r'"^AlA: 
D, = d^-ig^eAl. 

A^j^ denotes the gluon fields with color a and Lorentz index /i. ipj and ipj are the anti- 
quark fields and quark fields, respectively, of flavor / = 1, ...,Nf and color a. For a SU(Ai'c) 
gauge group a = 1, A''^ — 1 and a = 1, ...,Nc. F^^ is the field strength tensor and 
the fundamental covariant derivative, t"' are the generators of the gauge group and /"^'^ its 
structure constants. The m/ are the masses corresponding to the quark flavor /, and gd is 
the gauge coupling where d is the space-time dimension. 7^ denotes the Dirac 7-matrices 

CI- 

The introduction of equilibrium thermodynamics in chapter IHl includes a Wick rotation 
[T!^ and hence the Euclidean version of 1)2.11) will be more important throughout. It is 
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given by ^1] 

L = ^F^^^F^'''^ -J^i^fl'D,^'} -rrif^'}^'}) (2.2) 

/ 

F;, = d,A:-d^A';,-g,r''^AlA: 
= d^ + igdf'AI. 

All further expressions will be in Euclidean space-time, if not otherwise noted. In the 
remainder of this work only the gauge-subsector of ()2.2|) will be investigated and the 
matter content discarded. This reduced theory is known as Yang-Mills theory PI . In 
section 12.31 it will be discussed to which extent this is justified. In section 17.21 a short 
assessment of the influence of quarks on the results presented will be given. 

In chapter |3J the Yang-Mills sector will be coupled to an adjoint scalar field. For the 
sake of completeness, its Lagrangian is given here in Euclidean space as 

£ = lF;,F^'^'" + l(D;V^-0^ + m^0>») + ^0>W (2.3) 
Df = 5'^'d, + g,r'^Al. 

(f)"' is the scalar field of color a, rrih its mass and h its self-coupling. For an SU(A^c) gauge 
group a = 1, A*"^ — 1. -DJ^^ is the adjoint covariant derivative. A three-Higgs coupling is 
not present due to the antisymmetry of the coupling constants and the color symmetry at 
tree-level. 

The next step is to quantize the Lagrangians ()2.2|) and ()2.3|) . 



2.1.1 Quantization, Gauge Fixing and the Gribov Problem 

Quantization can be performed using either canonical or path integral methods. The latter 
will be used here, because they give a somewhat more sophisticated access jTHI. In this 
case, the generating functional is given by 

Z[j;] = I VA^^exp (1 d'^x (-£ + A-^{x)f''^{x))^ =: e^t^^^l, (2.4) 

which depends on the classical sources for the gauge field; the 'free energy' W is the 
generating functional of connected correlation functions. The integration is over the set 
of all configurations of the gluon field AJ^. This set is termed gauge space. The classical 
sources have to be set to at the end of all calculations. 

The prescription ()2.4j] has the problem that not all possible fields contribute inde- 
pendently to the path integral, since Yang-Mills theory is a gauge theory. £ is invariant 
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under local gauge transformations 

a; ^ A;, + D;'6e' (2.5) 

r ^ r + 9dr''^'S9', (2.6) 

where the second transformation is only relevant for ()2.Hj) . 66°' are independent functions. 
Field configurations, which are equivalent up to a gauge transformation ()2.5p or ()2.6p are 
said to lie on the same gauge orbit. This gauge freedom leads to over-counting in the 
path integral ()2.4j) . It is hence necessary to take the gauge freedom into account when 
quantizing. 

The most powerful, albeit technically complicated method to circumvent these prob- 
lems is stochastic quantization [T7j. It is in generally well suited to calculate gauge in- 
variant quantities, including all observables. Furthermore it is possible to obtain gauge 
dependent objects like gluon propagators corresponding to a conventional gauge choice 
[inillHl- For Landau gauge, this process yields the same equations describing the gluons 
[T^ as the approach followed here within the approximation scheme employed. 

The conventional approach followed here is to fix the gauge prior to quantization [T^ . 
The aim is to include only one configuration on each gauge orbit. This is performed by 
introducing an appropriate 5-function in the path integral ()2.4|1 . The argument of this 
function is a local functional of the fields of the form C"[74|^,x], and its equality to 
defines the gauge. The remaining steps are a standard procedure and will not be detailed 
further [l]. 

There is still one point to be mentioned. The gauge fixing usually used is an algebraic 
or differential condition, such as 

d,A; = (2.7) 

for Landau gauge or 

VA" = (2.8) 

for Coulomb gauge. Evidently, ()2.7|1 and ()2.8|1 are not complete gauge-fixings. Any har- 
monic gauge transformation in ()2.7|) and any transformation depending only on time in 
()2.8|) are still allowed. It is possible to fix this residual classical gauge freedom. However, 
the conditions ()2.7p and ()2.8|) are still not unique. In non-abelian gauge theories there 
is more than one gauge-equivalent solution to them, i.e. more than one configuration on 
each gauge orbit satisfies the gauge condition. E.g. in the case of Coulomb gauge ()2.8|1 
a gauge-equivalent solution to the vacuum A" = is a hedgehog configuration [20] • The 
consequence is over-counting in ()2.4j) . This is known as the Gribov problem [20] and it can 
be shown to extend to a large class of local gauges [21] • As the additional copies in general 
involve large gauge field fluctuations, this problem does not appear in perturbation theory. 
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It is not yet known whether there exists any local gauge condition which resolves 
this problem^. With regard to the problems induced by non-local gauge conditions, it is 
worthwhile to investigate if there is a way to circumvent this problem. Indeed it has been 
shown that the zeros of the Faddeev-Popov determinant 



define non-intersecting convex and compact regions of gauge-space. Each region is inter- 
sected at least once by each gauge orbit. The one enclosing the origin and thus perturbation 
theory is called the first Gribov horizon [201 1^- Within it, ()2.9|) is positive. It is possible 
in the approach followed here to ensure M > and thus to be inside the first Gribov hori- 
zon. This will be discussed in section imi However, this condition is not sufficient [23], 
and gauge copies are still present. The copy-free region contained inside the first Gribov 
horizon, called the fundamental region, can be defined using a non-local minimalization 
condition [21]. This in turn implies that there is again no local condition like the Gribov 
horizon to eliminate the copies, and a full solution is still missing. 

Nevertheless, it has been argued that a bounded volume in an infinite dimensional 
space, like gauge space, is dominated by its boundary. Therefore only this region con- 
tributes to objects constructed from a finite number of operators^ [12]. Using this entropy 
argument, the common boundary of the fundamental region and the first Gribov horizon 
is the only region contributing and the Gribov horizon condition is sufficient. For the 
remaining part of this work, this argument will be accepted as an assumption and hence 
the Gribov horizon condition will be used to restrict the space of possible solutions. 

The impact of Gribov copies can be studied by lattice calculations. A relatively weak 
effect on the gluon propagator and a stronger quantitative effect on the ghost propagator, 
to be discussed shortly, is found [23]. It is therefore likely that even if the assumption is 
incorrect, the results found here are still qualitatively reliable. 

2.1.2 Landau Gauge 

Most calculations based on perturbation theory restrict the gauge condition as weakly as 
possible and use the requirement of gauge independence to check the results for errors. 
In a well-defined scheme such as perturbation theory, this is clearly advantageous. In the 

^There are arguments that Landau gauge is still well-defined, if a summation is made over all signed 
Gribov copies |22j . It is however probable that the truncations introduced in section 13.31 will interfere 
with these cancellations. Thus other techniques are necessary here. 

^So the following argument does not necessarily apply to the Polyakov loop or other exponentials of 
operators. 




(2.9) 
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case of non-perturbative calculations no such rigorous scheme exists up to now. Hence, 
the requirement of approximate gauge invariance will actually be used to investigate the 
quality of the solution. To this end, Landau gauge ()2.7|) turns out to be well suited for 
reasons to be described shortly. 

Landau gauge belongs to the class of covariant gauges. Using the standard prescription 
for gauge fixing P], the generating functional ()2.4|) becomes 



where is an arbitrary constant, the gauge constant. Landau gauge is obtained by 
the limit — ^ 0, in which case the last term in ()2.10|) is dismissed. The Faddeev- Popov 
determinant ()2.9j) emerges from an intermediate change of variables as a Jacobian. Faddeev 
and Popov showed [22] that it can be rewritten as a path integral over scalar Grassmann 
fields, leading in Landau gauge to 



Z[j'^] = / VAVcVcexp - d'^x{C + e{x)d^D^'^''\\x) - Al{x)f^''{x)) , (2.11) 



where c°' and c" are the (Faddeev-Popov-) ghost and anti-ghost fields. As they are anti- 
commuting scalars, they may not appear in final states, since they would violate the 
CPT-theorem. This can be guaranteed in both, perturbative calculations and also in 
the non-perturbative calculations in this work, as will be detailed in subsection 12.1.31 and 
section Eini Consequently, no external sources are associated with the ghosts at this level. 
They will be introduced in chapter |Hl when off-shell ghosts will be investigated in more 
detail. Similar to time-like gluons, ghosts contribute indefinite norm states to the Hilbert- 
space already in perturbation theory. This is directly visible from the fact that ghosts act 
as 'negative degrees of freedom' in scattering processes 

The ghost fields have a new global symmetry, the ghost number symmetry. Rescaling 
the ghost fields by a scale transformation exp(s) and its anti-field by exp(— s), leaving all 
other fields unchanged, is a symmetry of the Lagrangian^. It gives rise to the conserved 
ghost number Qg-, in analogy to the fermion number. As the ghosts are the only fields 
carrying them, it is necessary that all observable final states must have ghost number 0. 

Note that the hermiticity assignment of the ghosts in 1)2.111) is different from the com- 
mon choice, which is refiected in the different sign for the ghost term in ()2.11)1 . This would 
lead to problems in general covariant gauges, but in ghost-anti-ghost symmetric gauges 
like Landau gauge this is permissible and for technical reasons advantageous |2Z1 • Of 
these gauges. Landau gauge turns out to be favorable as it is less singular than other 
•^It is a left-over from the original local symmetry, which was broken by gauge fixing. 
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gauges |2ZI. This manifests itself in the non-renormahzation of the ghost-gluon vertex, 
which will be discussed below. A further advantage of Landau gauge is that as long as 
multiplicative renormalizability holds^, it is a fixed point of the gauge parameter, since 
the latter is exactly 0. 



2.1.3 BRST Symmetry 

With the introduction of ghosts in ()2.1H) a further global symmetry arises, which is also 
a residual of the local gauge symmetry. It is named BRST-symmetry after its discoverers 
Becchi, Rouet, Stora, and Tyutin [21]. Consequently it is possible to define the BRST 
charge Qbrst, which has ghost number 1. 

The corresponding symmetry transformations are 

^BRST^^Ka:) = 5\Dfc\x) (2.12) 
5bkstc^{x) = -5\]^gdr''c\x)c%x) (2.13) 

5brstc"(x) = 5\^d,Al{x) (2.14) 
= 5Xg,r'''<P%x)c\x), (2.15) 

where for the Landau gauge the appropriate limit has to be taken. The transformation 
rule for the adjoint field 0" was added for completeness. 6X is an infinitesimal constant 
Grassmann parameter. This defines the BRST-operator s as 



jbrstJ 



6XsF, 



by its action on any field F. As it is defined as the left-derivative of the transformed field 
with respect to 6X, it directly obeys the generalized Leibniz rule. Since it is Grassmann 
in nature, as can be seen from the fact that it changes the number of Grassmann fields in 
the transformations ()2.12II2.1'H|) by 1, this rule reads 

s{FG) = {sF)G T F{sG). (2.16) 

The sign depends on whether F is Grassmann or not. It further carries ghost number 0, 
as the transformation parameter 6X has to carry ghost number -1 since the BRST charge 
carries ghost number 1. As ()2.12|) is a global Grassmann valued gauge transformation, 
the gauge part of ()2.2|) is invariant on its own. The combination of the gauge fixing part 

"^Note that multiplicative renormalizability of Yang-Mills theory with and without matter fields is only 
proven perturbatively order by order. It is unknown if it also holds non-perturbatively. In the results 
presented here it does hold. 
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2.1. Formulation 



and the ghost contribution is also invariant, where in Landau gauge the hmit ^9 — > has 
to be taken. It is possible to linearize the transformation rules (12.12112.1'^ by introducing 
an auxiliary field, the Nakanishi-Lautrup field 5" |2H]- Without going into details, this 
establishes manifestly the nil-potency of the BRST transformation^ 

Sl^ST = 0. (2.17) 

This establishes a closed algebra 

{Qbrst, Qbrst} = 

[iQc, Qbrst] = Qbrst 
[^Qg,Qg] = 

of the residual local gauge symmetry. 

A well-defined nilpotent charge directly splits the state space into three disjoint parts 
[23 IHOl 1^- The states which are not annihilated by the BRST-transformation form 
a subspace Qi, carrying BRST-charge. By acting on these states, daughter states in a 
subspace Q2 are generated which are annihilated by the BRST-charge. The last possibility 
are states which are also annihilated by the BRST-charge but are not generated from 
parent states. These form a subspace Qq. Physical states must be gauge invariant and 
are therefore annihilated by Qbrst 32j. In addition, any states in Q2 do not contribute 
to matrix elements. Therefore the physical subspace is 

Hphys = Ker^BRST/IrnQBRST = Qo- 

It is this subspace in which the perturbatively physical transverse gauge bosons exist, while 
forward polarized gluons and anti-ghosts belong to Qi and backward polarized gluons and 
ghosts belong to Q2. This can be seen directly using the Nakanishi-Lautrup formulation 
of the gauge-fixed Lagrangian [2Hj. Due to the relation of Q2 and Qi, the unphysical 
degrees of freedom are connected by BRST transformations and are thus metric partners. 
They are said to be confined by the quartet mechanism . Hence in perturbation theory 
the physical subspace Qq contains only transverse gluons, and perturbatively unphysical 
degrees of freedom are confined^. One of the confinement mechanisms proposed, the 
Kugo-Ojima scenario discussed in subsection I2.2.2[ requires also transverse gauge bosons 
to belong to either Q2 or Qi and thus provides confinement. 

^Without this field, the nil-potency would only be manifest on-shell PU], 

^In principle it is possible to have states in Qo with non-vanishing ghost number, which would render 
the theory ill-defined (^Jj- This seems not to be the case for Yang-Mills theories. 
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2.1.4 Slavnov- Taylor Identities 

It is possible to directly construct identities relating different Green's functions by calcu- 
lating the BRST transform of an operator expression, using that the vacuum belongs to 
Qq. These are the Slavnov- Taylor identities (STI) jSHlEll- These identities are a result of 
gauge invariance. A failure in fulfilling them therefore indicates a violation of gauge in- 
variance. This property makes them an important technical tool to check the consistency 
of calculations and they will be used in this way extensively in chapters El to El 

In this work two STIs are of particular importance. The first is obtained when forming 
the expectation value of the BRST transform of (fd^A^^ and yields PP after usage of the 
equation of motions 

PuDfM = -^^gP.S''' (2.18) 

with the gluon propagator D^iy. By virtue of Lorentz invariance, the Landau gauge limit 
hence requires D^i, to be of the form 

D,Ap) = P,Ap)^ (2.19) 
P,Ap) = 5,. - (2.20) 

i.e. to be transverse. Z is the associated scalar dressing function. This simple Lorentz 
structure is one of the advantages of Landau gauge. The other prominent feature of 
()2.18|) is its independence of higher n-point Green's functions. In general, the identity for 
an n-point Green's function depends on n'-point Green's functions with n' > n. Hence 
an infinite set of coupled equations arises. Therefore, their use is limited in the case 
of non-perturbative calculations^, as there is not yet any possibility to assess a-priori the 
contribution of the n'-point Green's functions. In special kinematic regions, however, those 
unknown contributions may drop out, providing relations between n-point functions only. 

A further result which can be obtained from ()2.18|) in connection with the equation of 
motion of the ghost is that the ghost-gluon vertex is undressed for a vanishing incoming 
ghost momentum in Landau gauge |3S1 EEj 

]i rf^Mc^p, q, -p~q) = ig.r'^qu. (2.21) 

It is thus not divergent and its renormalization constant Zi can and will be set to 1 here. 
This is maybe the most important property of Landau gauge. It makes it much more 



^In perturbative calculations such contributions can be neglected as they are of higher order in the 
expansion parameter. 
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2.2. Confinement 



tractable than other gauges j2Zj- The ghost-gluon vertex will be further investigated in 
subsection 12 .2. HI and used in chapters 01 to El 

The second identity is concerned with this ghost-gluon vertex. It can be obtained from 
the BRST-transform of c°'c^c'^ jHSl IHZI and reads 



ik,T-^^\p, -q)G{q') + zg,r^^(p, -k),G{k^) = p'Z^ ' ^iJ^ + V- (2.22) 



Tj summarizes contributions from n'-point Green's functions with n' > 3. Color indices 
have been removed by assuming a tree-level^ color structure. This assumption is discussed 
in more detail in chapter |21 G is the ghost dressing function, defined via its propagator as 



D^^{p^) = (2.23) 



The identity ()2.22|1 is consistent with the bareness of the ghost-gluon vertex 1)2.211) if the 
contributions from the higher Green's functions vanish in this limit. 



2.2 Confinement 

As one of the main observables for this work is the absence or presence of confinement 
it is necessary to detect it. Two fundamentally different approaches to this question will 
be discussed. The first is concerned with the mere statement of confinement. This will 
be investigated in subsection 12.2.11 The other approach consists of criteria deduced from 
possible confinement mechanisms. Three of them will be discussed in sections 12.2.21 to 
12.2.41 The number of proposed mechanisms for confinement is large. Only those which 
generate criteria that can be tested using the objects obtained in this work will be discussed 
here. A more general overview can be found in [T^ . 

In any case, cluster decomposition [HHllSni must be violated for colored objects to allow 
for a long range force. This necessarily implies the existence of a massless excitation in 
the complete state space. On the other hand since all experimental results show validity 
of cluster decomposition, the physical mass spectrum of QCD must have a mass gap for 
colorless objects. Otherwise it would be possible to scatter a colorless object into far apart 
colored objects, the so-called "behind-the-moon" problem. However, the existence of a 
massless excitation alone does not suffice for confinement of color, as it is necessary to 
show that it is not part of the physical spectrum. 

^This is exact in perturbation theory. 
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2.2.1 Criteria for Confinement 

There are two criteria which will be used here. One implies the other. The basic object 
is in both cases the spectral density p jlO] of a given particle. For a particle to exist as a 
physical final state having a Kallen-Lehmann representation, it has to have a positive semi- 
definite spectral function^. This is known as the Osterwalder-Schrader axiom of refiection 
positivity [SHI EI]- A violation of positivity can be tested using the Schwinger function to 
be discussed in chapter El 

A stronger condition is violation of the Oehme-Zimmermann super-convergence relation 
02] • The spectral representation of a particle of mass itiq is pUj 

oo 

with m > mo being the threshold mass for contributions above the single-particle pole, 
the multi-particle threshold. Z is the positive overlap between an in-state and the field 
described by the propagator D. If the propagator vanishes at 0, 

\imp'^D{p) = 0, (2.25) 

then the spectral function must be at least partly negative, thus implying the first condition 
above. Secondly, a vanishing propagator aX p = implies the absence of a Kallen-Lehmann 
representation. The particle can thus be not a physical particle anymore: The particle is 
confined. This gives the second condition for confinement. 

For massless particles the interpretation of ()2.25p is more direct. As = is just the 
statement of the particle being on-shell, it implies the vanishing of the on-shell propagator: 
The particle does not propagate and is thus confined. 

Note here a subtle difference. It is possible to think of confined particles as either 
confined due to dynamic effects or due to being unphysical. For example magnetic con- 
finement of neutrons in a magnetic field is such a case: The particles are physical but they 
are bound in a way disallowing separation. The confinement e.g. of ghosts differs from 
this. These are unphysical particles and are thus confined in a different sense. Note that 
in an unbroken non-abelian gauge theory, charged currents are gauge-dependent and can 
thus not be observed directly. Observing a quark thus corresponds e.g. to an observation 
of a colorless object with non-integer electric charge. 

^In a theory with indefinite metric, as Yang-Mills theory, for unstable particles such that the width 
exceeds the mass the spectral function is also not necessarily positive semi-definite. They do not occur in 
final states, since they decay when letting t ^ oo. 
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One of the major questions concerning confinement is which of both possibihties apphes 
to colored objects. The Kugo-Ojima criterion described in the next subsection favors the 
latter option, as colored objects will turn out to be automatically BRST-charged. The 
more common attitude is the first option. 

2.2.2 Kugo-Ojima Scenario 

The Kugo-Ojima confinement scenario [Hi] puts forward the idea that all colored objects 
form BRST-quartets and therefore do not belong to the physical state space. The metric 
partners of transverse gluons would be ghost-gluon bound states. 

Thus, the confinement mechanism is essentially the same as in the case of ghosts in 
perturbation theory. This scenario is based on a rigorous derivation and requires three 
preconditions. One of them is an unbroken BRST charge also in the non-perturbative 
regime. Whether this is the case is not known^°. The second is the failure of the cluster 
decomposition theorem and hence the existence of a massless excitation. This is also 
unknown and can not yet be attacked in the approach used here. The third ingredient 
is an unbroken global color charge. In Landau gauge, this condition can be put into the 
form 03] 

lim p^Dg{p^) oo, (2.26) 

where Dq is the propagator of the Faddeev- Popov ghost ()2.23|) . This scenario also neces- 
sarily implies that ()2.25|1 holds for the gluon. Both of these conditions ()2.25|1 and ()2.26|1 
will be checked. 

Note that this scenario does not imply that colored objects do not have an asymptotic 
field. However, as they are BRST charged a configuration with arbitrarily many colored 
objects belongs to one equivalence class, and only the colorless objects contribute to matrix 
elements. 

2.2.3 Zwanziger-Gribov Scenario 

The central idea of the Zwanziger-Gribov scenario fOl I2II1 IMl 1^ is that zero-modes at 
the common boundary of the first Gribov horizon and the fundamental region dominate 
the infrared properties and thus generate confinement. Both regions necessarily have a 
common boundary. This boundary is convex, compact, and includes the origin 

As gauge space is infinite-dimensional, all the volume will be concentrated at the 
boundary, and the system is dominated by it. Since the Faddeev-Popov-determinant ()2.9|) 

^°It is even unknown how to define a non-perturbative BRST charge. Due to the non-trivial topology 
of the gauge group this necessarily has to be done patch-wise in gauge space. 
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vanishes there, the corresponding excitations have to be long-range. Using stochastic 
quantization and performing a Landau gauge hmit under certain assumptions, it can be 
shown that ()2.26|) follows. In addition, it follows that the infrared limit is dominated by the 
ghost-term of 1)2.111) alone P^. As this term can be written as a BRST-transform, Yang- 
Mills theory in this limit is a topological field theory of Schwarz type with no propagating 
modes pUj. Thus colored objects do not appear, implying confinement. Using these results, 
()2.25|) for the gluon is also obtained, thus leading to the same criteria as the Kugo-Ojima 
scenario. This analogy, if it is more than mere coincidence, is not yet understood. 

It is also intuitively clear that a strongly divergent ghost propagator at zero momentum 
can mediate confinement. After Fourier-transformation such an infrared divergence relates 
to long-ranged spatial correlations. These are stronger than the ones induced by a Coulomb 
force since the divergence in momentum space is stronger than that of a massless particle. 

A bare ghost-gluon vertex in the infrared is sufficient to generate this behavior. Such 
a vertex would also be consistent with the perturbative renormalization group J^. This 
Zwanziger-hypothesis has been checked numerically in context with this work jUj and 
turns out to be well fulfilled. 

Recently, investigations in Coulomb gauge indicate, that a similar connection between 
the dynamics on the Gribov horizon and confinement of gluons holds there as well |1H] ■ 

2.2.4 Gribov-Stingl Scenario 

The last scenario investigates a manifestation of confinement different from the two pre- 
vious ones. The Gribov-Stingl scenario [201 EHl HD] puts forward the idea that confined 
particles have one or more pairs of complex conjugate poles, and thus cannot be physical 
states. In general, such a pole structure leads to violation of causality |2H]- However, for a 
special structure of the propagators it can be shown that it is possible to reconcile such a 
pole structure with causality on the level of the ^-matrix [IHl- The essential argument is 
that due to the absence of real poles, application of the LSZ-reduction formula [2H] leads 
to vanishing matrix elements for states with colored objects [IHI. Thus colored objects 
can only exist as short-hved quantum excitations. In addition, colorless objects only have 
real poles but no continuum, if they cannot decay into colorless objects 

2.3 Interplay with Quarks 

As quarks will be neglected throughout this work, a short comment to which extent the 
results may be affected by the presence of quarks is in order. Calculations in the vacuum 
show that the presence of quarks does not qualitatively alter the gauge propagators, as 
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long as there are less than 5 light flavors [SI]. This is not expected to change at flnite 
temperature, as quarks acquire an effective mass increasing with temperature and thus 
behave more perturbatively [S^. Especially they should not contribute to the inflnite 
temperature limit of the gluon propagator. Thus the results presented will probably not 
be changed qualitatively by quarks; this is only a conjecture, though. Especially in the 
vicinity of the phase transition quarks are relevant. Lattice calculations indicate a possible 
change of the order of the phase transition when including quarks [H] , but this issue is still 
under debate. 

The impact on quarks by the results found here may be more signiflcant but also more 
intricate to understand. Lattice calculations show a drastic change in the properties of 
quarks above the phase transition. The free energy of two static quarks changes sub- 
stantially and may change from a linear rise with distance to a logarithmic or even flat 
shape |53j. However, as the mechanism of quark conflnement is not understood yet ;54j . 
assessing the consequences of this change is not simple. The most prominent change is the 
restoration of chiral symmetry [Oj, which surprisingly occurs at the same temperature as 
the changes in the gauge sector. 

One possible chain of arguments to understand this coincidence is based on evidence 
from lattice calculations that the underlying degrees of freedom confining quarks and 
gluons are topological objects [SS]. These carry topological charge, and can thus be brought 
into connection with zero modes of the Dirac-operator by the Atiyah-Singer index theorem 
|56j . The density of such modes again is related to the chiral condensate by the Banks- 
Casher formula ^7\. If this chain of arguments is correct, it gives a connection between 
confinement and chiral symmetry breaking. As will turn out in this work, at least part of 
the gluons are confined even at high temperatures. It is therefore still unclear what the 
relation above the critical temperature is and how the quark properties are affected. 



Chapter 3 

Derivation of the Dyson-Schwinger 
Equations 



Knowledge of all of the Green's functions would grant complete knowledge of a theory 
The inverse 2-point Green's functions are the propagators. These are the main 
objects of interest here, as they carry information concerning confinement and other non- 
perturbative properties. The equations determining these Green's functions are the Dyson- 
Schwinger equations [SHj (DSEs), which can be obtained using the functional equations of 
motion. This will be described in section 13.11 Since the aim of this work is to investigate 
the equilibrium high temperature phase, temperature will be introduced into the DSEs 
in section 13.21 As there are an infinite number of coupled DSEs, it is generally not 
possible to solve these simultaneously. To obtain approximate solutions, truncations are 
necessary and these are discussed in section 13.31 The connections to perturbation theory 
and renormalization are investigated in sections 13.41 and 13.51 

This work is based on a calculational scheme which has been applied successfully in 
the vacuum to both Yang-Mills theory and full QCD, as will be described in section 

3.1 Vacuum Formulation 

The most straightforward way to derive the DSEs for a generic field 0" is by using the fact 
that the integral of a total derivative vanishes [Hj 
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3.1. Vacuum Formulation 




Figure 3.1: The Dyson- Schwinger equations for the 2-point functions described by the 
Lagrangian ()2.3|1 . Dotted hues are ghosts, wiggly hues are gluons and dashed hues are 
Higgs. Lines with a dot are full propagators, vertices with a small black dot are bare 
vertices and white circles are full 3- and 4-point functions. 

5* is the action, j"" is the source of 0° and the integral is over full field space. Performing 
the derivative and pulling the resultant factor out of the integral by replacing 0" with 
S/Sj"", the prescription to calculate the full one-point Green's function is obtained as 

Further derivatives with respect to the fields generate the Green's functions of arbitrarily 
high order. They form an infinite set of coupled non-linear integral equations. For Yang- 
Mills theory and full QCD, these are known for the propagators, see e.g. ^3]. For the 
Lagrangian ()2.3p these are derived for all 2-point Green's functions in appendix El A 
graphical representation of these is given in figure 13.11 

In the case of Yang-Mills theory or QCD, it is not necessary to integrate over all field- 
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space for Green's functions. As the Faddeev-Popov determinant ()2.9|1 appears explicitly in 
()2.1()|1 . it is sufficient to integrate only over the first Gribov horizon, since the determinant 
vanishes on its boundary |19j. Hence the DSEs have the same form whether the integration 
is over all of gauge-space or only over the first Gribov horizon. Therefore it is necessary 
to restrict the space of solutions to those from inside the first Gribov horizon. This is 
discussed in section IH.HI 



3.2 Dyson-Schwinger Equations at Finite Tempera- 
ture 

The starting point of the analysis are the DSEs of Yang-Mills theory described by ()2.2|) 
with the quark fields set to zero. Besides neglecting all equations for n-point Green's 
functions with > 2, all genuine full two-loop graphs in figure ITT] are neglected, too. The 
consequences of this and further truncations will be discussed in detail in section l!01 The 
remaining system is then represented in figure IH.2I and given by 

+ / -^fT^^'^'^-^.P. ^ - P)D%{P - q)D'jiq)Tf^''^fi-p, q,p- g)(3.2) 
D:t'\p) = 5''\5,y-P,Pu)+T^.'' 

_g_rti,c.A.ca^_^ - g, q,p)D^^{q)D'j{p + g)rf ^'^^"(-g, p + q, -p) 

+ \ I -0f%f'''\PA-p.-q)Dll{q)D%{p-q)V^^^^^^^ 

(3.3) 

Here T^'^^ is the full ghost-gluon vertex and F^^ is the full three-gluon vertex. The respec- 
tive tree-level quantities are denoted by a superscript 'tl' and given in ()B.14|) and ()B.15|) . 
T^^ is the tadpole term. The full vertices will be discussed further in section 13.31 Note 
that these are in general Minkowski-space equations, but equations ()3.2|1 and ()3.3|1 have 
already been rotated to Euclidean space. 

To obtain the equilibrium Green's functions at a temperature T, the Matsubara or 
imaginary time formalism is implemented EH] • This amounts to a compactification 
of time and entails a genuine Euclidean formulation. All objects depend on the three- 
momenta p and the Fourier-component po separately. Both ghosts and gluons have to 
obey periodic boundary conditions jHO], thus 

po = 27mT, neZ. (3.4) 
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3.2. Dyson-Schwinger Equations at Finite Temperature 




Figure 3.2: The truncated Dyson-Schwinger equations for the 2-point functions of Yang- 
Mills theory at T = 0. Dotted lines are ghosts and wiggly lines are gluons. Lines with a 
dot are full propagators, vertices with a small black dot are bare vertices and white circles 
are full vertices. 

Hence in both soft mode exists, defined by Euclidean p'^ <^ T, i.e. those with 

Po = in contrast to the hard modes with pq ^ 0. Note that although Lorentz invariance 
is no longer manifest in this formulation, it is not lost [HI]. Furthermore, as the gluon 
is a vector-particle, its propagator exhibits two independent tensor-structures and two 
independent dressing functions^ T!^ instead of one as in ()2.2()|1 . 

DM = PtAv)^^^ + PlAv)^^^^ (3.5) 
PtAp) = S,.-^^ + S,o^ + 6o.^-6,oSoJl + 4) (3.6) 

pZ pZ pZ \ P J 

PlAp) = PAp)~PtM- (3-7) 

The projectors Pt^iv and Pl^v are transverse and longitudinal with respect to the heat bath 
or alternatively with respect to the three spatial dimensions. Both are four-dimensionally 
transverse, so ()3.5|) still satisfies the STI ()2.18|) . This is necessary, as STIs are still valid 
at finite temperature [Sn]- For T = 0, Z = H and ()2.2()|1 is recovered. Note that ()3.6|1 
projects solely on the space-space components and the pq = 0-component of ()3.7p solely on 
the 00-component of the propagator. Hence the soft mode of Z is purely chromomagnetic, 
while the one of H is purely chromoelectric. For the ghost, being a scalar, one dressing 
function is sufficient at finite temperature. 

In general, all full Green's functions obtain a much more complicated tensor structure 
at finite temperature, e.g. the ghost-gluon vertex obtains four instead of two tensor struc- 
tures. This generates a significant amount of technical problems. As the main assumption 
of the truncation scheme presented in section is the negligibility of dressings other than 

^In general, three dressing functions exist, but only two are independent due to the STI H2.18|l . 
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Figure 3.3: The truncated Dyson-Schwinger equations for the 2-point functions of Yang- 
Mills theory at T 7^ 0. Dotted lines are ghosts, wiggly lines are 3d-transverse gluons and 
dashed hues are 3d-longitudinal gluons. Lines with a dot are full propagators, vertices 
with a small black dot are bare vertices and white circles are full vertices. 



those of propagators, the irrelevance of such effects is assumed as well. Hence the tree-level 
tensor structure is used for all full vertices. The last issue to be addressed before writing 
down the equations is the one determining the gluon. It is a matrix equation, but only 
two components are independent. The most direct way to obtain two scalar equations for 
the two dressing functions is to contract the gluon equation once with ()3.6|1 and once with 
fl3.7|) . The more convenient way is to contract with 

P^u = CPt,u + (C - 1)At,. (3.8) 
pLu = ^Pl,u + - 1)Al,u. (3.9) 

Varying the parameters ( and ^ allows to investigate the amount of gauge invariance 
violation, as discussed in section 13.31 The tensors ()3.8j) and fl3.9|l reduce to ()3.fij) and ()3.7|1 
at C = ^ = 1- The selection of At/l^iu is governed by practical aspects and will be different 
for the high temperature case and the finite temperature case. Thus writing down the final 
equations will have to await chapters IH and respectively. Graphically the equations for 
the three scalar functions G, Z and H at finite temperature are given in figure 13.31 



3.3 Truncations and Constraints 

As already indicated, the DSEs to be solved in the following chapters are truncated. 
Although in the following several arguments will be made to support the truncations, there 
is no method (yet) known which permits an a priori controlled truncation. In addition, no 
possibility is yet known to conserve local internal symmetries at least on the level of the 
truncation. Even for global symmetries this leads to enormous technical complications, see 
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e.g. [021. Furthermore any such truncation necessarily violates unitarity. Also no scheme 
is yet known which at the same time conserves energy and agrees with perturbation theory 
in the far ultraviolet. Hence it is currently only possible to assume that all these effects 
do not contribute significantly. This can only be done by comparing a posteriori either to 
experiments or different methods. As the objects treated here are gauge- variant quantities, 
only other calculational methods are available for comparison. Only lattice calculations 
have been performed in the region of interest up to now and only for SU(A''c) gauge theories. 
It will turn out that the results agree remarkably well, considering the drastic assumptions 
made. By systematically assessing the error due to gauge invariance violation, it is found 
that the effects of the truncation are of a quantitative nature only. Hence the truncation 
scheme to be described seems to be applicable. 

The most extreme truncation is the neglect of the equations of all ri-point Green's 
functions with n > 2. This truncation is well justified in the ultraviolet due to asymptotic 
freedom. In the infrared, assuming the correctness of the Zwanziger-Gribov scenario, 
this is also justified and the results support this assumption self-consistently. At mid- 
momenta, significant deviations are to be expected. Indeed this is the region where the 
largest deviation from lattice results will be found. The further truncation is to neglect all 
genuine full 2-loop contributions within the equations for the 2-point Green's functions. 
Largely the same arguments apply here: The 2-loop contributions are sub-dominant in 
the ultraviolet and most probably sub-leading in the infrared. The latter was checked and 
found to be correct for tree- level instead of full vertices ^|. In addition it was found that 
only considerable fine-tuning of the vertices allowed the 2-loop contributions to become as 
leading as the 1-loop contributions in the infrared . 

The next assumption is that the color structure is the same as in perturbation theory. 
All investigations concerning this point have supported this assumption [Hll . This also au- 
tomatically entails the vanishing of all 1-point Green's functions due to the antisymmetry 
of the color structure of the vertices. 

The last ingredient of the truncation is the construction of the remaining full ver- 
tices. The ghost-gluon- vertex will be kept at its tree- level form ()B.14|1 . motivated by the 
Zwanziger-Gribov scenario. This is exact for vanishing incoming ghost momenta as dis- 
cussed in subsection 12.1.41 A bare vertex is also supported by numerical studies in three 
and four dimensions jlZj and lattice calculations in four dimensions jSHI. Furthermore, 
at least in the vacuum, the qualitative nature of the infrared solution is independent of 
the detailed structure of the ghost-gluon vertex to a large extent [SS]. It has also been 
shown that, under weak assumptions, the qualitative infrared solution for the ghost is 
independent of the truncation . 



Chapter 3. Derivation of the Dyson-Schwinger Equations 



25 



The various three-gluon vertices for 3d-longitudinal and 3d-transverse gluons are not 
fixed yet. These will be constructed by requiring minimal gauge invariance violation. This 
will be addressed in chapters |3] and completing the truncation scheme. 

Concerning the artifacts of the truncation, it is not possible to solve the problem of 
unitarity violation. Unitarity will always be violated as long as the full infinite system is 
not solved^. The problem of energy conservation will be addressed in section in21 

The last point to be discussed is the violation of gauge invariance. There are two 
aspects to be treated. The first is the problem of Gribov copies. Due to the arguments 
given in section 12.1.11 this can be resolved by requiring 

G{p^)>0, Z{p')>0, H{p')>0, (3.10) 

since this guarantees to stay within the first Gribov horizon. Note that by condition ()3.10p . 
the ghost propagator is negative definite and cannot have a positive semidefinite spectral 
function. 

The second problem is much harder to address. Even provided the Gribov problem is 
solved, gauge invariance is violated, as the STIs are no longer fulfilled. Indeed it is not 
even possible to test whether the STIs are fulfilled, as in general the STIs for an n-point 
Green's function contain contributions from n'-point Green's functions with n' > n. These 
are negligible in perturbation theory because they are of higher order in the expansion 
parameter, but this is not true in general in non-perturbative calculations. In principle it 
would be possible to test the STIs when truncating them to the same level as the Green's 
functions, but it turns out that this leads to inconsistencies, see e.g. Nonetheless 
trying to construct vertices which as best as possible fulfill the truncated STIs, it is found 
that the results are only weakly affected compared to tree- level vertices [33]. This gives 
confidence that these violations are small. 

There are two consequences of these gauge invariance violations. The first is the 
appearance of spurious divergences in cases where the degree of divergence is lowered by 
gauge invariance compared to naive power counting. This is the case for the gluon self- 
energy. These spurious divergences have to be removed. This will be done using the tadpole 
terms, which are not left as free parts of the equations, but are chosen to compensate such 
spurious divergences and other artifacts of the truncation, thus mimicking their role in 
perturbation theory. 

The second consequence is the violation due to finite contributions of the STIs. In the 

remainder of the work, the main equation to test the amount of gauge invariance violation 

^In perturbation theory, these violations are of higher order in the expansion parameter and can thus 
be neglected. The same applies to the results here in the realm of applicability of perturbation theory, as 
the same results are obtained in this domain. 
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will be the STI for the gluon propagator, ()2.18j) . It is for this reason the projectors ()<18|1 
and are used instead of and ()H.7|1 . If ()2.18|1 was exactly fulfilled, the results 
obtained would be independent of C and C,- If the results do only weakly depend on these 
parameters, gauge invariance violations are most likely small [HHl IHH HO] ■ The variational 
range for ( and ^ in the case of a violation cannot be extremely large, since otherwise the 
projection will be primarily on the gauge-violating longitudinal part and will not give rise 
to further useful information. 

3.4 Perturbation Theory 

The DSE approach followed here aims at the full 2-point Green's functions. Hence it is 
necessary that perturbative results are embedded in the final results. For sufficiently large 
momenta the Green's functions must reduce to their perturbative counterparts, due to 
the asymptotic freedom of Yang-Mills theories. As the DSEs are truncated at one-loop 
level, this reduction can only be correct to leading order (LO) in g'^/p^~'^-, where p is 
the momentum scale. Subleading contributions will necessarily deviate from perturbation 
theory. This also guarantees that the violation of gauge invariance will be not worse than 
in LO perturbation theory, and thus at least the same level of gauge invariance is achieved. 

The high temperature limit in chapter |3] will give an explicit example of this. In case 
of the finite temperature corrections in chapter the results are restricted to momenta of 
the order of T, see section 15.51 Therefore, perturbation theory will only be reproduced if 
these momenta are already in the perturbative regime. 

In the case of 4d vacuum calculations, agreement with LO resummed perturbation 
theory turns out to be a significant task and requires modifications of the three-gluon 
vertex EH] • Comparison to perturbation theory is at the current level of truncation 
the only possibility to compare to experiment, and thus an important constraint. 

3.5 Renormalization 

As in perturbation theory, the usual ultraviolet divergences of Yang-Mills theory are en- 
countered and must be regularized and renormalized jTj. A wide variety of possibilities 
to deal with the divergences at the perturbative level exist, especially with dimensional 
regularization for gauge theories. However, most of these concepts, including the latter, 
are not applicable to non-perturbative calculations [7TJ[72j. 

Due to the lack of a symmetry conserving regularization scheme for DSEs which is of 
technically acceptable complexity, an alternative route is chosen. It is always possible to 
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use gauge non-invariant regularization prescriptions if appropriate compensating counter- 
terms are chosen [21]. This approach will be employed in chapter El There the effects of 
a non-gauge-invariant regularization will be absorbed into the tadpoles. 

The remaining divergencies are then regularized and renormalized using counter-terms. 
In general a counter-term for e.g. renormalizing the wave function is introduced in a 
Lagrangian by performing the replacement 

L = (9^0)2 ^Cr = (9^0)2 + 6Z{d^<p') =: Z{d^ct>f. 

Here cf) is an arbitrary field and the counter-term 5Z has been chosen such as to cancel 
the divergences jHj. This guarantees multiplicative renormalizability even in the non- 
perturbative regime by explicitly constructing the wave-function renormalization constant 
Z. At finite temperature no new divergences arise compared to the vacuum divergence 
structure of the phase the theory is in [20]. Nevertheless, the finite parts of the counter- 
terms may depend on temperature. 

When treating finite temperature effects in chapter El a few subtleties concerning renor- 
malization in the employed truncation scheme arise. These are treated in section 15.5.21 

3.6 Solutions in the Vacuum 

Solutions of the DSEs in the vacuum ^1 ESI EE] have already convincingly demonstrated 
the non-triviality of the infrared regime and also showed good agreement with lattice 
calculations. This supports the applicability of the truncation scheme. In this section 
these results will be described briefly to put this work in the appropriate context and 
demonstrate that the method is sufficiently stable for an extension to finite temperature. 

3.6.1 Yang-Mills Theory 

For the pure Yang-Mills sector, results have been obtained with increasing precision over 
time^. The DSE results satisfy ()2.25|) and ()2.26|) and thus exhibit manifest gluon con- 
finement, in accordance with the Kugo-Ojima and Zwanziger-Gribov scenario ^1 ESI EH]. 
Also, the analytical structure has been understood to some extent [54J- The results jHE] 
are shown in figure EiH compared to lattice results [7^ ITUj . 

The lattice points farthest in the infrared suffer from finite volume effects and are 
expected to bend down for larger lattice volumes. In 3d-calculations, where significantly 

■^The earliest attempt neglected the ghost contribution |[74||. This "Mandelstam approximation" leads 
to results contradicting recent lattice results and is therefore dismissed today. 
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Figure 3.4: The left panel shows the ghost dressing function and the right panel the gluon 
propagator in vacuum Yang-Mills theory [BH] • The results are compared to lattice data of 
the ghost dressing function [73] and the gluon propagator TU'. The indicated errors are 
statistical only. 

larger lattices can be used, this is indeed the case (2^ and will be seen when comparing 
to lattice results in section 14.61 Recently, the results have also been confirmed by exact 
renormalization group methods |78j. These results give confidence that the method can 
be apphed to the finite temperature well. 

A further result of the vacuum studies is that the quantity 

a{fi) = a{s)G{s, ix^Z^s, fi), (3-11) 

where /i is the renormalization scale and s the subtraction point, is a renormalization group 
invariant, and agrees with the running coupling in the perturbative regime. Although it 
is under debate what the non-perturbative definition of a coupling constant is, if any, it 
is possible to investigate this quantity. It does not exhibit a Landau pole and has an 
infrared fixed point of a{0)Nc = 8.915 jHE]- This result has also been confirmed by exact 
renormalization group calculations [175] . 



3.6.2 Full QCD 



DSEs have also been applied to full QCD in two different ways. A more phenomenological 
ansatz has employed model gluon propagators. This has led to extensive and successful 
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Wave function renormalization | Mass function | 




Figure 3.5: The quark propagator in the chiral hmit fET\ compared to lattice results using 
the asqtad action ^j. The left panel shows the quark wave-function renormalization 
compared to lattice results. The right panel shows the quark mass function compared to 
chirally extrapolated lattice results. The indicated errors are statistical only. 

investigations of hadron phenomenology. For a review see e.g. ^UZHj- As the approach 
followed here is more bottom-up, this aspect will not be discussed further. The other 
approach couples the Yang-Mills sector described above to the quarks and solves the 
corresponding DSEs in a similar truncation scheme [SI]. The general structure of the 
Euclidean quark propagator is 

with the wave-function renormalization Z and the mass function M. For less than 4 light 
or massless quarks (as it is realized in nature), the Yang-Mills solutions remain essentially 
the same, besides the according changes in the ultraviolet anomalous dimensions. Thus 
even the effect of massless quarks is small. The two independent dressing functions Z and 
M are shown in figure EUl Chiral symmetry breaking is manifest. However, the amount 
of symmetry breaking is sensitive to the specific quark-gluon vertex, as the quark results 
are in general. In the present case a modified Curtis- Pennington vertex has been employed 
[SI]. Also, the problem of quark confinement is not yet understood in this ansatz, as the 
information obtainable from the quark propagator alone are not conclusive 



Chapter 4 

Infinite- Temperature Limit 



The infinite-temperature limit of Yang-Mills theory is presented in this chapter^. Although 
the limit itself is purely academic, it is valuable not only due to technical simplifications. 
Firstly, properties surviving in this limit will also be present at lower temperatures. Es- 
pecially persisting non-trivial effects are genuine features of the high-temperature phase. 
Secondly, it is found in the finite-temperature calculations of chapter El as well as in lattice 
calculations [HSj that the propagators are already close to their asymptotic values for quite 
small temperatures, as low few times Tc- 

The chapter starts with the derivation of the infinite-temperature limit in section 14.11 
The emerging theory is a 3d- Yang-Mills theory with an additional adjoint Higgs field. The 
infrared and ultraviolet properties of this theory will be inspected in section ^21 It is then 
studied numerically in three truncation schemes. The first is the ghost-loop-only scheme, 
to test the Zwanziger-Gribov scenario, in section 14.31 The next schemes address the pure 
Yang-Mills theory in section 14.41 and finally the full theory in section 14.51 The numerical 
method deployed is described in appendix |Fl 

A comparison to lattice results will be made in section 14.61 The 3d-theory is not 
only relevant to the high temperature behavior of Yang-Mills theory. As an example, a 
recently obtained relation between Landau gauge in 3d and Coulomb gauge in 4d [HHI will 
be addressed in section 14.71 

As this chapter deals nearly exclusively with a 3d-theory, the notation p = |p| is used, 
if not noted otherwise. 

4.1 From 4d to 3d 

To obtain the infinite-temperature limit, temperature is introduced into the vacuum equa- 
tions ()3.2j) and ()3.3|1 as detailed in the previous chapter. To obtain scalar equations for the 
^Most of the results presented are published in ^T]. 
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dressing functions Z and if, the additional tensor structures in ()H.8|1 and are chosen 
most conveniently in 4d as 

At^. = - (l + ^) <^mo^o.) (4.1) 

Al^^u = (^l + ^^6^oSo.. (4.2) 

In [HI] it was demonstrated that, in zeroth order, the infinite-temperature limit can be 
found by neglecting any contributions from Matsubara frequencies different from zero. 
This yields an effective 3d- Yang-Mills theory with an additional adjoint Higgs. The Higgs 
field is the Aq field of the 4d-theory, and therefore the number of degrees of freedom is 
conserved in this process. 

By considering the structure of the projectors ()3.8j) and ()3.9|1 in the case of po = 0, 
which is the zeroth Matsubara frequency, it is also possible to find the connection of the 
3d and 4d degrees of freedom. At po = 0, -P|^^ becomes ^oo; independent of ^. It projects 
out the time-time component of the propagator, which belongs to the Aq field. Therefore 
the 3d-longitudinal part of the 4d gluon propagator corresponds to the Higgs propagator. 
On the other hand, P^^^, becomes 

4 = *.-C^ (4-3) 

with zero time-time and time-space components. At = 3 this is the Brown-Pennington 
projector of the 3d-theory jHSj. f)4.3|) projects onto the 3d-subspace and thus establishes 
the connection between the 3d-transverse gluon and the gluon of the 3d-theory. Therefore 
the dressing functions H and Z of the gluon propagator ()3.5j) describe the Higgs and the 
3d gluon, respectively. 

This amounts to integrating out the hard modes at tree-level. In general, this is not 
sufficient jHE]. Higher order effects of the hard modes can potentially still influence the 
interactions of the soft modes. Therefore the general prescription to obtain the effective 
3d-theory is to write down the most general Lagrangian allowed and match the parameters 
by comparing to the 4d-theory This can be done e.g. by lattice calculations |H21 and 
perturbation theory (HH] . In the present case a tree- level mass for the Higgs and a 4-Higgs 
coupling are additionally present. This also modifles the Dyson- Schwinger equations, and 
they therefore have to be rederived. As an explicit example of this process the generation 
of the tree- level mass will be demonstrated in section EISl 

Hence the Lagrangian ()2.3|) governs the 3d-theory and describes a Yang-Mills fleld 
coupled to an adjoint scalar fleld ^^ISSj. AH occurring constants are effective constants, 
which arise by integrating out the hard modes. The effective constants can only be derived 
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by calculating the full theory. Therefore the values obtained from lattice and perturbative 
calculations will be used here, as listed in |H2.j. In most results these constants are irrelevant 
and will drop out, except for the objects discussed in chapter IHl 

The only exception is the 4-Higgs coupling constant h, which can be uniquely deter- 
mined already in the 3d-theory. In principle, the Higgs self-energy may contain linear 
divergences, if this theory stood on its own. However, the Higgs field is only a component 
of the 4d gluon field, and thus should not contain a linear divergence due to the STIs of 
the 4d-theory. Implementation of this requirement, in leading-order perturbation theory, 
fixes h as 



1X1 Note that by ()4.4j] exact t'Hooft scaling [87^ is also maintained, which would otherwise 
be broken by the Higgs-tadpoles. 

Hence the 3d-theory is finite and therefore all renormalization constants can be set to 
1, i.e. no counter-terms are necessary. The divergences have not disappeared, though. 
By sending the temperature to infinity while maintaining for the renormalization scale 
fi ^ T, renormalization takes place at p — oo and can therefore be neglected at finite 
momenta. This intuitive argument is shown to be correct in chapter where the limit is 
taken explicitly. 

The DSEs for the Yang-Mills sector are already known, see e.g. I79j . Adding the 
Higgs, the ghost equation will not be modified compared to pure Yang-Mills theory, since 
no tree-level Higgs-ghost- coupling is present. The remaining alterations of the equations 
due to the Higgs are derived in appendix iBl They lead to the DSEs 
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Figure 4.1: The system under study in the infinite-temperature hmit. Wiggly fines rep- 
resent 3d-transverse gluons, dasfied fines Higgs and dotted fines gfiosts. Lines witfi a dot 
are fufi propagators. Black dots are bare vertices and wfiite circles are full vertices, wfiicfi 
fiave to be reconstructed in tfie truncation discussed. 

wfiere T^^ are tfie tadpole contributions. Tfie first index gives tfie equation wfiere tfie 
tadpole contributes, G for gluon and H for Higgs, and tfie second index tfie type of 
tadpole appearing. Tree- level quantities are again denoted by a superscript 'tl', and can 
be found in equations ()B.14IIB.19|l in appendix |Bl In tfiese equations already tfie genuine 
two-loop contributions in tfie gluon and Higgs equations fiave been neglected. Tfie ansatze 
for tfie 3-gluon and Higgs-gluon vertices are motivated by tecfinical considerations and are 
tfius deferred to sections 14.41 and (4. 51 Tfie grapfiical representation of tfiis set of truncated 
equations is sfiown in figure 14.11 

Tfie Higgs propagator Dh is linked to tfie dressing function H by 

Dh{p') = (4.8) 

Replacing tfie propagators in ()4.5ll4.7j) by tfieir respective dressing functions, equations for 
tfie latter are obtained. Tfiese are divided by p"^ to make tfiem dimensionless. To obtain 
a scalar equation for tfie gluon dressing function, equation ()4.7|) is contracted witfi ()4.3|) 
and divided by 2. Tfiis results in 

W) ^ ^ + ^2j^^dqAT{p,q)G{q)Z{p-q) (4.9) 

^ = X _|_ _|_ J'^G _|_ rpHH _|_ 

H{p) 

+ j dedq[N,{p,q)H{q)Z{p + q) + N2{p,q)H{p + q)Z{q)) (4.10) 
l^T^H^r^GG^a^^j dedq(^R{p, q)G{q)G{p + q) 

+Ml{p, q)H{q)H{p + q) + Mrip, q)Z{q)Z{p + q)) . (4.11) 



Z{p) 
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Here 9 is the angle between p and q. At-, Ni, N2, R, Ml and Mt are the integral kernels 
for the employed truncation and are given in appendix 1(^.11 The tadpoles in the gluon 
equation are now also contracted. Although A^^i and N2 can be rearranged in one kernel 
by an integral transformation, they are left separated for a better comparison to the finite 
temperature case. Ca is the second Casimir of the adjoint representation of the gauge 
group. It only appears in the combination q^Ca, thus any change of the gauge group can 
be cast into a change of g^. Especially 't Hooft-scaling is therefore manifest. 

4.2 Asymptotic Analysis 

In this section, equations 1)4.9114.11]) will be solved analytically for asymptotically small 
and large momenta. In all these calculations, a massive Higgs will be assumed, as this is 
the case for Yang-Mills theories. However, a massless Higgs offers a rich infrared phase 
structure. Since this is somewhat out of the main line of interest it is deferred to appendix 



4.2.1 Ultraviolet Analysis 

Since asymptotic freedom is expected to hold also in the 3d case , the dressing functions 
and vertices should reduce to their one-loop counterpart and ultimatively to their tree-level 
value for sufficiently large momenta. 

As has dimension of mass and does not enter in the loop integrals in ()4. 9114. 11]) . 
by dimensional analysis the 1-loop contributions must be proportional to g'^/p- All loop- 
integrals are therefore sub-leading in the ultraviolet compared to the tree-level contribu- 
tion. In the case of the Higgs, this is only true if also p ^ rrih. This has been verified by 
explicit calculations, see appendix |El The 3d-theory gives therefore a very vivid example 
of asymptotic freedom. Thus all dressing functions D acquire their tree-level value in the 
ultraviolet. 



10:2! 



lim D{p) = 1. 



(4.12) 



Stated otherwise, to leading order for p ^ g^, rrih 



D{j>) ^ 1 + 



cglCA 



(4.13) 



P 



The constant c turns out to be positive for all dressing functions. Thus, all dressing 
functions approach the tree-level behavior from above. 
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4.2.2 Infrared Analysis 

To obtain analytical solutions in the infrared, the first step is to note that all integral 
kernels contain a term (p ± g)^ or (f in the denominator. Therefore the integrands are 
strongly peaked for p ^ q. To obtain analytic infrared solutions it is hence permissible to 
replace the dressing functions by their asymptotic infrared form. Since the infrared limit 
is a critical limit of the theory, the ansatz of power-laws is justified. The ansatze for the 
dressing functions are 

G{p) = Agp-^^ (4.14) 
Z{p)=A,p-^' (4.15) 
H{p) = Ahp-^K (4.16) 

In principle, three possibilities for the infrared behavior can be distinguished: Dominance 
of tree- level-, tadpole- and loop-terms. 

Tree-level dominance is motivated by the naive idea of a free gas of gluons in the 
asymptotic temperature limit. A first tempting ansatz is g — t — and I — —1 due 
to the explicit mass term for the Higgs. This corresponds to a purely perturbative or 
Coulomb phase. This is already impossible in the ghost equation, since in this case the 
ghost self-energy would diverge as 1/p in the infrared and hence superseding the tree-level 
term. The only possibility would be a ghost-gluon vertex suppressed at least like p in the 
infrared, which would be unexpected for tree-level propagators. Similar arguments apply 
to the gluon equation. Hence no Coulomb phase exists in 3d. 

If the gluon would behave as a massive particle in the infrared then the ghost equation 
would allow for a tree-level ghost. This is expected in a Higgs-phase, which would generate 
a magnetic screening mass. Such a term could only be generated by a tadpole diagram, 
since otherwise any of the possible vertices would have to diverge at least as strong as 1/p 
in the infrared (from the ghost-loop) or like 1/p^ (from the gluon- or Higgs- loop). This 
is not very probable with massive or tree-level particles. This should then also have to 
be true for arbitrary projections of the gluon equation, especially for ^ = 3. Since the 
tadpoles, even without truncation, drop out identically, a contradiction arises, and this 
solution is excluded. 

Attempting to save this option by using the ghost-loop to generate a mass also fails, 
since the required ghost exponent would lead to an non-renormalizable ultraviolet di- 
vergence of the integral. The only way to generate such a Higgs phase^ would be by 
appropriate fine-tuning of the ghost-gluon vertex and the ghost. This could of course 

^Besides spontaneous breaking of the gauge symmetry on the level of the Lagrangian by giving the 
Higgs a vacuum expectation value. 
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happen, but is unlikely. Besides, it is hard to see how the vacuum corresponding to such 
propagators avoids the violation of Elitzur's theorem [ZS]- Hence the tadpole must at least 
be compensated in the gluon equation and t < — 1, but integral convergence will require 
t < — 1, see appendix iDl 

Thus, the remaining option is loop dominance. Motivated by the reasoning of Zwanziger 
[in], ghost dominance is assumed, and the solutions indeed satisfy it. Hence neglecting 
all contributions without ghost lines, it is then only necessary to specify the ghost-gluon 
vertex. As discussed in section 13.31 a bare ghost-gluon vertex is chosen. As it is already 
required that 

t < -1, (4.17) 

confinement is present due to the Oehme-Zimmermann super-convergence relation ()2.25p . 

It is then possible to calculate the infrared limit of the DSEs (14.9114.11]) . This is done 
in appendix iDl and leads to 



A, 
A.. 



y 



-^'-'^'^-'-'iGTig.m.A, (4.18) 



l + y-i^-'^)/^-^HGG{9X)Al (4.19) 



4- = l + <±^, (4.20) 
Ah y ' 

where y = and d is the dimension. A subtraction in ()4.18p has been performed and only 
the finite part has been retained. The expressions for Igt and Igg stemming from the 
ghost-self energy and the ghost-loop, are calculated in appendix |Dl In the Higgs equation, 
a finite renormalization of the mass has been allowed for. The mass renormalization will 
be discussed in subsection 14.51 and is given in equation ()4.40|) . Equation ()4.2U|) is then 
solved immediately by setting Z = — 1 and 

m\ + 5m\ ' ^ ^ 

since 1 can be neglected for p <^ rrih- This already indicates that a qualitative change 
occurs in the high temperature limit, as in the vacuum t = I. This also immediately shows 
that the Higgs particle decouples, at least in the infrared, from the Yang- Mills sector. This 
agrees with corresponding findings on the lattice |H21- 

By dimensional consistency in the ghost equation ()4.18|1 . a relation between g and t 
follows directly jSj as 

ly ,4,2) 



2 V 2 / 2 V 2 
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Ghost exponent for different dimensions 



Ghost exponent as function of dimension 




Figure 4.2: The left panel shows the two solution branches of equation ()4.24j) as a function 
of ( within the region allowed by ()3.10p and integral convergence, see equation ()D.2|) and 
|13j. This excludes the second branch ind = 2. Dashed is 2d, solid is 3d and dashed-dotted 
is 4d. Note that at d = 4 there seems to be only one solution branch. It has been argued 
that the second branch is g = 1 [44 . The right panel shows the solution of ()4.24|) as a 
function of d for C = 1- 



The result at asymptotic temperature is hence different from the relation 

g = -2t 



(4.23) 



which is found at zero temperature. The additional power of p introduced by this change 
compensates the dimension of the effective coupling constant in 3d or the temperature in 
general finite-temperature 4d calculations. 

Due to ()4.17|) it then follows directly, that g is less than 0, and thus any solution 
will automatically satisfy ()2.2(i|l and therefore the Kugo-Ojima and the Zwanziger-Gribov 
conditions are both fulfilled. Also the tree-level contribution in the gluon equation in ()4.19|) 
can now be neglected, since the ghost-loop diverges in the infrared. Hence dominance of 
the gauge-fixing term as predicted by the Zwanziger-Gribov scenario is confirmed. 

Dividing equations ()4.19|1 and ()4.18|1 and eliminating t using ()4.22j) . a conditional equa- 
tion for g is obtained as 



22.-<i0F(| + ^)(2 + d{C - 2) - Ag{C - 1) - C) CSC (^^) sin(7r^^)r(f + g) 



((d-iygT{'-±^)Ti2g)) 



(4.24) 
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This equation has solutions in any dimension, see figure 14 .21 For = 3 it simphfies to 

32g{l-g){l-cot^{gn)) 



il + 2g)i3 + 2g)i2 + 2g{C-l)-Cy 
This equation has two solution branches. One is C-independent and yields 



(4.25) 



i9,t) = il,-l) (4.26) 

while the other one varies with (. In the special case of ( = 1, the other branch yields 

{g, t) ^ (0.3976, -1.2952) (4.27) 

thus reproducing the results of a previous analysis [44J which only regarded ( = 1. The 
second solution is a function varying significantly with (. It is therefore necessary to fix 
the allowed range of g. By virtue of the Gribov condition ()3.10|) and ()4.17|) as well as 
convergence of the integrals in the infrared, the allowed range for g in 3d is 

13 , , 

4 < ^ < 4- (4.28) 

It is further restricted by the requirement that a Fourier transform of the ghost propagator 
should exist, at least in the sense of a distribution, requiring g < 1/2. This restricts the 
range of allowed (^-values for the varying branch to 

1/4 < C < 3. (4.29) 

At the lower boundary both solutions merge into one. Note that the position space ghost 
propagator for (7 = 1/2 is a half-sided distribution, since its Fourier-transformed exists 
only in the sense of a limiting procedure. This is analogous to older expectations for the 
gluon propagator in 4d |H3] . 

Comparing the solution ()4.26|1 with the 4d-case 125] where 

l<9<h (4.30) 

the ghost exponent is only very weakly different from the 4d case, and the difference is 
even less pronounced for the gluon exponent. 

In equations ()4.18p and ()4.19|) the coefficients Ag and only appear in the product 
A'^Az- Therefore only this combination is determined by the infrared analysis and one of 
both coefficients has to be determined during the numerical solution procedure. This is a 
non-trivial problem, see appendix |Fl Expressing A^ by Ag yields 



1 AlCAgj 2<g"i)(2 + 2g{(: - 1) - C)r(2 - 2g) sm\ng) 
A. (47r)f cosi27Tg){g-l)g^T{l-2g) 



(4.31) 
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This result is used to check whether a correct solution is found during the numerical 
calculations. It also emphasizes that there is only one free parameter of the infrared 
solution. This one degree of freedom is necessary when using the full solution to perform 
the exact cancellation of the tree-level term in the ghost equation, which is the single 
crucial point in the derivation. 

With these results and appropriate regularization, it is possible to calculate the other 
contributions in all equations in the infrared limit as well. No contradiction is found. It 
would therefore be safe to proceed and solve the equations at all momenta numerically. 

However, both solutions can be realized at all momenta. Therefore the question arises 
which is the correct infrared solution or are both physical, characterizing different phases? 
In principle, the solution having the lower thermodynamic potential should be the correct 
one. This is investigated in chapterlHland prefers the varying branch aX ( = 1. Likewise the 
comparison with lattice calculations in subsection 14.61 prefer the varying branch at ( = 1. 
Furthermore, under the assumption of continuity in d, the solution in 3d and 4d should 
lie on the same solution branch. From the right panel of figure it can be inferred that 
again the varying solution lies on the same branch as the solution found in 4d. Therefore 
is is probable that the varying branch is the physical solution ai ( = 1. However, it is 
possible, that the presence of one of the solutions is only an artifact of the truncation. 
This cannot be ruled out in the current approximation scheme. 

4.3 Ghost-Loop-Only Truncation 

The ghost-loop-only truncation retains only tree-level expressions and loops with explicit 
ghost-lines. It thus implements the hypotheses of Zwanziger of dominance of the gauge- 
fixing contribution at all momenta. Calculations in 4d with this truncation [89] give 
qualitatively the same results as with the gluon-loop included f68j . 

The result of the calculations for both the ghost and the gluon at C = 3 are shown in 
figure l^ni In this case only the g = 1/2 infrared solution exists. The gluon propagator, 
for kinematical reasons, exhibits a significant maximum, with a center around p/g^ = 
0.25. The Higgs is in this case trivially tree-level, and therefore not shown. The infrared 
coefficients for all presented truncations are typically of order Agg^"^^ = O{10~^) and 
A^gf' = 0(10). 

Note that the 3d-theory is finite, and therefore no dimensional transmutation occurs. 
Hence the coupling constant is the only dimensionful quantity left. Measuring then all di- 
mensionful quantities in units of g^, the theory is without scale and should be independent 
of the value of the coupling constant. In this truncation scheme, this is the case, which 
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Ghost and gluon dressing function | 



Ghost and gluon propagator 



-g3G(p)/p' 
--g3Z(P)/p' 




Figure 4.3: The solutions of the DSEs in the ghost-loop-only truncation scheme at C = 3. 
The left panel shows the dressing functions and the right panel the propagators. The 
dashed line refers to Z and the solid line to G. All quantities have been made dimensionless 
by dividing out appropriate powers of g\. 



can been seen by rescaling the integration momenta in 1)4.9114.11]) appropriately. Thus, the 
plots shown in figure 14.31 are the results for any positive value of the coupling constant. 
Even for very small values strong non-perturbative effects for momenta smaller than g\ 
appear clearly. Hence this theory is never purely perturbative. 

For values of C different from 3, spurious divergences are present, and it is therefore 
necessary to include and adjust the tadpole term T'^^ such that these terms are canceled. 
To identify the correct contribution, the integration kernel R is split as 



^(C) = + (C - 3)i?3 + Rd{0. (4.32) 

where i?o and i?3 are independent of C and do not give rise to divergences, and Rd contains 
the divergent part for C 7^ 3. In perturbative calculations, the last term would be exactly 
canceled by the tadpole. However, it is not possible just to remove Rb for two reasons. 
First, Ry) contributes dominantly to the infrared, and the subtraction should not alter the 
infrared. Secondly, Rd is necessary to compensate an infrared divergence in the integration 
(not in the external momenta) of i?3 at C 7^ 3. Setting 

= 1 dqdQRn^V. q) (2G(|pl + |g1) - 1) (4.33) 



Chapter 4. Infinite- Temperature Limit 



41 



would be sufficient. This quantity contains no finite part for C = 3 and contributes only a 
small finite part compared to Rq and -R3 otherwise. 

However, even for a small deviation from = 3, no further solution can be found. This 
is not unexpected, since the ghost-loop should only dominate in the infrared and is only 
there restricted to be purely transverse. It is not transverse in the ultraviolet, as pertur- 
bative calculations show. However, as the corresponding compensating parts proportional 
to S^i, of the gluon-loop are neglected in this truncation, a failure was likely. Therefore 
using ()4.33p does not permit to go beyond C = 3 in the ghost-loop-only truncation. This 
behavior is not altered when using the second solution branch of the infrared. It is neces- 
sary to use a different subtraction scheme which removes these terms at finite momenta. 
This is possible by using 

R{p,q)G{q)G{p + q)^ 
Ro{p, q)G{q)G{p + q) + ((C - 3)Rs{p, q) + Rd{p, q))Alq-^3^p + q^^ 

+ (C-3)i?3(p,g) (4.34) 

or equivalently setting 

= j dqde{Rn{p,q){G{q)G{p + q)-Alq-'^{p + q)-'^) 

+ (C - 3)/?3(p, q) {G{q)G{p + g) - A^q-'^ip + q)-'^ - l)) . (4.35) 

The last term in ()4.34|) ensures the correct ultraviolet behavior to reproduce the 1-loop 
result. 

Note that by using ()4.34j) instead of ()4.33j) . also the solution at C = 3 is changed, 
since the second term in fl4.34|) misses the correct mid-momenta behavior. This change is 
however negligible, and the results would hardly be distinguishable from the plots shown in 
figure 14.31 This is the expected trade-off of the truncation, a deficiency at mid-momenta. 

Both solutions a.t ( = 1 are shown in figure 14.41 The main features visible are the 
following: The ghost dressing function is nearly featureless. It approaches its tree-level 
behavior in the ultraviolet from above. The gluon dressing function has a (non-visible) 
peak, and it approaches its tree-level behavior also from above. 

It is now possible to study different values of (. This gives a systematic error for the 
solution, since it corresponds to varying the amount of gauge symmetry violation. This 
is shown in figure 14.51 The variation for the g = 1/2 branch is obtained by varying ( 
between and 4. It is possible to obtain solutions for other (, but at some point the 
projection would be on the gauge violating longitudinal parts only and would hence give 
no information anymore on the transverse part. For the other branch ( is varied from 
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4.3. Ghost-Loop-Only Truncation 




Figure 4.4: The ghost and gluon dressing function at ^ = 1. The left panel shows the 
ghost dressing function, the right panel shows the gluon dressing function. The solid line 
denotes the solution for g = 1/2 and the dashed line gives the other solution branch with 
g PS 0.4. 



1/4 up to 2.55, the latter corresponding to g = 0.2708 and t = —1.042. In principle, it 
should be possible to extend the variation over the complete range ()4.29|1 . However, it 
was not possible numerically to go to larger (, since the available numerical precision of 
16 significant digits was not sufficient. Once higher precision is available the remaining 
solutions should be obtainable. The reason is that as t comes close to —1, the dressing 
function more and more resembles a massive solution. This corresponds to a jump from 
to a finite quantity in the corresponding propagator. Hence, the zero is only reached far in 
the infrared, and it is numerically necessary to have at the same time very large momenta 
from the integration and very small ones from the infrared, a combination limited by 
precision. This is also the case when including the gluon- or Higgs-loop. 

Looking now in detail at the ^-dependence, it can be observed that there are only 
quantitative, but no qualitative changes. The vanishing of the peak in the gluon dressing 
function for ( > 2 is a purely perturbative phenomenon, since the contribution of the 
ghost-loop changes sign at C = 2 in the ultraviolet. The changes in the infrared are partly 
due to changes of the coefficients Ag and Az- The variable branch shows a much larger 
variation as expected. The deviation is comparatively small when going to smaller (. It 
increases significantly when increasing (, since the gluon propagator becomes more and 
more like the one of a massive particle. The propagator always vanishes at zero, though. 
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Figure 4.5: The ghost and gluon dressing function in ghost-loop-only approximation for 
different values of (. The left panel shows the ghost dressing function, the right panel 
the gluon dressing function. The solid lines represent the solutions with g = 1/2, and the 
dashed lines for the varying branch. At the peak in the gluon, the middle lines correspond 
to solutions at C = 1- The upper and lower line at mid- momenta give the solutions at 
C = and 4 for the g = 1/2-branch and at C = 1/4 and 2.55 for the other branch. Note 
the linear scale for Z in the right panel. 



4.4 Yang- Mills Theory 

The next step is to restrict the calculation to the Yang-Mills subsector without Higgs. 
This is then a theory defined on its own. Starting again at C = 3 no solution can be found. 
Even trying a significant diversity of subtraction schemes to do calculations at other values 
of ( does not help, nor does employing the other solution branch. 

Closer inspection reveals the source for these problems. The ghost- and gluon-loop 
contribute with opposite sign. This is correct in the ultraviolet and irrelevant in the 
infrared, where the gluon-loop is subleading. However, at momenta of the order of g^, the 
gluon loop is able to dominate the ghost-loop and the tree-level term, and thus leads to 
a violation of the Gribov condition (jH.lOj) . Therefore, within this truncation scheme, the 
tree-level 3-gluon vertex is not an acceptable truncation. The reason is that either the bare 
three-gluon vertex overestimates the true vertex at mid-momenta, the bare ghost-gluon 
vertex underestimates the true one at mid-momenta, or the contributions of the two-loop 
graphs are important. A finite contribution from the tadpole would not help, since the 
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problem persists even at ^ = 3 where its contribution vanishes identically^. This is clearly 
an artifact of the truncation. Any of these options deserves further attention, and as a 
first step, currently the ghost-gluon vertex is under closer investigation [TFj . 

For the purpose of this work, the specific reason for this mid-momenta deficiency is 
of minor importance. For further calculations it is only necessary to improve the vertex 
such that it provides a sufficient suppression at mid-momenta. This situation is similar to 
the case of 4d, where without a change of either vertex, the correct ultraviolet behavior 
is not reproduced j3T]. Since the only relevant effect is at mid-momenta, such a vertex 
construction has to fulfill two criteria: It has to be tree-level in the ultraviolet and not 
more strongly divergent than on either leg in the infrared. 

The simplest solution is to multiply the bare 3-gluon vertex with an additional function 
delivering this suppression. This leads to the ansatz 



This ansatz is motivated by the vertex used in 4d jHHj, with the anomalous dimensions 
set to 0, as appropriate for 3d. Then an appropriate structure has been selected so that 
A is constant in the infrared and tree- level in the ultraviolet. The additional parameter S 
smoothly interpolates between a tree-level vertex and a large suppression by selecting it 
from [0,oo). This ansatz conserves Bose symmetry, as opposed to the case of 4d, where 
this was not possible without spoiling the ultraviolet properties. 

Alternative concepts have been studied. The simplest possibility to remedy the situa- 
tion, albeit at the expense of an incorrect ultraviolet behavior, is to multiply the vertex just 
with a constant, thus introducing a vertex renormalization Zi. Qualitatively the same and 
quantitatively similar solutions have been found for Zi ^ 0.6. The Bose non-symmetric 
ansatz A{q, G, Z)A{q + p, G, Z), similar to the 4d-case, has also been made. Although it 
delivers a smaller suppression at the same value of 6 as ()4.3(j|l . the solution is possible 
for sufficiently large values of 6. It does not differ qualitatively from the results obtained 
using ()4.36|) . Even the quantitative changes are small. The last option investigated was 
modeling the two- loop terms by adding a term which is Gaussian in momentum to 1)4.111) . 
As long as its strength around was sufficient, the system was solvable, again with only 
minor modifications. Thus it is likely that the exact solution does not change the results 

■^Albeit a scaleless integral also in 3d, it is not necessarily a purely divergent quantity. In 4d, it is 
possible to absorb any finite contributions in the renormalization constant. Lacking a renormalization 
process, this is not possible in 3d. 



itl, A^,abc 

l3(Tfl 



i-q, q + p, -p) iA{q, G, Z)A{q + p, G, Z)A{p, G, Z)Y 



A{q,G,Z) = 



(4.36) 
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Figure 4.6: The vertex and the gluon loop with and without modifications. The left panel 
shows the function A from ()4.36p . The solid line is 5 = 0, the tree-level vertex. The dashed 
line is the g = 1/2 solution and the dotted line the g ~ 0.4 solution. From bottom to top 
these are S = 1, 6 = 1/4 and in the last case S = 0.14 and 6 = 0.12, respectively. The 
right panel shows at C = 3 the full gluon loop compared to the gluon loop suppressed by 
the vertex ()4.36|) and subtracted by the gluon part of the tadpole ()4.38|) . again for both 
solutions with the same type of line. The larger values are always without suppression. 



qualitatively. The approach taken with ()4.3fij) is supported by lattice calculations of the 
dressing functions, see section I4.fjl 

An interesting possibility to assess the effect of different vertices is to change 6. While 
6=1 corresponds to a significant suppression, the smallest values with still a stable solu- 
tion are 0.131 for the g = 1/2 solution and 0.114 for the other one. As a convenient value 
which still provides very stable numerics and at the same time not too large suppression, 
5 = 1/4 is arbitrarily chosen in the following. The structure of the vertex and its impact 
on the gluon loop are studied for different values of 6 in figure I4.fj[ 

It is now possible to solve the equations at ( = 3. The results are shown in comparison 
to the ghost-loop-only results in figure 14.71 The only qualitative difference compared to 
the ghost-loop-only truncation is that a peak in the gluon dressing function is now present 
for all ( due to the now correctly reproduced LO perturbation theory for p ^ g^- 

After this first view at the full Yang-Mills solution, it is again interesting to study its 
dependence on (. In principle, it would be sufficient again to split off the divergent part 
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Figure 4.7: Solution of the Yang-Mills sector at ( = 3. The left and right panel shows 
the ghost and gluon dressing function, respectively. The dashed curve gives the ghost- 
loop-only solution for comparison. The solid curve gives the full Yang-Mills solution at 
6 = 1/4. Note the linear scale for Z in the right panel. 

Mtd of Mt and subtract it in the same way as in ()4.33|) . by setting 

= -^Jdqd9[Rnip,q){2Gm + m-^) 

+ {A{q, G, Z)A{q + p, G, Z)A{p, G, Z)YMtd{p, q) {2Z{\p\ + - 1) ) .(4.37) 

The naive expectation is that after the inclusion of the gluon loop the spurious terms 
would be ehminated. This does not occur, since by introducing an ad-hoc suppression, 
these cancellations are affected. Hence the same problem arises again, and it is necessary 
to correct both loops individually. This can be done by 

= ~^ldqde[Rnip,q)iGiq)G{p + q)-Alq''^ip + q)-'<^) 

+ (C - 3)Rs{k, q) {G{q)G{p + q) - A^q^'^ip + q)"'^ - l) 

+ {A{q, G, Z)A{q + p, G, Z)A{p, G, Z)YMtd{p, q)Z{q)Z{p + g)) , (4.38) 

i.e. by completely removing the term Mtd and keeping the same subtraction as before 
for the ghost-loop. Removing is not harmful to the infrared, since the gluon-loop 

is subdominant there. It is also irrelevant in the ultraviolet, where it contributes only a 
pure divergence and would therefore anyhow be canceled by the tadpole. Thus, only the 
mid-momenta behavior is altered. 
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Figure 4.8: Solution of the Yang-Mills sector at C = 1. The left panel shows the ghost 
dressing function and the right panel the gluon dressing function. The solid curve denotes 
the (7 = 1/2 solution, while the dashed curve displays the g ~ 0.4 solution. Both are at 
5 = 1/4. 



It is now possible to solve for any (. Both infrared solutions are found to connect to 
full solutions and a comparison at ( = 1 and 5 = 1 /4 is shown in figure 14.81 The difference 
is again only quantitative in nature, and even this difference is quite small. 

It remains to assess the dependence on the parameters ( and 5 as a systematic error 
estimate. While the dependence on ( summarizes the impact of violation of gauge sym- 
metry, the dependence on 6 exhibits the uncertainty in the correct mid-momenta behavior 
of the vertices. Both problems are linked, and do not cause independent errors. However, 
both generate different changes when varied. Figure 14.91 and 14.101 display the dependence 
on ( and 6, respectively. 

The uncertainty in the structure of the vertices is much more significant at mid- 
momenta than the violation of the STI ()2.18|1 . On the other hand, the uncertainty in 
the vertex only amounts to a shift in the infrared but not to a change of the infrared 
exponents as in the case of the violation of the STI. 

The height and the position of the peak depends weakly on the vertex construction. 
The height of the peak can be traced to the amount of suppression the vertex gives. When 
the suppression is weak enough, the peak tends to become a divergence with sign change. 
Therefore, this truncation scheme can give no definite answer to the size of the peak. It 
turns out that the peak only starts to get very large when 6 gets close to the critical value. 
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Figure 4.9: The ghost and gluon dressing functions of Yang-Mills theory for different 
values of (. The left panel shows the ghost dressing function, and the right panel the 
gluon dressing function. The solid line gives the solution for g = 1/2 and the dashed line 
is for the other solution branch. At the peak in the gluon dressing function, the middle 
lines represent the solution at ( = 1. The upper and lower lines at mid- momenta give the 
solutions at C = and = 4 for the g = 1/2-branch and at ( = 1/4 and ( = 2.45 for the 
other solution branch. 



Before that, it is stable at height around 1.5 and below. Hence, it is to be expected that 
the correct suppression would also not deliver a significantly larger peak. This is confirmed 
by lattice results, as will be shown in subsection 14.61 

Comparing these results with those of the 4d-theory [22] , little difference is found. On 
qualitative grounds, both are the same, besides renormalization. They only differ in their 
quantitative nature. Thus the 3d Yang-Mills theory is a strongly interacting, confining 
theory, confirming the expectation |88j . 



4.5 Full 3d-Limit 

Finally, by adding the Higgs, the full system of equations ()4.9M.ll]l is implemented. Al- 
though the contribution of the Higgs-loop in ()4.9j) is positive, it is not sufficient to allow 
for a solution with a bare 3-gluon vertex. Hence the vertex ansatz and also the tadpole 
construction will be kept the same as in the Yang-Mills case. 

In addition, as no further input is available, the tree-level gluon-Higgs vertex ()B.17jl 
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Figure 4.10: The ghost and gluon dressing functions of Yang-Mills theory for different 
values of 6 at ( = 1. The left panel shows the ghost dressing function, and the right panel 
the gluon dressing function. The solid line gives the solution for g = 1/2 and the dashed 
line for g ^ 0.4. At the peak in the gluon dressing function, the middle lines represent 
the solution at 6 = 1/4. The upper and lower lines at mid- momenta give the solutions at 
S = 0.131 and 6 = 1 for the g = 1/2-branch and at 5 = 0.114 and 6 = 1 for the other 
solution branch. 



is employed. A closer inspection of the consequences of this choice is made in subsection 

o 

The ratio ruh/gl remains to be fixed. In the pure 3d-theory, nih is a parameter and 
can therefore not be determined a priori. It stems from integrating out the hard modes 
when making the 3d-approximation |SS1- By taking these modes into account in chapter 
it will be possible to investigate it more closely rather than making an ansatz. For the 
investigation in this chapter, the origin (and also the exact value) of the Higgs mass is 
not of direct importance, and hence will be fixed to mh/ g1 = 0.8808, a value extracted by 
lattice calculations and perturbative fitting |H2!- 

By adding the Higgs, three more tadpole contributions arise, as can be seen from figure 
13. 31 and in appendix |B| Since the Higgs has a tree- level mass, two of the tadpoles can have 
a non- vanishing finite part already in leading-order perturbation theory, see appendix lEl 
Hence, a more thorough discussion is needed here. 

The first additional tadpole appears in the gluon equation. On the one hand, all 
diverging contributions must cancel again due to the STI, hence the structure of diver- 
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Figure 4.11: The solution of the Yang-Mills system compared to the full system at C = 1- 
The solid line denotes the full solution at = 1/2 while the dotted line is the corresponding 
Yang-Mills solution. The dashed line gives the full solution for g ^ OA while the dashed- 
dotted line is the corresponding Yang-Mills solution. 

gencies must be the same for all graphs. On the other hand, when projecting with the 
Brown- Pennington projector ()4.3|1 . the tadpole identically drops out. Therefore all finite 
contributions of the tadpole must also be contained in the other loops as well. For a 
vanishing tree-level mass, the tadpole is canceled by divergent loop contributions. It is 
hence a good assumption that any finite part of the tadpole will also be canceled by the 
loop contributions. Therefore the approximation concerning the Higgs-tadpole will be 

T^"" = J dqde{MLD{Kq)H{q)H{k + q)) (4.39) 

where Mld is the divergent part of M^. If this were not justified, at ^ = 3 a trace of 
the finite contribution should be visible. Any such term would behave like 1/p^, a clear 
signature. The gluon- and Higgs-loop both are at most constant in the infrared and can 
thus not induce such a contribution. Only the ghost loop with its strongly divergent 
behavior could contain it. However, it seems unlikely that the finite part of the tadpole 
should be completely moved to the ghost loop, hence the assumption described above will 
be made. 

The Higgs equation with two tadpoles is more problematic, as the argument using the 
projection is not available. As the Higgs is a component of a 4d gluon, its self-energy 
has to be finite. Thus the divergent parts of the tadpoles have to cancel. This fixes the 
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Figure 4.12: The left panel shows the Higgs propagator and the right panel its dressing 
function, both at ^ = 1. The solid line gives the solution for g = 1/2, the dashed line 
for g ^ 0.4, the dashed-dotted line denotes the leading-order perturbative result and the 
dotted line the tree- level behavior. 

tree-level coupling h at leading-order perturbation theory, as given by ()4.4|1 . The only 
remaining question is that of a possible finite contribution. Such a contribution would 
only change the mass of the Higgs, and this occurs already at leading-order perturbation 
theory. However, the subtraction of two divergent quantities is not simple, and due to the 
error induced by the truncation somewhat arbitrary. Hence the mass-shift is essentially 
unknown. Nonetheless, the Higgs dressing function behaves similar to perturbation theory. 
Therefore assuming the finite part of the tadpole to be the same as in leading-order 
perturbation theory seems to be justified. Thus, the tadpoles are fixed to be 



p2 p2 



defining the finite mass renormalization 5m^ of the bare Higgs mass. This amounts to a 
change of the mass of about 10%. 

The impact of the Higgs on the Yang-Mills sector at C = 1 is shown in figure H^TTl The 
infrared behavior is only weakly affected and the situation is similar in the ultraviolet. 
The largest impact on the gluon dressing function is on the peak height and the impact 
on the peak position is significantly weaker. Altogether, the effect of the additional Higgs 
field is only quantitative, as expected from lattice calculations [H2]- 

The Higgs propagator and dressing function are shown in figure 14.121 Both show a 
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Figure 4.13: The left panel shows the gluon dressing function and the right panel the Higgs 
propagator as a function of rrih. The gluon peak increases with decreasing mass, while 
the Higgs propagator decreases with increasing mass. Solid is the g = 1/2 and dashed the 
g ^ 0.4 solution. The masses are 2(yf|, O.SSg"^ for both solutions, and 0.75g'^ and 0.68gl, 
respectively. 



behavior very similar to the leading-order perturbative result, and the difference to a 
tree-level propagator is mainly due to the mass-shift by the tadpole. 

Although the mass is fixed in the current setting, the question of the dependence on 
the mass is still interesting. Making the tree-level mass larger, the solution approaches 
smoothly the Yang-Mills solution, with which it merges as — oo. On the other hand, 
the mass cannot be made arbitrarily small. When the mass is reduced below mh/g"^ ^0.7, 
no solution can be found without changing 6. The reason is that the Higgs self-energy 
is negative, and only while it is sufficiently small compared to the renormalized tree-level 
contribution, it is possible to maintain the Gribov condition ()3.10|) at mid-momenta. This 
additional constraint only stems from the 4d origin of the theory and the fact that the Higgs 
is part of a 4d gluon, and is not inherent to the 3d-theory. In the 3d-theory, a negative 
dressing function of the Higgs would be allowed, and smaller masses could be reached. The 
height of the gluon peak directly corresponds to the size of the Higgs-self-energy and hence 
to the lowest masses achievable without violating ()3.1U|) . Thus, increasing 5 compensates 
to some extent a lowering of the mass, but this has only been pursued up to 5 = 1 or 
a lowest mass of rrih/g^ ~ 0.6. The dependence of the Higgs propagator on the mass is 
depicted in figure 14.131 An effect at small Higgs masses had been anticipated from the 
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Figure 4.14: The Higgs dressing function and propagator as a function of 6 at ( = 1. The 
sohd hne is the solution for g = 1/2, the dashed hne for g ^ 0.4. The higher the peak, the 
larger 6. The lowest peak corresponds to 5 = 1, the middle one to 6 = 1/4 and the largest 
one to 6 = 0.0862 and 6 = 0.061, respectively. 



effect of a massless Higgs on the infrared solutions, see appendix ID. 21 

The dependence of the Higgs on ( is very small, on the order of a few percent at 
mid-momenta. The (^-dependence of the Yang-Mills sector is also only affected to this 
extent by the presence of the Higgs. The dependence on 6, which influences the size of 
the self-energy, is more pronounced. The latter is shown in figure 14.141 Similar to the 
gluon dressing function, in the Higgs dressing function the peak increases with decreasing 
S, at some point also generating a peak in the propagator. This is due to the direct 
link of the Higgs self-energy to the peak volume of the gluon propagator. As the latter 
increases, the Higgs-self-energy increases. Since it is negative, it starts to compensate 
the tree-level term at mid-momenta, where the effect of the peak in the gluon dressing 
function is most pronounced and thus increases the peak in the Higgs dressing function. 
At some point a further effect sets in, when the system tries to compensate the increasing 
negative contribution from the gluon loop by an increasing positive contribution in the 
Higgs-loop. This is most easily achieved by an enhancement of the Higgs dressing function 
at mid-momenta, thus also increasing the peak in the Higgs dressing function. This process 
slows the growth and for sufficiently small values of 6 even decreases the peak in the gluon 
dressing function. Therefore it becomes likely that a violation of the Gribov condition 
()H.10|1 will occur due to the Higgs rather than due to the gluon. This behavior is much 
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Figure 4.15: The dependence of the infrared coefficient Ag on ( in the left panel and on 6 
in the right panel. Solid lines are the g = 1/2 solution and dashed lines the other solution 
branch. 

more pronounced in the g = 1/2 solution. 

Apart from the dependence on ( of the exponents in the case of the second branch, 
also the dependence of the infrared coefficients on ( and 6 is of interest. Ah is fixed by 
the renormalized mass and A^ depends uniquely on Ag, hence only the dependence of Ag 
is relevant. Its dependence is shown in figure I^?TH1 The dependence on ( is rather weak. 
Both solution branches evolve smoothly out of one solution at ( = 1/4. The dependence 
on 6 is much more pronounced, where the trend to become singular is visible in the strong 
increase at low values of 6. The bending-over is the aforementioned effect of the peak 
growth switching from the gluon to the Higgs. 



4.6 Comparison to Lattice Results 

Currently, only lattice calculations for the gluon propagator [771 IH2] and the Higgs propa- 
gator jH2l are available, but not for the ghost propagator. Furthermore, lattice calculations 
do not yet see a significant difference in the gluon propagator with or without the Higgs 
[H21- Hence the lattice calculations for the gluon propagator are done in pure 3d Yang- 
Mills theory. Therefore, these results have to be compared to the results from section 14.41 
rather than to the full calculations in section 14.51 The comparison is made in figure 14.161 
Comparing the gluon propagator to the lattice results, the assumed suppression using 
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Gluon dressing functiori~t 



Higgs dressing function~~| 




Figure 4.16: The left panels show the gluon propagator and dressing function in Yang- 
Mills theory from lattice calculations and from section 14.41 The continuum-extrapolated 
values are from jHZI, the others from [ZZj. The right panel shows the Higgs propagator 
and dressing function in the full theory from lattice calculations jH21 and from section 1^31 
The errors indicated are statistical. 
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5 = 1/4 seems to fit the data quite well. The dressing function is better fitted by the 
(7 = 1/2 solution while the propagator and especially the infrared favors g ~ 0.4. At 
momenta p ^ g^, the truncation scheme is not trustworthy and the better agreement of 
the g = 1/2 solution with the lattice gluon dressing function is not conclusive. Although 
the lattice data for the propagator are quite scattered, they lie pretty close in case of the 
dressing function. The propagator is hence more sensitive than the dressing function. A 
better fit could be achieved by varying 6. This is not the aim of this work, but may be of 
interest for phenomenological studies. 

The lattice result for the Higgs propagator differ more from the DSE result than for the 
gluon propagator. Especially, the lattice Higgs propagator is much more similar to the tree- 
level result than to the leading-order perturbative result. The reason for the discrepancy 
between the current calculation and the lattice data may be twofold. The value at zero 
momentum indicates that ()4.4()j) overestimates the mass renormalization and the mass 
found is too large by roughly 5%. More severe is that the bare vertex used here seems to 
significantly overestimate the gluon-Higgs interaction compared to the fully dressed vertex 
of the lattice result. The self-energy must be significantly suppressed even for large values 
of p to allow for such a deviation from the leading-order perturbative result, and must be 
nearly negligible. In principle, such an effect could be modeled by a construction similar 
to ()4.36|) for the Higgs-gluon vertex. Again, this is not the aim of this work. However, 
this finding entails that strong subleading or even genuinely non-perturbative effects are 
also present in the Higgs sector. Even for quite large momenta of the order of lOg^, where 
the gluon propagator has already become perturbative, this is not the case for the Higgs 
propagator, which will only become perturbative at much larger momenta of the order 
of 100(73 . Indeed, the overestimation of the mass already shifts the result towards the 
tree-level propagator, and a smaller mass would enhance the discrepancy. This is indeed 
quite an interesting effect, and it will be investigated in some more detail in section IfTTl 



4.7 Coulomb Gauge Instantaneous Potential 

Recently, a connection between Coulomb-gauge in 4d and Landau-gauge in 3d has been 
established [HHl 1^ . It was shown jHS] that under certain approximations, the static 
quark- ant i- quark potential, which is linked to the 00-component of the Coulomb gauge 
gluon propagator, is given in four dimensions by 

V,oM = ^^ + VUp)- (4.41) 
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Vcon designates contributions which stem from connected parts of expectation values of 
the Faddeev-Popov operator. Neglecting this part, the potential is solely given by the 
function G. It can be shown to be identical to the ghost dressing function of the 3d 
Landau gauge Yang-Mills theory jH3J. While the solution g ^ 0.4 generates a potential 
which behaves as ~ and thus a little less than linear, the second branch generates 

a solution proportional to for small momenta, thus generating the behavior expected 
for a linear confining potential. 

In principle it is expected that after Fourier-transformation ()4.4H) would give the lead- 
ing contribution of the static potential. However, it is not possible in the given approx- 
imation scheme to obtain the Fourier transform. The Fourier integral is quadratically 
divergent in the ultraviolet. The reason is that G goes to 1 in the ultraviolet, thus to- 
gether with the measure it generates a contribution proportional to p. Using a prescription 
as will be used in section it may still be possible to extract some information, but since 
not a constant but a function needs to be fitted, this would turn out to be significantly 
more complicated. It is improbable that the numerical accuracy would be sufficient to ex- 
tract relevant information, as it is already hard to do so in the case of a constant quantity 
like the thermodynamic potential treated in section This is due to the approximations 
involved when obtaining the infrared limit ()4.4H1 . 



Chapter 5 

Finite- Temperature Effects 



Lowering the temperature down from infinity also the hard modes must be included. The 
corresponding DSEs are derived in section IKTTl The precise implementation of the limiting 
process needs some deeper analysis of the running coupling, which will be conducted in 
section The infrared properties of the dressing functions will be investigated in section 
15.31 The truncation will be discussed in section In general, an infinite number of Mat- 
subara frequencies contributes at finite temperatures. At sufficiently high temperatures, 
most of them can be neglected, however. This corresponds to a small-momentum approx- 
imation, to be defined in section 1^31 Also the renormalization process will be discussed 
there. Results on the ghost-loop-only truncation and the full theory will finally be shown 
in sections l5 . 61 and l5T| respectively. The chapter will be closed in section ESI by presenting 
results for low temperatures, which have been obtained in the context of this work. There, 
also a comparison to the high-temperature case will be drawn. 

For convenience, the 3d-longitudinal and 3d-transverse gluons will be referred to as 
longitudinal and transverse gluons. In this chapter is again = p^+fP' = vl+pi- External 
momenta are referred to as p while internal loop momenta and Matsubara summation 
momenta are referred to by g. 



5.1 Finite- Temperature Dyson-Schwinger Equations 

As described in section 13. 2[ the equations are obtained from the vacuum equations ()3.2|1 
and ()3.3|) by application of the Matsubara formalism J^. To obtain scalar equations 
for the (infinite) set of dressing functions Z and H of the gluon propagator ()3.5p . the 
gluon equation is contracted with the generalized projectors ()3.8|1 and ()3.9j) . Due to the 
kinematical singularity present, the choices 1)4.11] and ()4.2|1 are extremely cumbersome. 
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Therefore 

ATfMu = S^u-5^o5ou (5.1) 

. r POPU . . P^PO ,^ ^^ 

AL^,v = V — — + ^ou — ^ (5.2) 
are used instead. This is a speciahzation of the most general projector 

PtiPu 



p2 



(5.3) 



where rj can be either ^ or ^. The choices of the constants t^, f, ft, Ct, and in ()5.1|) 
and ()5.2|1 were made to obtain a well-defined 3d-limit. They are not unique, indicating 
a larger range of possibilities than just one parameter rj in general. Nevertheless, as the 
choices fulfill all requirements they are sufficient. 
Employing the projectors ()3.8p and ()3.9|) yields 



7 ^ glTCA 

^3 + 



2 



oo „ 

J dedq(^AT{p,q)G{q)Z{p-q) 



G{p) " i2n) 

+ALip,q)Giq)Hip-q)) (5.4) 

= e^3^ + T^<? + T^^ + |^ j d9dq{p{p,q)G{q)G{p+q) 

+Nl{p, q)Z{q)Z{p + q) + N,{p, q)H{q)Z{p + q) + N^ip, q)H{p + q)Z{q) 
+NT{p,q)H{q)H{p + q)^ (5.5) 

-1- = Z3T + T«^ + T«^ + |% f; [dedq(R{p,q)G{q)G{p + q) 

+Mi(p, q)H{q)H{p + q) + M^{p, q)H{q)Z{p + q) + M^^p, q)H{p + q)Z{q) 
+Mrip, q)Z{q)Z{p + q)) + ^Z,, - ^) . (5.6) 

The C and ^ dependence is acquired by using ()5.H) and ()5.2|) . At C = ^ = 1 the original 
form of the equations is recovered. For the correct solutions, the appearance of these 
terms would be compensated by corresponding structures in the loop contributions, and 
the result would be independent of ( and ^. For po = 0, equation ()5.5|1 is only superficially 
dependent on ^. As all integral kernels are proportional to ^, the dependence can be 
divided out for ^ 7^ 0. Then, only the implicit dependence due to the interactions with 
the hard modes remains. The latter effect vanishes in the 3d-limit and the Higgs-equation 
()4.6|1 is again independent of the projection and thus unique. The finite-temperature 4d- 
theory is no longer finite, but renormalizable. Thus explicit wave-function renormalization 
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constants Z3, Z^l, and Z^t have been introduced and will be discussed in section 15.5.21 
and l5.5.:-il 

The summation is over all Matsubara frequencies go = ^irTn. Z\ = \ has been em- 
ployed Z\ cannot be calculated in the present approach, but it is finite and set to 
1 IHHllZni- The kernels At, Al, R, Mt, Mi, M2, Ml, P, Nt, Ni, N2, and Nl are listed 
in appendix IC.2I The tadpoles T*-' are again used to cancel spurious divergences in much 
the same way as in the previous section. However, also in the soft equations, they can 
possess finite parts at finite temperature. In case of the longitudinal equation, these can be 
absorbed in the mass renormalization discussed below. In the transverse equation, these 
are at ( = 3 completely contained in the loop terms, and their continuation away from 
C = 3 is thus arbitrary. Also these contributions scale at best only as and are thus 
irrelevant to the infrared. Therefore, they are dropped, especially as no simple prescription 
as ()5.1()|1 below can remove the related spurious divergences 

Note the considerably higher symmetry between the equations for H and Z than in 
equations ()4.9M.ll]l . Close inspection of the equations reveals further that the dressing 
functions can only depend on |po| and The corresponding symmetry, under po —po, 
is used to reduce the number of equations significantly. 

5.2 Temperature Dependence of the Coupling 

The equations ()5.4II5.6|) depend on the coupling constant as the only parameter. In turn 
g4 depends on the renormalization scale /x, which can be chosen arbitrarily at any fixed 
temperature T. It directly enters into the definition of the infinite-temperature limit, as 

depends on g^. In the simplest case, g^ ~ gl{fi)T. Comparing equations ()4.9ll4.1ljl and 
fj5.4H5.6p . the constant of proportionality has to be 1 in this truncation scheme to obtain 
a smooth infinite-temperature limit. Hence it only remains to choose the temperature 
dependence of /x. There are two options to consider. 

If /i is fixed, the 3d-coupling grows without limit as T grows. This does not necessarily 
pose a problem, as the infinite-temperature propagators are independent of 5^3 as long 
as expressed as a function of p/g^- Under such circumstances, all momenta below T are 
effectively infrared, and the non-perturbative regime would extend to all momenta. 

Alternatively, it is possible to use the limit prescription gl{fj,(T))T = Coo, where Coo is 
an arbitrary constant, effectively performing a renormalization group transformation when 
changing the temperature. This fixes g^ to be proportional to Aqcd, the dynamical scale 
of QCD jTj, while the 4d-coupling vanishes like 1/T. This defines a t'Hooft-like scaling 
in T. Note that —> for T — 00 corresponds to the conventional arguments of a 
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vanishing coupling in the high-temperature phase [HI]. However here, as in the large- 
limit, t'Hooft-scaling is performed, generating a well-defined theory. Hence this possibility 
defines a smooth 3d-limit with a finite 3d-coupling constant. Here Cqo = 1 is chosen. 
Nonetheless, the case of fixed /i is also of interest, as it permits to compare the results for 
different temperatures directly. 

These possibilities all correspond to the 4d-situation, where it is possible to shift the 
onset of non-perturbative physics by renormalization group transformations. Under all 
circumstances, for any finite temperature, there is always a non-perturbative regime at 
sufficiently small momenta. 

Furthermore, it would be useful to give explicit units for the temperature scale. In 
the vacuum the scale is fixed via a comparison of the running coupling to perturbation 
theory [HHl- This is not possible here because the Matsubara sum is truncated and the 
coupling cannot be calculated for momenta of the order of 2nNT, where is the number 
of Matsubara frequencies included^ 

The simplest procedure is to compare to lattice calculations. There the effective 
3d-coupling is found to be g1{2Tc) /2Tc = 2.83 82 . The phase transition temperature 
is Tc = 269 ± 1 MeV fP. Albeit this will not be exactly gfT^ due to the truncation, 
a first approximation is to require the same temperature scale at 2Tc. Using the fixed 
gfT = Coo prescription, then 1 = Coo/(2Tc2.83). Since here Coo = 1 in internal units, 
internal units have to be multiplied by 1.5 GeV to yield physical units. This temperature 
scale will be used in sections 15.61 and 15.71 The ratios of temperatures are independent of 
this prescription. 

5.3 Infrared Properties 

Regarding the infrared properties, asymptotic freedom turns out to be advantageous also 
for the infrared. Just at the phase transition, the n = 1 Matsubara frequency has already 
an effective 'mass' Pq = 2ttTc ~ 1.7 GeV. By the Appelquist-Carazzone theorem the 
hard modes are suppressed by powers of |p|/po in the infrared. Thus hard modes do not 
contribute significantly. Studying each equation in detail, it turns out that this quite 
general statement is implemented very differently in each equation. 

In case of the hard modes, there is no pure soft contribution due to momentum conser- 
vation at the vertices. Thus they decouple and a self-consistent solution is that all hard 
mode dressing functions are constant in the infrared. 

In the case of the soft ghost equation ()5.4j) . all contributions become also constant in 

^Always, only the frequencies with po > are counted. 
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the infrared. Thus they can be canceled by the (renormahzed) tree-level value. Only the 
subleading behavior remains, as in the case of the 3d- and 4d-theory. As the subleading 
contribution of the purely soft term dominates the hard terms, the same behavior as in 
the infinite-temperature limit emerges. 

In case of the transverse equation ()5.6|) . the hard mode contributions give rise at best 
to mass-like terms. These terms are subleading compared to the soft ghost-loop. 
Thus the infrared is dominated by the latter and the same infrared solution as in the 
infinite-temperature limit is found. 

The longitudinal equation ()5.5p is finally quite different from its equivalent in the 3d- 
or 4d-case. In the prior case, it was dominated by its tree-level mass, while in the latter it 
is identical to the transverse one. In the present case, the absence of pure soft interactions 
with ghosts, which could generate divergent contributions as in the transverse equation, 
leads to dominance of the hard modes^. These provide mass-like contributions in the 
infrared. Inspecting e.g. P at ,^ = 1 



and using the fact that the hard dressing functions are constant in the infrared, directly 
leads to a mass-like behavior due to the explicit factor and the finiteness of the 

remaining expression. Thus the 3d-mass of the Higgs is generated spontaneously by the 
interaction with the hard modes. 

Hence, the infrared behavior of all soft mode dressing functions does not depend on 
the temperature, expect for changes in the infrared coefficients. 

Note that these observations imply that by using a finite number of Matsubara fre- 
quencies it will not be possible to obtain the 4d-vacuum solutions. Any finite number 
of such modes will not be able to generate a divergence stronger than mass-like in the 
longitudinal equation to compensate the missing soft interactions with ghosts. Hence, the 
4d-behavior can be established only by infinite summation. In this case the generated 
mass has to diverge even after renormalization at p = to obtain the vacuum solution, 
providing over-screening instead of screening. The possibility that the ghost-gluon-vertex 
changes due to temperature to couple soft ghost modes to soft longitudinal gluons at low 
temperature and not at high temperature is unlikely, since a bare ghost-gluon vertex is 
already sufficient at low temperatures and in the vacuum ES] • 

^The pure soft interactions represented by the kernels A^i and N2 only generate a constant term in the 
infrared. Without a tree-level mass, the soft tadpoles generate a contribution which can be removed by 
renormalization. 



^(0,go,<?,p) 



(5.7) 
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5.4 Truncation and Spurious Divergences 

Also for large but still finite T the 3d-limit is to be made explicit. Thus all the problems 
encountered in the 3d-theory persist, including the necessity of a modified soft 3-gluon 
vertex. Hence all modifications of the pure soft terms will be left as in the case of chapter 
m except for the tadpoles in the Higgs equation ()4.10|) . Due to the spontaneously gener- 
ated mass, it will be necessary to alter this behavior. There are also additional spurious 
divergences due to the hard mode contributions. As they occur at mid-momenta, when 
viewed from the 4d-perspective using = p^+p^, it is expected that due to the truncation 
ansatz, they are much harder to compensate. This is indeed the case. 

The spurious divergences in contributions of hard modes due to the integration are dealt 
with in exactly the same way as in 3d by adjusting the tadpole terms. Each integration 
kernel K is split as 

K = Ko + Ko (5.8) 

Ko is finite and Kd divergent upon integration. Each kernel Ko will be compensated by 
a corresponding tadpole as for Mtd and M^d in the 3d-theory in equations ()4.38|) and 
()4.39|) . respectively. 

However, a new kind of spurious divergences appear when performing the Matsubara 
sum. In principle this is no problem when using a finite number of Matsubara frequencies, 
since then all expressions are finite. However, when including more and more Matsubara 
frequencies, it is found that their contribution in the gluon equations ()5.5|1 and ()5.6|) at 
\p\ < max(go) scales as go, thus behaving as a quadratic divergence. This is an artifact of 
using a finite number of Matsubara frequencies. It does not vanish as long as the number 
of frequencies is finite, no matter how large the number. This behavior is spurious and 
must be removed. 

The ghost equation does not contain spurious divergences^ and can thus serve as a 
starting point to trace the origin of those in the gluon equations. Analyzing the |p| = 
limit for po 7^ 0, it is directly found that if the equation looks like 

then the integral I must scale like 1/po because of dimensional consistency. Here 
is the wave-function renormalization of the ghost. For any finite number of Matsubara 
frequencies, is finite. Hence the infrared behavior of the hard ghost dressing function is 

■^The problem is also not present in the mixed contributions Mi, M2, Ni and N2 in the gluon equations. 
Hence these do not need to be subtracted. 
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1/(^3 + aT/po) with a suitable constant a. For sufficiently large po this can be expanded 
to yield I/Z3 — aT/{poZ^) in leading order. 

The effect is most directly seen in the po = longitudinal equation in the ghost-loop. 
For dimensional reasons, in the infrared the terms in the Matsubara sum must scale as 
Tgo/|p|^^g('?o)- Hence the divergence is directly visible due to the leading constant term in 
Ag{qQ). The situation is analogous, but not as directly visible in the remaining equations. 

This spurious divergence can be removed when removing the leading I/Z3 term in the 
ghost dressing function, i.e. by the replacement 

G{q,qo)G{p + q,qo+po) 
^ (Giq,qo)-j^) [G{q + p,qo+Po)-j;), (5.10) 

and correspondingly for the gluon loops. 

Including a finite number of Matsubara frequencies only is equivalent to studying a 
3d-theory with 1 massless particle and 2N massive particles for each 4d particle species. 
Therefore the ultraviolet asymptotic analysis of subsection 14.2.11 applies, and I in ()5.9|1 
will again fall off polynomially in the ultraviolet. Hence ()5.10p amounts to subtracting 
the asymptotic constant to which G evolves and thus corrects for being finite. 1/^3 
must become when including an infinite number of Matsubara frequencies, as this is the 
correct perturbative 4d behavior for p ^ T. Hence the spurious divergences are in this 
case only an artifact of the finite number of Matsubara frequencies. 

In the the longitudinal equation the Brown-Pennington value ^ = is pathologic. It 
removes all terms in the soft equation and so is not applicable on its own. Thus the 
replacement ()5.10p has to be applied to the complete loop contributions P, Nt and N^. 
In the transverse equations all spurious divergences including the one of the Matsubara 
sum are again removed at C = 3. Therefore the replacement ()5.10|1 will only be necessary 
in those contributions in R, Mt, and Ml which are proportional to (C — 3). Hence the 
integral kernels are split differently than ()5.8|) as 

similar to ()4.32j) . Here Kq and are finite and independent of (, and Kd contains all 
divergences upon integration. The corresponding subtraction is then performed by the 
replacement 

KD{q)D{p + q)^ K,D{q)D{p + g) + (C - 3)7^3 [D{q) - ^-^ [d{p + " ^) 
where D is a generic dressing function and Z3 its wave-function renormalization. 
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In the soft transverse equation it is additionally necessary to remove tadpole-like struc- 
tures at 7^ 3 in K^^. They contribute a logarithmically divergent mass-term, which is 
irrelevant in the infrared and in the ultraviolet. It cannot be renormalized, as a counter- 
term would be required, which is forbidden by gauge symmetry. Thus these contributions 
are absorbed by the tadpoles as well. Then, in the transverse equation for the soft-soft 
interactions, the subtractions are as in chapter |3 The soft-hard contributions, i?. My, 
and M/) are subtracted as 

KD{q)D{p + q)^ K^D{q)D{p + q) 
+ (C-3)(ir3-^(limp^op'i^3)) [D{q)-ir) [D{p + q)-^^. (5.11) 

5.5 Small- Momentum Approximation and Renormal- 
ization 

Performing the subtractions ()5.1()|1 and 1)5.1111 . the resulting equations are still logarith- 
mically divergent for oo. This also includes the ghost equations. The integrals 
themselves are convergent quantities, but the Matsubara sum is not, as the terms scale 
like 1/go- Thus, for an infinite number of Matsubara frequencies, the sum diverges. This 
is the 4d logarithmic divergence and thus the usual one of Yang-Mills theory pp. 

In the case of a finite number of Matsubara frequencies, the sum is finite. Nonetheless, 
there are a few aspects to be dealt with in the following subsections. 

5.5.1 Small-Momentum Approximation 

There are three regions of momentum to be distinguished: 

At p < 27rT, all Matsubara sums behave essentially the same due to the Appelquist- 
Carrazone decoupling theorem. It is valid as the hard modes behave essentially tree-level- 
like. As the contributions of the hard modes behave as l/go? the final result will depend 
on the number of Matsubara frequencies included, and by adding an arbitrary number, 
any quantitative (but not qualitative) result can be generated. Thus, it is necessary to 
renormalize in order to be independent of the cutoff. The cutoff is in this case imposed 
by the number of Matsubara frequencies included. The renormalization procedure will be 
implemented in the next subsections. By this approach, the results can be made quite 
reliable in this regime. 

At 27rT < p < 2ttT{N — 1), where is the number of Matsubara frequencies included, 
more and more Matsubara terms depart from their 1/qo behavior to a 1/p behavior. As 
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the external momentum p becomes large compared to the effective mass go of a hard mode, 
the mode becomes dynamical and behaves like a massive 3d-particle. The results are still 
quite reliable in this region as also in the case N —>■ oo, for any finite momentum p only 
a finite number of Matsubara frequencies are dynamical. It becomes less and less reliable 
when approaching the upper limit 2tiTN . 

At p > 2ttT{N — 1), the situation changes drastically. Opposite to the case N = oo, 
all Matsubara modes are dynamical and their contributions will scale as 1/p. Thus the 
number of Matsubara modes will now enter linearly instead of logarithmically. This is 
an artifact of cutting off the Matsubara sum. In this region the results are not reliable. 
However, as the sum is still finite and suppressed by 1/p, the contribution is subleading 
with respect to the tree-level term and thus the system of equations can still be closed 
consistently^. 

By renormalization, these artifacts can be reduced, if not completely removed at suf- 
ficiently small momenta. In that sense the finite Matsubara sum approximation is indeed 
a small-(3-)momentum approximation. 

5.5.2 Renormalization of a Truncated Matsubara Sum 

There are a few subtleties involved concerning the renormalization of a truncated Mat- 
subara sum. These are discussed here. 

Firstly, in the vacuum case ^HEHl, renormalization was performed using a momentum 
subtraction scheme (MOM), i.e. by subtracting from the DSEs themselves at a fixed 
subtraction point s. This leads to a generic form of 



where D is the corresponding dressing function and / is the self-energy contribution. As 
the highest divergences encountered are logarithmic P , ()5.12j) is finite and can be solved if 
a proper value of D{s) is known. This can be achieved by choosing s inside the perturbative 
region. This approach fails if the large momentum asymptotic value of D is a constant 
different from 0, since ()5.12j) is ambiguous with respect to such a constant. In the case 
of a truncated Matsubara sum, the ultraviolet behavior is that of a massive 3d-theory. 
Thus the self-energy contributions vanish as 1/p in the ultraviolet, see subsection 14.2.11 
Therefore the dressing functions D are dominated by the tree-level term 




(5.12) 



'^Otherwise the corresponding finite 3d-theory would be ill-defined, which is not the case. 




(5.13) 
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with Z3 their wave-function renormahzation. For a finite number of Matsubara frequencies 
Z3 is finite and thus D goes to a constant in the ultraviolet and ()5.12|1 is not applicable. 
Therefore the renormahzation is performed by explicit counter-terms. This is discussed in 
the next subsection. 

A second point is the mass and mass renormahzation necessary for the soft longitudinal 
mode. The soft mode with frequency po = of the Aq component of the gauge field 
transforms homogeneously instead of inhomogeneously under gauge transformations ()2.5p . 
Therefore gauge symmetry permits to add a term 



to the Lagrangian. In the vacuum, however, such a term is forbidden by manifest Lorentz 
invariance. At finite temperature in the Matsubara formalism, this is no longer the case, 
and such a term could in principle be present. This term replaces the p^A^ term of the 
hard modes, which stems from the AqOqAq term in the Lagrangian ()2.2|) for po 7^ 0. 

Concerning the counter-terms, the wave-function renormahzation is performed by adding 
the counter-term 



to the Lagrangian. In the case po = 0, the first term is not present for the longitudinal 
mode Aq. Its place can be taken by a counter-term for the mass term ()5.14|) . This implies 
a relation between the wave- function counter-term 5Z^ and the mass counter-term 5m?, 
which cannot be exploited here due to the truncation of the Matsubara sum. Therefore 
an independent renormahzation of the wave function and the mass of the soft longitudinal 
mode is necessary and will be performed in the next subsection. 

A last point concerns the implications for the counter-terms due to the truncation of 
the Matsubara sum. As the divergence structure must be the same as in the vacuum, 
the counter-terms must be the same for all frequencies. This is no longer the case when 
the Matsubara sum is truncated, as can be seen directly by counting. In the sum for 
the soft mode, po = 0, contributions from 2A^ — 1 Matsubara modes are present for each 
loop. For the hard mode with po = 27rT(A^ — 1), only contributions are present, as 
bo + Q'ol ^ 2'n-{N — 1)T. Thus different numbers of modes contribute and the counter- 
terms cannot be the same. This is always the C3iSG clS long as A^ < 00. To surpass the 
problem in a constructive manner, each mode will be renormalized independently. 

Therefore, a counter-term Lagrangian is added, given by 



m'Al{0,p) 



(5.14) 



5Zs{A;d,d.A^^ - A^d'A;) 




10 



+ 




(5.15) 
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where A'^{qo) are the modes of the gluon field and c(go) and c(go) are the modes of the 
ghost and anti-ghost field, respectively. At/l/mu are the appropriate tensor structures 
of derivatives. Consequently, the corresponding wave-function renormalization constants 
in the subtractions ()5.10|) and 1)5.111) have to be replaced by their frequency-dependent 
versions. 

This shortcoming is a consequence of the small momentum approximation yielding the 
incorrect ultraviolet properties: The treated system is for large momenta equivalent to a 
3d-theory of 2N — 1 particles per 4d particle species. In such a theory all fields can be 
renormalized independently. This problem has to be investigated further to obtain reliable 
results on the large momentum behavior of the dressing functions. 



5.5.3 Implementation of Renormalization 

In all equations but the one for H{0,p), using the counter-terms of ()5.15|) amounts to 
replacing the tree- level term 1 by 1 -|- This way explicitly multiplicative renor- 

malization is obtained. In the equation for H{0,p), the tree-level term is replaced^ by 
1 -|- 6Z3L + 5m? /p^. Thus, it also generates a mass-renormalization. By the arguments 
following equation ()5.7|) just a renormalization is indeed necessary, and the mass is not 
generated purely by this procedure. 

The last ingredient is the renormalization prescription. To investigate both cases dis- 
cussed in section 15. 2| either a fixed subtraction point s = sq or to establish a well-defined 
3d-limit s = T is employed. Thus, the following prescription is applied to all dressing 
functions D, except H{0,p): 

D{s) = 1. 

This also ensures G'(s)^Z(s) = 1 as is required in the 4d-theory [SS]- The soft ghost mode 
is not renormalized at as in the vacuum calculation pHI. 

In case of H{0,p), two prescriptions are necessary. The first is^ 

^ 1 1 



lim \p\-'H{0,p) - „ 



where is the tadpole-improved mass of the 3d-theory. It depends on r = ruh/g^, 
which is again taken to be the same value as in the previous chapter. For fixed (73, this 
mass is independent of the temperature while for fixed /i it scales with temperature as 



^Note that when using the projector 1)5. 2|l . also Sm"^ is multiplied by 

^For numerical reasons, it is actually performed at 2Si, where Si is the numerical IR-cutoff of the 
integration. This is only a marginal difference. 
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expected. This again indicates the necessity of control over the infinite-temperature hmit 

of the theory in terms of the effective couphng constant. In addition, 

1 ^ . msd 
H{0, s)~ s2 

is required to fix the mass renormahzation. 

It should be noted that, by renormalizing at T, it is guaranteed that the correct 3d-limit 

is obtained. This prescription requires that atT^oo all dressing functions approach 1 at 

infinity. This yields the 3d-results^ as in chapter HI It would be easily possible to readjust 

the requirements on H{0,p) to comply with the lattice data in the infrared and ultraviolet. 

Here the exact correspondence to the 3d results in chapter 0] is preferred, as this will make 

the 3d-limit explicit. 

Therefore the explicit implementation of the renormahzation prescription for the DSE 
of a dressing function D with self-energy contributions / 

1 



is then 



and in the case of H{0,p) 



Dip) -'^'^'^ 
<5Zs = -/(s) 



Sm^ = ml — \irap^I(p) 
limp_op'^I{p) 



6ZsL = -lis) 



S2 



where mr = m^d is the renormalized mass. 



5.6 Ghost-Loop-Only Truncation 

The first step is again the calculation of the ghost-loop-only solutions. Performing the 

calculations turns out to be numerically more complicated than in 3d due to the recursive 

definition of the elimination of the spurious divergences and ()5.11|) . This is discussed 

in more detail in appendix iFl Nonetheless, both solutions can be found. They are shown 

in figure 15.11 with the largest number of Matsubara frequencies currently available^ at 

^This is only correct, if g|T <C T. This is immanent, if g\T is fixed. If /i is fixed, this is only possible 
if /i is taken sufficiently large for a fixed temperature. Still deviations from 1 will be found for \p\ 3> /i, T . 
^The number of frequencies included is only a matter of CPU-time. 
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Figure 5.1: The left panels show for the soft modes from top to bottom the dressing 
functions of the ghost and transverse gluon and the propagator of the longitudinal gluon. 
The right panels show the hard mode dressing functions. The less the deviation from 1, 
the harder the mode is. The sohd hues are the g — 1/2 solution and the dashed hues the 
g ~ 0.4 solution. Both are at T = 1.5 GeV and N — 20. Note the different momentum 
scale in the left and the right panels. 
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( = ^ = 1. In general, all hard mode dressing functions are similar, but behave more and 
more tree-level like with increasing frequency po = 27mT. The different asymptotic values 
of the dressing functions for each mode are a direct manifestation of the artifacts of the 
truncation of the Matsubara sum, as discussed in section 15.5.21 

The small deviation from tree-level for the hard longitudinal modes is most likely due 
to the severe subtraction ()5.1()|1 . Within the current ansatz however, this is the best 
achievable. Note that due to = Po + P^^ infrared of the hard modes is actually 
mid-momenta with respect to 4d. Thus, these problems do not come unexpected, as 
the truncation is tailored to the infrared and thus to the soft modes. In addition, the 
transverse gluon wave- function renormalization is smaller than 1 for n = 0. Thus the 
dressing function would diverge for large momenta in the limit N ^ oo. This is also an 
artifact of the truncation, as the ghost loops alone have the wrong sign at large momenta. 
This will be remedied by adding the gluon loops in the next section, and can be ignored 
for now. 

The dependence on the number of Matsubara frequencies is shown in figure 15.21 As 
expected, in the infrared region, the soft mode dressing functions are not affected by the 
number of Matsubara frequencies. The onset of deviation of the hard mode dressing func- 
tions calculated with N frequencies from the solutions with M > N frequencies increases 
in 3-momenta with A^. This is a manifestation of the small momentum approximation. 

Naturally, the dependence on the temperature is the next subject. Using the t'Hooft- 
like scaling, the results shown in figure 15.31 are obtained. While the soft mode dressing 
functions are relatively unaffected by temperature, significant changes are found for the 
hard mode dressing functions. These become increasingly tree-level with increasing tem- 
perature and thus the hard modes cease to interact. At low temperature they exhibit 
significantly more dynamics. 

The case for fixed fi is shown in figure 15.41 The soft longitudinal gluon propagator 
is scaled by the interaction strength, since the mass of the soft longitudinal gluon is 
a dynamical effect. The results represent, how the system reacts on a change of s/T, 
which are the only dimensionful quantities involved as long as (74 is fixed. The changes 
are governed by the renormalization prescription, and are qualitatively different from the 
previous case. As expected, the higher the temperature, the further out in momentum 
the non-perturbative effects in the soft mode dressing functions extend, only limited by 
the requirement to obey the renormalization prescription. At the same time the dynamics 
of the hard modes is pushed to higher momenta, although the dressing functions do not 
decrease. Indeed, this scenario describes rather a change of the relevant scale for the onset 
of non-perturbative effects than the effects of temperature. 
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Figure 5.2: The dependence of the solutions on at T = 1.5 GeV. The left panels show 
for the soft modes from top to bottom the dressing functions of the ghost and transverse 
gluon and the propagator of the longitudinal gluon. The right panels show the dressing 
functions for the n = 1 hard mode. The red lines show the solution g = 1/2 and the black 
lines g ~ 0.4. Solid is = 2, dashed is = 5, dotted is = 12 and dashed-dotted is 
A^ = 20. 
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Figure 5.3: The dependence of the solutions on T at fixed gfT and subtraction at T at 
= 5. The left panels show for the soft modes from top to bottom the dressing functions 
of the ghost and transverse gluon and the propagator of the longitudinal gluon. The right 
panels show the dressing functions for the n = 1 hard mode. The red and black lines 
are the g = 1/2 and the g ^ 0.4 solution, respectively. Solid is T = 0.8 GeV, dashed is 
T = 1.5 GeV, and dotted is T = 137 GeV. 
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5.7. Full Theory 



At fixed fi, it was possible^ to follow tlie solution down to T = 91 MeV. Lattice 
calculations find a phase transition temperature of (269 ± 1) MeV [9,- Thus super-cooling 
is possible, provided the temperature estimate of section 15.21 is correct within a factor 
of 2 and the phase transition temperature in the approximation scheme used here is not 
significantly lower than in lattice calculations. 

The remaining issue is the dependence of the solutions on the projection and thus 
on ( and ^. It is shown in figure 15.51 The largest effect is seen in the infrared region 
of the soft solutions of the alternating branch when varying (, similar to the infinite- 
temperature limit. The hard modes are significantly less affected, with the largest impact 
on the hard longitudinal dressing functions when varying ^. For the smallest ^ values, the 
self-energy of the hard longitudinal mode seems to switch sign, generating a wave-function 
renormalization smaller than 1. This effect has not been observed in the full system and 
is thus likely an artifact of the ghost-loop-only truncation. Thus, the dependence on the 
projection is of similar extent as in the infinite-temperature limit. 

5.7 Full Theory 

Finally, the full system of equations ()5.4ll5.(i|l is implemented. For the largest available 
number of Matsubara frequencies, the result is shown in figure 15.61 There are several 
observations. First of all, the soft modes are again nearly unaffected in the infrared by 
the presence of the hard modes. The latter show a significant modification, compared to 
tree-level, although still only of the order of 30%. Also, all wave-function renormalization 
constants are now larger than one, as is required. The hard mode dressing functions exhibit 
some structure. There are maxima in all dressing functions. This is most pronounced in the 
case of the dressing functions of the longitudinal gluon. These structures do not translate 
into a corresponding structure in the propagators, which are monotonically decreasing 
from a constant of order 1/pq in the infrared to in the ultraviolet. It is also nearly 
irrelevant for the hard mode dressing functions to which soft infrared solution they are 
coupled to. 

Figure 15.71 displays the dependence of the full solutions on the number of Matsubara 
frequencies. While the soft mode dressing functions are nearly unaffected, apart from the 
value of the renormalization constants, the effect on the hard mode dressing functions is 
significant. The ghost is quite insensitive, except for its wave-function renormalization. 
This is not the case for the gluons. The peaks at mid-momentum are sensitive to the 
number of Matsubara frequencies included. The effect is largest for the longitudinal gluon 

^In the case of fixed g^T, it is likely not yet found due to limitation in CPU time. 
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Figure 5.4: The dependence of the solutions on T at fixed /i and s at = 5. The left 
panels show for the soft modes from top to bottom the dressing functions of the ghost 
and transverse gluon and the propagator of the longitudinal gluon. The right panels show 
the dressing functions for the n = 1 hard mode. The red and black lines are the g = 1/2 
and the g ~ 0.4 solution, respectively. Solid is T = 0.3 GeV, dashed is T = 1.5 GeV, and 
dotted is T = 5.0 GeV. 
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Figure 5.5: The dependence of the solutions on ( and ^ at = 5. The left panels show 
for the soft modes from top to bottom the dressing functions of the ghost and transverse 
gluon and the propagator of the longitudinal gluon. The right panels show the dressing 
functions for the n = 1 hard mode. The red and black lines are the g = 1/2 and the 
g ^ 0.4 solution, respectively. Solid is (gX^O = (1/2,7/2,9/4) and (0.4,2,5/2), dashed 
is (1/2, 1, 1) and (0.4, 1, 1), dotted is (1/2, 0, 0.1) and (0.4, 1/4, 0.1). 
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Figure 5.7: The dependence of the solutions on at T = 1.5 GeV. The left panels show 
for the soft modes from top to bottom the dressing functions of the ghost and transverse 
gluon and the propagator of the longitudinal gluon. The right panels show the dressing 
functions for the n = 1 hard mode. The red lines show the solution g = 1/2 and the black 
lines g ~ 0.4. Solid is = 2, dashed is = 5, dotted is = 8 and dashed-dotted is 
A^ = 12. 
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dressing functions. From the available number of Matsubara frequencies, it is hard to 
estimate whether the peak grows to a finite value when — > oo or not. It cannot be 
excluded that the hard longitudinal gluon dressing functions violate the Gribov condition 
()3.10|) once sufficiently many Matsubara frequencies are included. This would be very 
similar to the case of the transverse gluon dressing function in the infinite-temperature 
limit and therefore would necessitate a similar vertex construction for the longitudinal- 
gluon-transverse-gluon vertex as for the soft transverse gluon vertex. If such a vertex 
would be necessary, the result would be a finite peak, very similar to the present situation. 
Thus, the results would be even quantitatively quite similar. 

The disappearance of the peaks in the gluon dressing functions at = 2 can be directly 
related to the vanishing of hard-mode couplings due to the restriction {po + go) / (2vrT) < 2. 
Only hard-soft-mode couplings contribute. Thus, the peak is generated due to pure hard- 
hard interactions alone. 

Again here and in the following only the n = 1 modes are presented. The n > 1 modes 
essentially follow their behavior, albeit much closer at tree-level. 

The dependence on temperature when using t'Hooft-like scaling is shown in figure EUl 
Here, a significant effect also on the soft modes can be seen. While a direct influence on the 
infrared coefficients for the ghost and transverse gluon dressing functions is visible, for the 
dressing function of the longitudinal gluon only a slight change at mid-momenta occurs. 
Therefore, the dressing function of the longitudinal gluon is dominated by its renormalized 
mass also when changing the temperature. As expected, the hard mode dressing functions 
become more and more tree-level like when increasing the temperature. The peaks in 
the hard mode dressing functions shift to higher momenta, owing to the renormalization 
condition, and become smaller due to the increase of the effective mass po of the hard 
modes. At low temperature, the opposite effect is observed. In general, the hard mode 
dressing functions become more dynamical as their effective mass decreases, albeit quite 
slowly. The insensitivity on the infrared solution of the soft sector is not changed when 
reducing the temperature, and it may need a significantly lower temperature to induce a 
change. 

Switching to a fixed renormalization scale and subtraction point, the solutions shown 
in figure 15.91 are found. The effect is quite the same as in the ghost-loop-only truncation. 
Changing the temperature at fixed 4d-coupling, and thus increasing the 3d-coupling with 
temperature, shifts the onset of non-perturbative effects to larger momenta. Obeying the 
renormalization at fixed momentum then requires the changes in the dressing functions 
seen in the figure. In this case, it was also again possible to super-cool down to at least 
T = 152 MeV. 
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Figure 5.8: The dependence of the solutions on T at fixed gfT and subtraction at T at 
= 5. The left panel shows the soft modes, where apart from the longitudinal gluon the 
dressing functions are shown. In the latter case the propagator is shown. The right panels 
show the dressing functions for the n = 1 hard mode. The red and black lines are the 
(7 = 1/2 and the g ^ 0.4 solution, respectively. Solid is T = 0.9 GeV, dashed is T = 1.5 
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Figure 5.9: The dependence of the solutions on T at fixed and s at = 5. The left panel 
shows the soft modes, where apart from the longitudinal gluon the dressing functions are 
shown. In the latter case the propagator is shown. The right panels show the dressing 
functions for the n = 1 hard mode. The red and black lines are the g = 1/2 and the 
g ^ 0.4 solution, respectively. Solid is T = 0.3 GeV, dashed is T = 1.5 GeV, and dotted 
is T = 3.2 GeV for ^ = 1/2 and T = 2.3 GeV for g ^ 0.4. 
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Figure 5.10: The dependence of the solutions on ( and ^ at = 5. The left panel shows 
the soft modes, where apart from the longitudinal gluon the dressing functions are shown. 
In the latter case the propagator is shown. The right panels show the dressing functions 
for the n = 1 hard mode. The red and black lines are the g = 1/2 and the g ~ 0.4 
solution, respectively. Sohd is {g,CyO = (1/2)3,2) and (0.4,1.7,7/4), dashed is (1/2,1,1) 
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The dependence on ( and C, is shown in figure EH Varying ( leads to similar variations 
for the soft mode dressing functions as in the 3d-case. The largest effect is seen in the 
infrared. As the hard mode dressing functions are insensitive to the infrared behavior of 
the soft modes due to their effective mass, there are only weak variations of them with (. 
Only at small ( an additional structure appears in the transverse gluon dressing function. 
This is likely due to the cross-term in equation ()5.fjj) . The effect of varying ( on the 
longitudinal sector is negligible, as in the 3d-case. Correspondingly, the effect of varying 
^ is only significant for the longitudinal sector. The sensitivity is much less pronounced 
then in the case of varying (. Especially the soft longitudinal dressing function is nearly 
unaffected by variation of ^. The reason is that the soft equation does not depend explicitly 
on ^, as it can be divided out of the equation for ^ 7^ 0. Therefore, any ^ dependence only 
enters indirectly by the weak dependence of the hard modes on ^. 

The dependence on 6 is qualitatively not different from the high temperature case for 
the soft modes. With decreasing temperature, the lowest attainable value of 6 increases 
and is around 6 ~ 0.11 for the g = 1/2 solution and 6 ~ 0.15 for the g ^ OA solution at 
T = 1.5 GeV. Surprisingly, the roles of the g = 1/2 and g ~ 0.4 solutions are interchanged: 
At low temperatures, it is possible to reach lower values of 6 in case of the half-integer 
solution. The hard mode dressing functions are nearly not affected by the change in 6. 
Thus, there is no significant change concerning the dependence on 6 at finite temperatures. 

At this point, a comparison can be made to a different approach to obtain the high- 
temperature gauge propagators. The usual continuum method is the semi-perturbative 
hard thermal loop (HTL) approach [01]. It is based on resumming the hard mode con- 
tributions in self-energy diagrams, hence its name. In the transverse infrared sector it is 
plagued by severe problems, due to its perturbative nature. The final result is a transverse 
gluon propagator with a particle like pole at p = 0, thus t = 0. This is in sharp contrast 
to the results found here and, as discussed in chapter El such a behavior is not likely. Also 
the lattice results tend to support the results found here. Concerning the soft longitudinal 
mode and the hard modes, HTLs and the ansatz presented here find qualitatively similar 
results on the level of the propagators. 

5.8 Solutions at Small Temperature 

As the temperature region of and above the phase transition is probed by experiments, 
the aim of the calculations presented here is to eventually describe the phase transition. 
This is not yet possible. Currently, the phase transition is approached from two sides. The 
first approach is from above, and is presented in this work. The second approach is from 



84 



5.8. Solutions at Small Temperature 



Ghost dressing function 




T=0 MeV 
T=100 MeV 
T=200 MeV 



p[GeV] 



Figure 5.11: The dressing function of the soft ghost mode at zero |^ and finite tempera- 
ture inn!- 

below and treated especially in and also to some extent in the context of this work 
[nSl- The main aim in these calculations is to investigate whether any qualitative changes 
are observed at small temperatures as compared to the vacuum solutions. Sufficient for 
this aim is to employ not a full continuum method but solving the DSEs approximately on 
a discretized space-time with a toroidal topology. This method had already been applied 
successfully to the vacuum [HHl HO] • 

The results for the soft mode dressing functions are shown in figure IKTTl for the ghost 
propagator and in figure 15.12! for the two independent dressing functions of the gluon 
propagator ()3.5|) . compared to the zero temperature solution. While the 3d-transverse 
part of the propagator becomes steeper in the infrared, the 3d-longitudinal part becomes 
more shallow. Thus, it experiences significant finite volume effects, as in the case of lattice 
calculations. Finite volume scaling shows that it is still confined, as in the case of the 
vacuum (HH]- Indeed it cannot be distinguished yet whether the exponents or merely the 
coefficients of the infrared solution change. 

Interestingly, the 4d-transverse projection of ()3.5|1 . yields that the quantity 2Z + H, is 
extremely temperature-independent [HS] • As this combination is the only one entering the 
thermodynamic potential to be discussed in section 16. 2[ and G in itself is rather tempera- 
ture independent, the corresponding thermodynamic potential is itself nearly independent 
of temperature and compatible with 0. This is what would be naively expected from a 



Chapter 5. Finite- Temperature Effects 



85 




Figure 5.12: The propagator of the soft gluon mode at zero jHEI and finite temperature 
P^ . The left panel shows the 3d-transverse part Z of ()3.5p and the right panel the 3d- 
longitudinal part H. 

confined gluon system, as bound states are not included in the description here. 

In addition, it is possible to trace the solution up to temperatures of 800 MeV [03], 
significantly above the expected phase transition temperature of (269 ± 1) MeV according 
to lattice results p. Therefore it is possible to super-heat the system. 

If the temperature scale used in [Hni and here are compatible within a factor of 4, then 
within a substantial temperature range, two qualitatively different solutions exist. One 
solution has an over-screened and one a screened soft longitudinal gluon. Together with 
the observation of super-heating and super-cooling when comparing to lattice results, this 
implies a first order phase transition. 

Compared to the high-temperature case presented before, the sector of ghosts and 
transverse gluons is not changed qualitatively. Quantitatively, especially the infrared ex- 
ponents changed, though. The main difference is in the longitudinal sector, where a change 
from over-screening to screening occurred. Comparing the hard modes in both cases, no 
qualitative change is found in any sector. Therefore, the main difference of the two phases 
is only in the longitudinal/chromoelectric sector. Therefore the nature of the phase tran- 
sition will likely be chromoelectric. This implication for the phase transition and the 
compatibility of these observations with lattice calculations will be discussed in section 



Chapter 6 
Derived Quantities 



The dressing functions obtained in chapters ID and El are only the first step from the basic 
Lagrangian towards observables. In this chapter, derivations towards such observables are 
performed. The screening masses and with them the analytic structure of the particles 
described will be investigated in sect ion IfTTl A second natural object is the thermodynamic 
potential discussed in section IFT^ 

If not stated otherwise, throughout this chapter Ca = Nc = 3, ( = ^ = 1, 5 = 1/4, 
and the lattice value for myj are used. 

6.1 Schwinger Functions and Analytic Structure 
6.1.1 Schwinger Functions 

Screening masses are most directly extracted from the analytic structure of the propaga- 
tors. To obtain access to these analytic properties, the Schwinger function related to the 
dressing function D is calculated. It is defined as [31] 



i.e. the Fourier transform of the propagator with respect to (Euclidean) time. Note that 
this definition is independent of the dimensionality of the underlying theory. Negative 
values for the Schwinger function can be traced to violations of positivity and therefore to 
absence of the particle represented by D from the physical spectrum [5^ . 

The numerical implementation of (jfj.lj) is non-trivial. To obtain sufficient precision, a 
FFT-algorithm using edge corrections and compensation for the infinite integration range 



oo 




(6.1) 
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Figure 6.1: The left panel shows the Schwinger function of the gluon. The right panel 
shows the comparison of the fit-function ()6.2j) to the full solutions of the dressing function. 
The solid line gives the numerical result for g = 1/2 while the dashed line is for g ^ 0.4. 
The dotted and dash-dotted line denote their respective fits using the ansatz ()6.2|) . 

of ()6.H) together with at least 512 or more frequencies^ was employed 

The time t in ()6.1|) in the 3d theory of the infinite temperature limit is actually a 
space-direction in 4d. To such a direction also the arguments concerning positivity can be 
applied in Euclidean space-time. Therefore, the result for the gluon, presented in figure IFTTI 
clearly exhibits positivity violations. This is in accordance with the Oehme-Zimmermann 
super-convergence relation ()2.25|) and thus has been expected. Furthermore, the position 
of the zeros can be interpreted as the confinement scale. For the g = 1/2 solution, the zero 
occurs at g^t ^ 3.29 and at g^t ~ 3.98 for the other branch. This result is in agreement 
with recent lattice results 

To be able to perform an analytic continuation of the gluon propagator into the complex 
p^-plane, the Schwinger function is fitted. This is performed by using the ansatz^ 



1 + / + /l,p4<;+i ^ 1 + jauV) (1 + a„/p)(l + / + A.p^^+i) 
for the dressing function. is the infrared coefficient determined previously, and is 



"'^Ten thousands to a million are much better. The results presented here have been obtained using 
roughly 5 • 10^ frequencies. 

^In |54| . a parameterization with only a branch cut and no isolated pole provided a successful fit. Due 
to the different asymptotic behavior in 4 and 3 dimensions, a different ansatz is used here. 
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Solution 




au9i 


/ 


5 = 1/2 


20.3 


64/27 


1.32511 


g « 0.4 


13.4 


64/27 


1.03148 



Table 6.1: The coefficients for the gluon fit ()6.2|) . 



the ultraviolet coefficient of leading-order resummed perturbation theory, as calculated in 
appendix|El The fit parameters for both solutions are given in table lU^Tl As demonstrated 
in figure lU^Tl the Schwinger function is fitted very well. The gluon propagator and dressing 
function are also reasonably well described by the fit. 

Similar to a meromorphic fit ansatz used in [31] , the Higgs propagator and its Schwinger 
function, the latter being analytically calculated from the former ^7\, are described as 

^^^^^ ^ p^ + 2w? cos (20)^2 + ^4' (6.3) 

^/W = TT^^TTTTVe"*"'"^'^^ fsin (0 + tmsin (0)) + ^ sin (0 - tmsin (0))^ . 

im-^ sm (20) \ e / 

As demonstrated in figure these fits already describe the Schwinger function quite well, 
but miss around 20% of the propagator at zero momentum. This indicates that further 
massive modes are present. Indeed, for the g = 1/2 solution, a further term of the form 
fl6.3j) has to be added. For the g ^ 0.4 solution, adding a term with one pole, 

= (6-4) 
A/W = (6.6) 

improves the fit also in this case. Both subleading fits are not very accurate for large t, 
and the results have to be taken with care. They indicate, nevertheless, the existence of 
subleading contributions due to the presence of further massive-particle-like contributions. 
The fit parameters of both solutions can be found in table 16.21 

As can be seen from figure 14.161 the result for the Higgs propagator deviates signifi- 
cantly from the lattice results, the self-energy being significantly overestimated. While the 
gluon Schwinger function is reasonably independent of the truncation, this turns out not 
to be the case for the Higgs Schwinger function. In order to show this, the solutions for a 
bare Higgs-gluon vertex suppressed via a scaling factor uj are obtained^. The ffist result is 



■^In the Landau gauge, due to the transversahty of the gluon, this is equivalent to modifying the tensor 
structure of the vertex. 
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Figure 6.2: The Higgs propagator is shown in the top panels and its Schwinger function 
in the middle panels compared to their fits. The bottom panels show the comparison of 
the numerical (sohd) and the fitted (dashed) subleading Higgs contribution, see text. The 
left and right side shows the g — 1/2 and g 0.4 solutions, respectively. The solid lines 
represent the numerical solution, the dashed lines give the leading contribution and the 
dashed-dotted lines the first subleading contribution. The dotted hues underneath the 
solid lines give the sum of the leading and subleading contribution. 
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Figure 6.3: The bottom panel shows the Higgs propagator for various suppression factors 
a; of the Higgs-gluon vertex. The top panel shows the corresponding Schwinger functions. 
The left panels contain the g = 1/2 solution while the right ones display the g ^ 0.4 
solution. The lattice data are the same as on the right side of figure 14.161 The solid 
line is for = 1/4, the dashed line for = 1/2, the dotted line is for a; = 3/4 and the 
dashed-dotted hne for = 1. 
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Solution 


e 


fgl 


4> 




^7 = 1/2 


4.0199 


0.3545 


-0.67078 


1.4998 


subleading 


9.4493 


0.6736 


-0.18387 


2.561 


g ^ 0.4 


4.0697 


0.37793 


-0.63975 


1.5045 


subleading 


0.81188 






1.9223 



Table 6.2: The coefficients for the Higgs fit ()6.3j) and the subleading one ()6.3|l and ()6.4j] . 
respectively. 

that the gluon Schwinger function, in contrast to the Higgs one, is not susceptible to such 
a change. 

Figure IFT^ displays corresponding results for different values of u. The position of the 
first zero tends to increase for decreasing uj. At ^ 1/4 any oscillation, at least within 
the available numerical precision, seems to be gone altogether. A fit to the Schwinger 
function using ()6.5p reveals again additional structure at very small t which cannot be 
captured by such a simple fit. As the fit also misses some strength at zero momentum for 
the propagator, this again indicates the presence of further massive contributions in the 
propagator. 

When adding the hard modes, the qualitative picture does not change, as is visible 
in figure 16.41 In this case the Schwinger function also depends on a 'space '-variable with 
respect to 4d and Pq = 2iTTn acts as an effective mass. The gluon soft mode shows an 
oscillation as expected. The same is still true for the Higgs, thus this property of the 
truncation used here is still present. At this rather low temperature, it exhibits only one 
oscillation, and is thus qualitatively much more similar to the gluon. In both cases, the 
oscillation frequency has decreased due to the additional interactions with the hard modes. 
The interpretation of the oscillations in the longitudinal gluon Schwinger function are still 
inconclusive, though. 

For the hard mode Schwinger functions no oscillations are seen. They cannot be 
excluded as they might probably not be resolvable numerically if the oscillation frequency 
is of the order of the mass. Up to the numerical accuracy, the hard modes present the 
behavior expected for a massive particle. This equally well applies to the higher modes, 
as is visible by the example of the next mode n = 2. 

Thus, the features found in the infinite temperature limit persist at finite temperatures. 
Note, however, that if the hard modes oscillate, this oscillation does not necessarily vanish 
when the T — > oo limit is taken. If the oscillation frequency only vanishes as fast as 
the interaction, then even for arbitrarily small interactions, oscillations remain. It will 
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n=0 transverse gluon Schwinger function 



n?^ transverse gluon Schwinger function 




Figure 6.4: The Schwinger functions at T = 1.5 GeV and = 5. The left top panel shows 
the soft transverse gluon Schwinger function, the right top panel the hard transverse gluon 
Schwinger function for n = 1 with a black line and for n = 2 with a red line. The bottom 
panels show the same for the longitudinal gluon Schwinger function. In all cases are solid 
lines the g = 1/2 solution and dashed lines the g ~ 0.4 solution. 
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take more sophistication to settle the question of whether the hard modes show positivity 
violation or not. 

6.1.2 Analytic Properties in the Infinite- Temperature Limit 

Although the high-temperature limit of the four- dimensional Minkowski theory is a gen- 
uinely Euclidian theory, it will be of interest for other applications to extract the analytic 
structure of the propagators investigated. 

The gluon propagator exhibits similar behavior for both solutions, but there are also 
some significant differences. The denominator of the ansatz ()fi.2j) contains two factors, 
both of which could possibly give rise to a non-trivial analytic structure. The first part 
stems from the fit of the perturbative tail necessary to generate the maximum in the gluon 
dressing function. Since all fit parameters are positive, this factor does not give rise to 
a pole on the first Riemann sheet. However, it generates a pole on the second Riemann 
sheet, which will occur at 1/auf, that is, at Euclidian momenta. This pole does not have 
a physical interpretation, and may well be an artifact of the fit, since ()6.2j) is tailored to 
generate the correct leading-order perturbative behavior. Thus, it generates most likely a 
structure which has the Landau pole of perturbation theory on the second Riemann sheet. 
It is to be expected that this pole vanishes when using a more sophisticated fit. 

The second factor generates a genuine isolated pole at {-{1 + f)/A,)-^/\ This expres- 
sion has only one value on the first Riemann sheet given by (—0.0933 + 0.1615i)g^ for the 
(7 = 1/2 and (—0.1462 + 0.1248^)5^3 for the other solution^. In both cases, a pole close to 
the origin is generated, with an imaginary part larger than the real part in one case. For 
the first solution, the pole is found for an angle significantly above 7r/4 while in the second 
case somewhat below. In addition, the first solution generates two more poles on two 
more Riemann sheets, while the second solution, with a (most likely) irrational exponent, 
generates an infinite number of further poles on an infinite number of Riemann sheets. In 
both cases, the residue is complex. In addition, there is a cut along the complete negative 
real axis starting at zero. In this way it is similar to the results in four dimensions 

This analysis infers that the gluon propagator is violating positivity, and thus satisfies 
the requirements of the Kugo-Ojima and Zwanziger-Gribov confinement scenario. 

On the other hand, the Higgs propagator very likely does not have a branch cut but 
a number of simple poles whose locations are given in table 16.31 The sensitivity to the 
Higgs-gluon vertex, however, necessitates further investigations before a firmer conclusion 
can be drawn. 



This is of the order of Aqcd when using the 't Hooft-like scahng of section 
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Solution 


Order 


Pole/^7| 


^? = l/2 


Leading 


-0.5112 ±2.191z 


^? = l/2 


Subleading 


-6.11972 ±2.35777i 


g ^ 0.4 


Leading 


-0.650078 ± 2. 16822i 


g ^ 0.4 


Subleading 


-3.69539 



Table 6.3: Location of the poles of the Higgs propagator for oo = 1. 



First zero crossing of the IHiggs Schwinger function 




0.5 0.6 0.7 0.8 0.9 1 



Figure 6.5: The dependence of the first zero of the Higgs Schwinger function on the vertex 
suppression u. Solid is the g = 1/2 solution and dotted its fit, dashed is the g ~ 0.4 
solution and dashed-dotted its fit. 

The dependence of the first zero of the Higgs-Schwinger function on u is shown in 
figure 16.51 The result can be fitted using 

a 

to = —r 

where to is the position of the first zero. The parameters are (a, 6) = (2.49,1.78) for 
the g = 1/2 solution and (2.40,1.73) for the g ^ 0.4 solution. Such a behavior would 
indicate a vanishing imaginary part and thus oscillations only when the Higgs ceases to 
interact. Thus, the oscillations would remain, albeit the poles move closer to the real 
axis. This might only be an indication of the possibility of the Higgs to decay, but the 
question remains into what, as gluons are confined. The existence of one or more complex 
conjugate poles on the same Riemann sheet can be interpreted in the framework of the 
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Gribov-Stingl scenario ^| as a signal of confinement for the Higgs as well. However, as 
the Higgs is sensitive to the truncation, this is not conclusive, until more studies have been 
performed. 



6.2 Thermodynamic Potential 

6.2.1 Soft Mode Contribution in the Infinite- Temperature Limit 

Knowledge of all Green's functions as functions of the temperature would permit to cal- 
culate the thermodynamic potential and therefore all thermodynamic quantities. Since 
not all of these are known, it is not possible to calculate an exact thermodynamic po- 
tential. Knowledge of the propagators is however sufficient to compute an approximation 
to the thermodynamic potential, using the Luttinger-Ward or Cornwall- Jackiw-Tomboulis 
(LW/CJT) effective action 

This effective action can be extended to abelian and non-abelian gauge theories, see e.g. 
|l()()j . Its calculation is prohibitively complicated when using constructed vertices instead 
of bare or exact vertices. Therefore, only some qualitative features using the simplest 
approximation to the thermodynamic potential will be extracted here. It only depends on 
the propagators, and reads 

2^^2-7 (2n)A \Z,{q)) \Z,{q) ) 

\G,{q)) \G,{q) ) 2 Uo(g); 2Uo(g) J 



(6.6) 



where d{G) = 5"" is the dimension of the gauge group. Zq, Gq, and Hq are the tree-level 
dressing functions 

Zoip) = Goip) = 1 
Hoip) 



Working in the infinite-temperature limit, the thermodynamic potential divided by 
becomes interesting. This motivates a rescaling of the integration momenta to obtain 

^ = ^« (6.7) 

where the dimensionless constant a depends only on the dimensionless ratio ruh/g^ and 
thus becomes independent of temperature in the limit T oo. Herein, the explicit value of 
gs enters. Using the t'Hooft-like scahng presented in section IK72| gl ~ Aqcd results. Hence 
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(If) .71) scales as and does not contribute to the thermodynamic pressure significantly 
compared to the hard modes. The expected Stefan-Boltzmann behavior must then be 
obtained from the hard modes alone, requiring 



where ah stems from the hard modes. However, the soft modes may still contribute 
significantly to thermodynamic properties, especially the trace anomaly 9 = e — 3p, near 
the phase transition flUl]. 

The calculation of ()6.6|) turns out to be plagued by spurious divergences, since the 
DSEs with the modified 3-gluon- vertex are no longer exact stationary solutions of ()6.6p . 
By close examination of and the fact that all dressing functions D can be expanded 
in the infinite temperature limit for p'^ g1 a.s 



the divergence structure can be isolated. The tree-level value drops out in each term 
independently. Expanding the logarithms also in g1/q directly shows that the LO cancels 
between the linear and the logarithmic terms as well. This is no longer the case in next- 
to- leading order (NLO), and this will only be possible by requiring a correct solution of 
the Dyson- Schwinger equations to NLO in perturbation theory and using the appropriate 
2-particle irreducible expression [,99^ for the thermodynamic potential. As the NLO term 
is proportional to {gl/qY-, it generates a linear divergence in Vl together with the measure. 
Although there is no a priori reason for this NLO contribution to be the same for both 
solutions, as this term is truncation dependent, still this is the case^. The NNLO will 
generate a logarithmic divergence and only from NNNLO all terms will be UV-finite. The 
general structure of the result is therefore 



where A is the ultraviolet cutoff when calculating (jfi.fij) . 

While the linear term is easily identifiable, this is not the case for the logarithmic one, 
especially as it is always possible to move contributions from c to a and back. This makes 
it in some sense arbitrary to calculate a. Assuming that the divergent contributions will 

^Since the two solutions only differ in their respective infrared behavior, this may have only a minor 
influence on the perturbative tail. 





9lT 



a + h log — \r dk 

c 



(6.8) 
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be small at sufficiently small cutoff, but a has already reached its final value, a first fit is 
performed using 




= a + dA. 



In a second step, with a held fixed, the full fit using ()(i.8|l is performed. This gives the 
final values. For this approach to be valid, the logarithmic term must be much smaller 
than a. At least for calculating the differences of both solutions, this is the case. Also c 
turns out to be always the same value for all fits, giving confidence in the fit procedure. 

In the case of the ghost-loop-only truncation, the difference is independent of the cutoff 
and is —0.001862(0). As all differences were performed by subtracting the g = 1/2 solution 
from the second solution, this result would favor the g = 1/2 solution. 

In all other cases the prescribed subtraction was necessary. Note that, although d 
should be when calculating the differences, this is only true up to numerical errors. 
Hence also for differences it was necessary to extract d and perform the corresponding 
fits^. 48 fit points have been included from the A interval [54763.70,4830021]. Within 
these values, the integrand and the integral did not show significant numerical fluctuations. 
The statistical error on a is obtained by calculating the mean square error of the full fit. 
Besides the statistical error, systematic errors were investigated using variations of the 
3-gluon vertex and the gluon-Higgs vertex by varying 6 and u. The results are listed in 
table ini3 Note that the calculated difference and the fitted difference are always the same 
within statistical errors. 

These results strongly indicate that the g ^ OA solution is the thermodynamically 
preferred one, as this result is stable with respect to variation of the parameters. This 
agrees with the comparison to lattice results in section However, as it still cannot be 
excluded that either of the solutions is a truncation artifact, this is not a final statement. 
Moreover, the value of a seems to be C(l). However, the systematic study of the vertex 
influence indicates that the value can only be determined within a factor of 2 or worse. 
The thermodynamic potential seems to be generated mainly in the chromoelectric sector, 
but also this is only an indication. 

6.2.2 Hard Mode Contribution 

To obtain the contribution from the hard modes, it is at least problematic to use ()6.6j) . 
The LW/CJT-action ()6.6|) considers only the interaction part. It vanishes identically for 
a free system, i.e. for a system containing only tree-level dressing functions. This free 

^The required accuracy for these calculations is very high. An expansion in up to 500 Chebychef 
polynomials, see appendix ^ was necessary. 
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Truncation 


Obiect 


Sol. 


a 


b 


/JO 




Ghost-loop 


An 




-0.001862(0) 








Ghost-loop 


Q 


1 


-0.2095(186) 


-0.258 


189 ■ 10^ 


-0.0025 


Ghost-loop 


Q 


2 


-0.2076(607) 


-0.258 


189 ■ 10^ 


-0.00250 


Yang-Mills 


AQ 




0.03942(0) 


0.00265 


189 ■ 10^ 


1.25 ■ 10"" 


Yang-Mills 


Q 


1 


1.735(135) 


1.92 


189 ■ 10^ 


0.0175 


Yang-Mills 


Q 


2 


1.696(135) 


1.92 


188 • 10^ 


0.0175 


Full 


AQ 




0.05276(0) 


0.00280 


189 • 10^ 


5.54 • 10-12 


Full 


Q 


1 


3.586(254) 


3.67 


189 • 10^ 


0.0322 


Full 


Q 


2 


3.533(254) 


3.67 


189 • 10^ 


0.322 


YM contribution 


AQ 




0.02083(0) 


0.00162 


189 ■ 10^ 


9.92 • 10-12 


YM contribution 


Q 


1 


0.8215(983) 


1.33 


189 ■ 10^ 


0.0108 


YM contribution 


Q 


2 


0.8007(983) 

V / 


1.33 


189 ■ 10^ 


0.0108 


6 = 1 


AQ 




0.08433(0) 

V / 


0.00910 


189 ■ 10^ 


4.90 ■ 10-11 


6 = 1 


Q 


1 


2.442(255) 

V / 


3.55 


189 ■ 10^ 


0.0322 


6 = 1 


Q 


2 


2.338(255) 


3.54 


189 • 10^ 


0.0322 


to = 1/4 


AQ 




0.03820(0) 


0.00258 


189 ■ 10^ 


9.37-10-12 


u = 1/4 


Q 


1 


1.708(136) 


1.93 


189 ■ 10^ 


0.0171 


u = 1/4 


Q 


2 


1.669(136) 


1.93 


189 ■ 10^ 


0.0171 


6=1, to =1/4: 


AQ 




0.08282(0) 


0.00999 


189 ■ 10^ 


6.9610-11 


5 = l,tu = 1/4 


Q 


1 


0.4518(4210) 


1.79 


189 ■ 10^ 


0.0171 


6 = 1,00 = 1/4 


Q 


2 


0.3689(1362) 


1.78 


189 ■ 10^ 


0.0171 



Table 6.4: The thermodynamic potential Vt and the difference between both solutions 
AVt in the different truncation schemes. The fit was performed using ()6.8|) . The quoted 
errors are statistical only. Solution 1 is (7 = 1/2 and solution 2 is (? ~ 0.4. Note that 
ln(4830021/189000) ^ 3.25. YM contribution is the Yang-Mills part of the full solution. 
The 6-, 10- , and ^cj- variations stem from the variation of the 3-gluon vertex, the gluon-Higgs 
vertex and the combined variations. 
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contribution has thus to be added exphcitly. Therefore the hard mode contribution can 
essentially be calculated to be the same as that of a non-interacting system of free gluons. 
The LW/CJT-expression fl6.6p does not capture the tree- level contribution of the soft 
modes, and it thus can be added here. This yields a thermodynamic potential of a gas of 
massless gluons with small corrections due to the explicit soft-mode contributions and the 
residual interactions of the hard modes. It is for Ai'c = 3 jH] 

n 167r2 ^ ^ 

lim — ^ . 6.9) 

Using the t'Hooft-scaling for the interaction strength, the corrections due to the soft contri- 
butions vanish like 1/T. The only remaining contribution in the infinite temperature limit 
is the Stefan-Boltzmann-limit ()fj.9j) . in agreement with results from lattice calculations [Hj. 
Therefore the thermodynamical properties of the high-temperature limit are governed by 
the hard modes, although they do not participate in the dynamical interactions. 



Chapter 7 

Conclusions, Summary, and Outlook 



7.1 Concluding Remarks 

Before discussing the implications of the results found in chapters |5 and IHl on the 
structure of the Yang-Mills (and QCD) phase diagram, it is worthwhile to take a step 
back to gather and assess the findings. Especially the reliability of the results deserves 
special attention. 

The primary objects investigated in chapters lU and El are the propagators of Yang- 
Mills theory in the high temperature phase. These can be grouped into three classes. 
One consists of the soft modes of the ghost and the 3d-transverse gluon, constituting 
the Yang-Mills sector of the limiting 3d-theory. The second is the soft 3d-longitudinal 
gluon, corresponding to the Higgs sector of the 3d-theory. From the point of view of 
the 4d-theory, these are the soft chromomagnetic and chromoelectric degrees of freedom, 
respectively. The last class contains all the hard modes. 

The first class on its own constitutes a 3d Yang-Mills theory. Its properties are found 
to be very similar to the corresponding 4d-theory. An infrared divergent ghost propagator 
is accompanied by an infrared vanishing gluon propagator. Therefore the gluon is confined 
and the results fulfill the Kugo-Ojima and Zwanziger-Gribov scenarios. Also the Schwinger 
function and analytic structure of the gluon propagator are similar to the 4d case. 

The systematic study of the errors due to the truncation find a very weak quantitative 
effect only. Especially the confining properties are extremely robust against variations in 
the truncation. Finally the comparison with lattice calculations on exceptionally large 
lattices shows a very good agreement. Thus it is very likely that these results on a 3d 
Yang-Mills theory are reliable. 

The situation is more subtle concerning the Higgs in the 3d-theory. Its propagator 
shows the typical behavior of a massive particle and it is dominated by its mass. However, 
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the comparison to lattice data indicates that to some extent subleading or probably non- 
perturbative effects are relevant. This is much more pronounced when regarding the 
Schwinger function. Here significant changes are found when modifying the Higgs-gluon 
vertex. Thus the qualitative behavior of the Higgs propagator seems to be under control 
and the presence of non-trivial effects is quite certain. The consequences of these effects 
on the other side are far less conclusive and require further investigation. Especially the 
question whether the Higgs is confined is of great interest. On the other hand, a quite 
firm statement is that it has only little effect on the Yang-Mills sector. 

Combining both classes, a qualitatively stable infinite temperature limit is found. Only 
the chromoelectric sector exhibits screening while at least in the chromomagnetic sector 
confinement prevails. Therefore the infinite temperature limit is non-trivial and the high 
temperature phase strongly interacting. Hence the simple picture described in the intro- 
duction in section li.jjl is not applicable. This and the detailed analysis of the infinite 
temperature limit are the central results of this work. 

Less reliable is the exploratory study of finite temperature effects in chapter E] and 
thus the properties of the hard modes. The success of the truncation scheme in the 
3d-theory is based on its exactness in the ultraviolet and its presumed exactness in the 
infrared. Thus the consequences of the deficiencies at intermediate momenta are strongly 
constrained. This possibility is lost by cutting off the Matsubara sum. As a consequence 
the ultraviolet properties of the hard modes are incorrect and a rooting in perturbation 
theory is prevented. This amplifies the second problem. The hard modes do not reach into 
the 4d-infrared due to their effective mass. Therefore the advantages of the truncation 
scheme are lost to a large extent and the equations for the hard modes become unreliable. 

Nonetheless, the results found are not irrelevant. At sufficiently large temperatures the 
truncation-dependent self-energies of the hard modes are suppressed by the large effective 
mass and they are dominated by the truncation-independent tree-level terms. Therefore 
the truncation artifacts vanish when the temperature goes to infinity, reproducing the 
reliable 3d-theory. The hard part is then to assess down to which temperature the results 
are still at least qualitatively reliable. 

As it is found that the system is only very weakly dependent on temperature, qualitative 
conclusions can probably be drawn even at quite small temperatures. This is due to the 
fact that the infrared properties of the soft modes are nearly independent of the hard 
modes. In general the hard modes effectively decouple, because of their effective mass, 
which is large compared to Aqcd even at the phase transition. The only exception is the 
generation of the screening mass for the soft 3d-longitudinal gluon, which is found on quite 
general grounds and thus can be considered also a reliable conclusion. It is also supported 
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by the systematic error estimations performed, which do not find quahtative effects but 
again only small quantitative ones. 

Thus, although the quantitative results may be subject to change, the qualitative 
result that, probably down to temperatures within the same order of magnitude as the phase 
transition, the high temperature phase consists of strongly interacting soft modes exhibiting 
confining properties and inert hard modes is expected to be quite reliable. This is also 
supported by lattice results, which find that the infinite-temperature limit is effectively 
reached at about 2Tc |H21 and the qualitative features of the infinite temperature limit 
up to the highest temperatures, for which propagators are available, of 6Tc This 
constitutes the second major result found in this work. 

Concerning as the last point the thermodynamic potential, no final conclusion can be 
drawn. The results indicate that the gross thermodynamic properties far away from the 
phase transition are completely dominated by the hard modes and that a Stefan-Boltzmann 
behavior is reached in the infinite temperature limit. Considering the crude assumptions 
made and the difficulties encountered, this result is indicative at best. Nonetheless, if it can 
be substantiated, it would allow to understand how a Stefan-Boltzmann-like behavior can 
emerge from a non-trivial theory. At the same time, the results also indicate that at low 
temperatures non-perturbative effects likely will play a role in the thermodynamics and 
will especially be relevant in the vicinity of the phase transition and thus to experiment. 

7.2 Implications for the Phase Structure of QCD 

Gathering the results found in chaptersEJEl andEland the solutions at zero and small tem- 
peratures presented in sections ESI and EHl a coherent picture emerges. Still, this picture 
is overshadowed by the problems of truncation artifacts, especially at finite temperatures. 

Combining these results has impact on the understanding of the phase transition. 
The main difference between the low-temperature and the high-temperature phase is not 
primarily one between a strongly interacting and confining system and one with only 
quasi-free particles. The chromoelectric gluons, whose infrared behavior change from 
over-screening to screening, come somewhat close to such a picture, and the hard modes 
certainly do. The latter become free due to the vacuum property of asymptotic freedom 
and the dependence of their energy on the temperature, thus in an expected although not 
entirely trivial manner. The chromomagnetic gluons remain over-screened in the infrared 
and are thus confined. At least in the transverse sector, the results provide strong evidence 
for the Zwanziger-Gribov and/or the Kugo-Ojima scenarios. 

The order parameter for the phase transition is then necessarily only a chromoelectric 
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one. This is in agreement with studies of Wilson loops in lattice calculations. There it 
is known that only the temporal (chromoelectric) Wilson lines show a behavior typical 
of deconfinement, while the spatial (chromomagnetic) ones do not |lU2j . Note that the 
order parameters used to study the deconfinement transition on the lattice are typically 
chromoelectric ones like the Polyakov lines [9^. This leads to the conjecture that if the 
corresponding chromomagnetic Polyakov lines were to be used almost no change would be 
found. 

Including the observation of super-heating in the low temperature regime [OBI and 
super-cooling in the high temperature regime, the scenario emerging is a chromoelectric 
phase transition of first order. This phase transition connects two different strongly in- 
teracting phases of Yang-Mills theory. Especially quantum fluctuations always dominate 
thermal fluctuations and at least part of the gluon spectrum is always confined. 

In the vicinity of the phase transition, the non-perturbative effects found are likely also 
relevant to thermodynamic properties, especially to the pressure and the trace anomaly 
This underlines the importance of non-perturbative effects at least for the temper- 
ature range relevant to experiment. 

It is still unknown which effects can be induced by quarks. In the infinite temperature 
limit, arguments exist that they will not contribute due to their fermionic nature. This 
nature requires anti-periodic boundary conditions and therefore generates effective 
'masses' of {2n + 1)ttT with n integer. Thus all quarks become infinitely massive in the 
infinite temperature limit, irrespective of their intrinsic mass and thus decouple due to 
asymptotic freedom. Lattice calculations indicate that the quarks possibly are able to 
change the nature of the phase transition up to a cross-over. It is hard to imagine how 
such a partially confined phase would behave, and thus the topic is of high interest. 

In addition, the observation of a drastic change in the inter-quark potential [22] and its 
connection to the residual confinement of gluons still needs to be understood. However, 
it cannot yet be firmly concluded that the inter-quark potential does not rise only loga- 
rithmically at high temperatures, as it would be natural to expect from a dimensionally 
reduced theory. On the other hand, quark confinement is still not understood at zero 
temperature, thus posing a more significant task. 

7.3 Summing Up and Looking Ahead 

Summarizing, in this work evidence is presented in agreement with lattice calculations that 
the finite temperature phase transition of Yang-Mills theory does not lead into a trivial 
phase. By investigating first the dimensionally reduced theory, thereby providing results 
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on genuine 3d-theories of Yang-Mills and Yang-Mills-Higgs type, the non-triviality of soft 
interactions was found. Extending the range to finite temperatures, this was confirmed and 
the decoupling of the hard modes explicitly demonstrated. This establishes the main fea- 
tures of the high temperature phase of Yang-Mills theories in Landau gauge. Furthermore, 
using the Schwinger functions, the absence of at least part of the gluons from the physical 
spectrum was found, making manifest the existence of residual confinement. The results 
found comply with Zwanziger-Gribov or Kugo-Ojima type confinement mechanisms. As 
a final point, a first step towards thermodynamic observables was performed. 

Taking a look ahead, two issues are of high interest. One from an experimental and 
one from a conceptual point of view. 

In heavy-ion collisions, not only finite temperature is of interest, but especially for 
low-energy collisions results at finite density are of importance. Such experiments are 
performed at the SPS at CERN and have been performed at the AGS at Brookhaven 
and will also be conducted at the planned FAIR facility at GSI. Finite density is also 
relevant for the physics of compact stellar objects as well. For sufficiently large chemical 
potentials (densities), at least the quark sector will likely be weakly interacting, although 
non-perturbatively, see e.g. jlU4j . However, it has not yet been convincingly demonstrated 
how the quarks can generate gluon deconfinement. In addition, those chemical potentials 
probably only occur in nature during black hole formation. For chemical potentials relevant 
to experiments and possibly neutron stars, strong non-perturbative effects are expected. 
Thus a framework as the one used previously for full QCD [5T] seems to be necessary. 
First steps in this direction have been performed |in5j and are pursued further, they, 
however, are not yet as advanced as the finite temperature calculations presented here. 
More results are available from studies using model gluon propagators [7^] and semi- 
perturbative methods |lU4j . 

On the conceptual sides the problems encountered in chapter El to obtain the vacuum 
solution indicate that it will probably be necessary to choose other means to investigate the 
phase transition directly, e.g. real-time methods jSH]- In addition, the inclusion of quarks 
is an urgent problem as their impact on thermodynamic properties is highly relevant for 
experiments, including RHIC at Brookhaven and ALICE, which is under construction 
at the LHC at CERN. In addition, understanding of quark confinement in the vacuum 
might be supported by understanding the high temperature behavior of quarks and the 
phase transition. This would also lead to a deeper understanding of the relationship 
of chiral symmetry and confinement of fundamental quarks, for which evidence exists. 
Indeed, their phase transitions seem to be interlinked jH]- This immediately leads to the 
problem of the connection to topological excitations, which lattice calculations indicate 
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to play a dominant role in the physics of confinement and chiral symmetry breaking. If 
they are relevant, these excitations must also be present in the current approach, as their 
consequences are. Therefore, it would be desirable to disentangle their contributions and 
identify their role. It is at this point also completely unclear, how their dynamics change 
above the phase transition and their connection to the at least partial confinement of 
gluons. 

To reach these aims, implementation of new methods and close cooperation with other 
approaches, especially lattice gauge theory, will be necessary. Nevertheless, over the last 
decade significant progress has been made in the understanding of low energy QCD, pro- 
viding hope that it will eventually be understood. 
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Appendix A 
Conventions 



Calculations are done in natural units, 

h = c = ks = I- 

The Minkowski-metric is 

77 = diag(l, -1,-1,-1). 
The Euclidean metric after factorization of -1 is 

r]e = diag(l, 1, 1, 1). 

Thus p"^ = Po+ P^- Sometimes \p\ = pa is used. The individual components of p do never 
appear exphcitly. 

The second Casimir Ca of a compact, semi-simple Lie-group G in the adjoint repre- 
sentation with structure constants /"'"^ is defined by 

j^acd j^acd Cj^ 

For SU(iV,) groups, Ca = N^. 

Fourier transformations are performed with all momenta incoming. 
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Appendix B 

Derivation of the Dyson- Schwinger 
Equations 

Starting from (jH.lj) . it is possible to derive the DSEs for the Lagrangian ()2.H|1 . In this 
context also the Legendre transform of W plays a role, where W is defined by the generating 
functional ()2.10|) . The Legendre transform can be obtained by 

W{f) = -r(0'^) + j d'^xf{x)(l)''{x) (B.l) 

which entails 



f - ^. 

In the case of Grassmann fields u like ghosts and fermions, two independent sources are 
necessary. This modifies the above to 

6W 6W 



Su^-^x) Su^-^x) 

where all derivatives with respect to Grassmann variables act in the direction of ordinary 
derivatives. 
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The general procedure to obtain the corresponding DSEs is to calculate equation 
and then take the derivative once more with respect to the field or with respect to the 
anti-field in case of anti-commuting fields. The additional source terms then yield the full 
propagators, while the right-hand-sides of the equations are found by the derivative of the 
action. 

Since in the course of the derivation, the source in equation (jH.lj) becomes the inverse 
full propagator, it makes sense to already rewrite p.ip as 

= , , , Z (B.2) 



at the sources set equal to in the end. This entails also that all single derivatives of the 
effective action, being classical fields, vanish at the end. 

The inverse propagators of the ghost, the gluon and the Higgs are defined as 



6^T 



5d>{y)6c''{x) 
_JH__ 

5^T 



The propagators are then given by 

5^W 



■■D^\x~y)-' (B.3) 
D'Hlix-y)-' (B.4) 
D'}^{x-y)-\ (B.5) 



5rf'{y)6ri"-{x) 
5^W 



Sk^{y)6k<^{x) 



: D^\x - y) (B.6) 

DfA^-y) (B.7) 

Dii{x-y). (B.8) 



7] is the source of the ghost field, t] of the anti-ghost field, j of the gluon, and k of the Higgs. 
flB.6IIB.8j) are inverse to ()B.3IIB.5|) which can be proven, e.g. for the ghost propagator by 



^ ^ 5r]''{z)Sf]''{x) 5d'{y)Sc''{z) J ^ ^ 5r]''{z) 5c^{y) 6&{y) ^ 
Note that all mixed propagators vanish when the sources are set to 0, e.g. 

5^V 



6c9{w)6AI{z^ 
due to conservation of ghost number and spin 



j=f]=ri=0 



0, 
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Higher n-point functions, the vertices, are defined accordingly as 

and correspondingly for the other vertices. 

There are a few useful identities in the course of the derivation. These are 

d zd w- 



6j-{x)6r]'i{x) J 6jl{z)Sj;{x) S&{w)5AI{z) Sr]9{w)6r]''{x) 

d^.J- = fd'z(df.- ^ ^ 



'^dri^i^x) J \ ^6t]''{x)6c''{z) J 6r]''{z)' 
Further, since 



d zd w 



SAliy)6f^{x)6f^{x) 6Al{y) J 5j^(a;)5j/(^) 6Aiiz)6Aliw) Sj^{w)6j-^{x) 

+ / d'^zd'^w- 



5Al{y)53-„{x)5l{x) J 6f^ix)6jliz) 6Aliy)54iz)6AUw) Sf4w)6f^{x) 

it follows that 

= — / a zd w 

6Ai{y)6j-^{x)6j-^{x) J Sj-,ix)5ji:iz) SAiiy)SA{,iz)SAUw) 5j^{w)5jl{x) ' 

(B.9) 

when the sources are set to 0. If not all sources and fields are of the same type, then it 
is necessary to sum over all possible intermediate fields uj"" and their sources where the 
index a includes also the field type. Thus 

- / d'^zd'^w ^^ , / X / xr, / X - (B.IO) 



6u'>{y)6l''{x)6l''{x) J 6l''{x)6lf{z)6uj^{y)6ujf{z)6uj9{w)6ls{w)6l''{x 
It is also important to note that 

5^W 



is not zero, despite the translational invariance of the two-point functions. The derivative 
of the propagator at does not necessarily vanish especially, if the propagator does only 
exist in the sense of a distribution. 

As the effective potential defined by (jB.lll is bosonic, it can only depend on even 
numbers of ghost and anti-ghost fields. Therefore all vertices with only one ghost or 
anti-ghost leg vanish, when the sources are set to zero. 

After a tedious calculation, the DSEs in position space are obtained. Adding suitable 
combinations of permuted diagrams, it is in all cases possible to make the appearance 
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of one tree- level vertex explicit. Performing a Fourier-transformation, and keeping the 
redundant momentum at the vertices explicit, the DSEs in momentum space are for the 
ghost ^ 

D'§{p)-^ = -S^Y 

+ I (0rr«'--^''^-(-g,p, q - p)D%{p - q)D''j{q)rf^'''9f{-p, q,p-q), (B.ll) 

the gluon 

- / (|^rJ^'^''^-(-p - g, q,p)D'i^{q)Di^{p + q)Tf^J-\-q,p + q, -p) 

+\ I (&r^i;f '-'^(p, q - p, -q)D^MDiUp - g)r£f ^(-p, Q,p-q) 



+ 1 1 Tlir"'^'^'''^'='(P> Q - P, -q)K'iq)D^H0f)ip - q)T^^"'''^i-p, q,p-q) 



{2-kY 

pi 

+ \ I (0r*l.f/^^'^(p, -p, g, -q)D^/^{q) 

■Df.ip - q - k)D%{q)DZ{k)Til'^^%p - q - k, -p, q, k) 
+i / '^^Xt'^^'^^P^ -q,q-k-p, k)D^\{k + p - g)!^^^/^-^, P + K -k) ■ 

■r^^f '(fc + p - g, g, -fc - p)Dll{p + A;)D^^JA:)D,- (g) 
+ \ I '^r^r'^p, -q,q-k-p, k)D^\{k + p - q)TfJ''^f{~p, -k, p + k) ■ 

+ p - g, g, -A: - p)Dgip + k)DZ{k)Dt^{q) 
+ \ I (0r*i;.^''^''"''^(P> -p. 1, -<l)D%{q) 
-\ I l^^ff (P^ q-p-k,-q, k) 
■Dl%p - q + k)Dj^{q)Df^{k)Vff^3'^\p - q + k, -p, g, -k) 
+ I ^^if'^^'^^'^iP, Q-p-k,-q, k)Dl%p + k - q)Df^{k)D%\k + p) ■ 

■D|(g)r^<^''fj'=(p + k-q,q,-k- p)T^^'''^'''{-p, k + p, -k) 
+ I ^r%^''^''^^''f{p, q-p-k,-q, k)Dl^^{p + k - q)D{iik)Dtiik + p) ■ 

■D%{q)V'^l^^'^={p + k~q,-k~p, g)r£^'^'^^(-p, k + p, -k) 
+ / 'i^^^lf^'^^^'^P. k,q-p-k, -q)D^^{p + k - q)Dl;{k)D'},'{k + p) ■ 

■Df{q)T^''3^^{p + k - q,q,-k - p)Tfp^'^'''{-p, -k, k + p) 
+ / 'i0^ff^"'^^''{P. k,q-p-k, -g)Df (p + k - q)Dgik)D^,tik + p) ■ 

■D%{q)Vt^'^^3^{-k -p,p + k-q, q)Tf^t''i-P^ -k. k+p), (B.12) 

and the Higgs 

D'§{p)-^ = 5''\p^ + ml) + 
^Note that some diagrams which vanish for tree-level vertices had not been included in |81| . 
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+ \ I ^^ff^'^-'^'i^^ -q^P^ -p)D%{q) 

-\ I 'i0^tf'^''''"'%-P -q + k, -k,p, q)DlUP + q - k)D%{k) ■ 
■Dfl{q)Vfj''''9^''{-k,p + q-k, -p, q) 
+\ I ^^if^''''''%-q, k,p,-p + q- k)D%{p - q + k)D;i{q)Dli{k) ■ 

.Y^<t>''9i''{q,p- q + k,-p- k)Df{j> + k)V^^^^^^^{-k, -p,p + k) 
+\ I ^^X'^''''''%-q, k,p,-p + q- k)D^^{p - q + k)D^^l{q)D;i{k) ■ 

■Ttl'^''''{q, -p-k,p-q + k)Dt{p + k)V^^^^^^\-k,p + k, -p) 
+\ I tS7^r^'x^'^'''''^1-?, q + k,p,-p- k)D'll{q)Dii{p + A;)Dg(g + k) ■ 

■^if\k, q, -q - k)D%{k)Tf'^^\-k, -p,p + k) 
+i / tST^r^'x^'^'''"^^-?. q + k,p,-p- k)D^^MDi}{p + k)Di^{q + A:) ■ 
■r^?'''''(g, -q - k, k)D%{k)V'*''''^^{-k, -p,p + k) 
+ \ I (0r*i'*^'^^^'^(p, -p, q, -q)D^^{q) - 

-\ I ^^''^'^"'"''Hp, -p + q-k,k, -q)D'^{p - q + ■ 
•L)^^(A;)r<^*''^^^-^(p -q + k,q, -p, ~k) 
+\ I gpr*i"^'''^^^^(p, -p-k + q,k, -q)D^^{p + k - q)D^^{q)Dfj^{k) ■ 

.Y4>^9dh^q^ p + A; - g, -A; - p)Df}{p + fc)r<^''^^^(p + k, -p, -k) 
+1 / tST^r*i"^'''^^^^(P, -p-k + q,k, -q)D^{p + k - q)D^^{q)Df^{k) ■ 

.YM\hgd^_l^ - p, g,p + A; - q)Df^{j> + k)T^'''''''^f {p + k, -p, -k). (B.13) 

Several of the appearing vertices do not have a tree- level counterpart. The tree-level 
vertices for the ghost-gluon-, 3-gluon-, 4-gluon-, 2-gluon-Higgs-, 2-gluon-2-Higgs- and 4- 
Higgs- interactions have been used. They are given by 

r^hc-cA,abc^p^^^f^^ = tgar'^ (B.14) 

^'ilf'^'^P^q^k) = -igdr"'{{q-k),5,p + {k-pU,, + {p-q),5,,) (B.15) 



^%fp''''\p,q.k^) 




^fjLp^ua) 






(B.16) 


Yt\M\ahc^^^ g, k) 


= ^g^r'^q-k). 


(B.17) 


v'l:f^"'''''\p.q.Ki) 




(B.18) 




= 2h{5ab5cd + 5ac5bd + SadSbc), 


(B.19) 



where the momentum-conserving 5-functions have been suppressed. The complete graph- 
ical representation of ()B.lHIB.13|l is shown in figure 13.11 
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C.l 3d Kernels 

The integral kernels in ()4. 9114 .11]) are obtained by using tree-level vertices and performing 
the corresponding contractions. The ghost kernel is given by 



(A;2 + g2_2A;gcOS^)2' 

The contributions in the Higgs equation are 

2g2 sin3(^) 



N2{k,q,9) 



{k^ + q^ + 2kqcos9y 
2sm^{9) 



k'^ + q^ + 2kq cos 9 
The kernels in the gluon equation are finally 

jDft. r, n\ - ((C-i)fcg cos(e)-gZ+Cg^ cos2(e)) sin e 

^y'^^'d^^) 2fc2(fc2+g2+2A:gcos6l) 
U (t. r, ft\ - ( (C- 1) (fc^ +4fcg cos e) -4g2 +Aq^ C cos^ (0)) sin 6 

^^^Lyi^,q,^) — 4fc2(fc2+g2+2fcgcose) 
Mxik, q,9) = 4fc2(fc2+^2+^2fegcos9) ■ 
■ + 2g2) ((( - 9)P - 4g2) + 8(C - 3){k^ + q^)kq cos 9 + (8Cg^ + (C + 7)A;^ 

+4(5C - l)A;2g2) cos2(e) + 4(4Cg^ + (C + 3)A;2) cos3(^) + ACk^q^ cos\9)^ . 
The modified gluon vertex is introduced by multiplying Mt with ()4.H(i|l . 
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C.2. Finite Temperature Kernels 



C.2 Finite Temperature Kernels 



At finite temperature, the kernels are obtained from tree-level vertices. The following 
abbreviations have been used: 

k ^ k 

2 I 2 

X = q +qo 

y = k^ + kl 

us = {q-kY 
Z3 = {q + kf 
U4 = (/co - go)^ + U3 
Z4 = {ko + qof + zz- 
The kernels in the ghost equations are 

k'^q^shi{ef 



Ariko, qo: k, q, 9) 



AL{ko,qo,k,q,9) 



xyusUi 

q^ sin 9{q^ko + k^qo — kq cos 9 {ko + qo)y 



xyusu^-" 

The kernels in the transverse gluon equations are 

q'{-q + k{-l + Ocos9 + qCcos{9f)sLn{9) 

n[Ko, qo, K, q, U) 
Mxiko, qo, k, q, 9) 



2xyz4 

f{-A + C)+k^{-29 + 5C) 



Mi{ko, qo, k, q, 9) = 



+2k'^q^ (-46 + 15C) ) sm(9) + 4A;g(|V (-3 + 2C) 

+e (-9 + 50 ) sin(2^^) + (Sq^ + k\7 + Q 

+k'^q'^ (-4 + 23C) ) sin(3^) + 2kq (4g^C + (3 + C)) sin(4^) 

+k^q^CsH59)^ 

^ll^ ( - (-4 + C) ^o' - k' (-9 + C) qo' + kV i^ko - qo) ■ 



■ {Ako + (-8 + C) go) - 2kqcos9(^ - IQq^ko^ 
+2 (8A;2 - g2 (_8 + Q) koqo + k^ (-13 + C) go') + cos(2^) • 
• (4g^C^o' - 4A;2g^ ((7 + C) A;o - 2go) go + A;^ (7 + C) go') + 
A;ggo (A;gCcos(4e)go + 2cos(3^) {-2q\ko + k^ {3 + C) go)) ) 
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2 

M2{ko,qo,k,q,e) = —A__(^-(^kVCsm{5e){ko-qof)+^kqsm{29)- 

■ (ko - go) {2q' (-4 + C)ko + k' ((-1 + C) - (-7 + C) go)) 
-2A;gsin(4^) {ko - go) {2q\ko + k' ((-1 + () A;o - (1 + C) go)) 
+ sin e (Aq"" (-8 + 3C) ^o' + k^(^- Sko'' + 3C(^o - go)' + 6A;ogo 

-23go^) + 2kY ((-10 + 70 k^' + (34 - 8C) ^ogo + (-4 + C) go') ) 

- sin(3^) (4gt A;o' + k^q^ ( (4 + 3C) ko^ -2(10 + () A;ogo 

- (-8 + C) go') + k^ ( (-1 + C) ko' - 2 (-1 + C) A;ogo 

+ (3 + go'' 



2 • n 

ML{ko,qo,k,q,9) = (fcV (-1 + C) - 4g'(x + fco' + Mo)' + 4A;VCcos^^ • 

•(2g2 + qo {-ko + go))' + 2k^ (2q^ (-1 + () ko' + 3q^ (-1 + () A;ogo 
+x (g2 (-3 + - 2go^) ) + kY (^' (-9 + C) + 4 (-1 + C) ko' 
+2x (-7 + 30 A;ogo + 12 (-1 + A;o'go + ko' {^q^ (-3 + C) 
+ (-21 + 130 9o')) + ^kqcos0^ (2g' + go {-ko + go)) • 

• (2g2C {x + ko' + koqo) + k^ {f (-2 + 40 + go (^o - Qko + go 
+C?o))) + 2kqcose(2q^ {x + ko' + koqo) {q' (-5 + 

+2 (-1 + C) ko' + (-1 + 30 Mo + (-3 + go') + k' (-1 + • 

• (4g2 + qo {-ko + go)) + A;^^ - 2 (-1 + A;o'go + (-1 + ko' ■ 
■ {W - go') + X (2g2 (-7 + 30 + (-5 + C) go') + A;ogo (g' • 

• (-9 + 130 + 2 (1 + go'))) + cos9^(4q\{x + ko' + A;ogo)' 
+AkV (g^ (-8 + 60-3 (-1 + C) A;o'go + Sq' (-1 + C) go' + 
(-1 + 20 go' + ko' {q' (-6 + 80 - (-2 + go') + 2A;ogo(g' • 

• (-1 + 40 + (1 + Qo') ) + k' (4g^ (-5 + 60 - M'qo (3 (-1 + 
■ko + go - 3Cgo) + go' ( (-1 + ko' - 2 (-1 + C) A;ogo 

+ (3 + 0?o'' 
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The kernels in the longitudinal gluon equations are 

sin Q 



P{ko, qo, k, q, 0) = ^ f - 2 - 1) klqo {ko + go) + k'^Qo {ko - ^ko + ^qo) + qko cos 6 

X U Za ^ 



(q^ cos eko + k(ko + 2qo -^{ko + 4go))) 



Nxiko, qo: k, q, 9) = ^ ^^^'^^ (SA;^ + + 8A;g cos O + k'^ cos{2d)) (ko {ko + 2qo) (k 



+ko {ko + 2qo) ) - e (A;o^ (-g + ko + 2go) (g + /cq + 2go) + (ko^ 

+2koqo - 2go^) ) + qko (q^ cos{29)ko + 2k cos 9 (ko + 2qo 

-e {ko + 4qo) 
q^ sin^ ^ 

^kW (go (A:o + 2go) (2g' + go (A;o + 2go)) + i {^q^ + g'go ( - 2A;o 

+go^ - go^(^o + 2go)^^^ + kqkoqo(kq^cos{29)koqo -2cos9 

■ {2q'^ko' + k^ [q^ {Ai - 2) + go {{i -l)ko- 2qo + 4ego) 
^ sin^ I 



N,{ko,qo,k,q,9) = (2gW + (2^'^ - 2^ - 1) Mo - - 1) V) 



N2{ko,qo,k,q,9) = (/c'/co (feo + go) [^k^ + ko {ko + 2qo) (2q^ + {ko + go) (fe 

+2go)) + k^ (6g2 + (3A;o + go) {ko + 2go)) ) + c(2A;' + 2g^A;o' 
+2A;^ (4g2 + A;o' + A;ogo + 2go') - kV ( - 4g^ + ko'' + 6A;o^go 
-g'go' + 4go' + 4A:ogo (g' + 3go') + ko^ (-Sg^ + 13go') ) + k" 
■ {2x' - ko" - Qko\o + ko^ (6g2 - 3go') + ko {-6q\o + 2go')) ) 
+kq(kq cos{29) {Ak^^ + ^ko^{ko - qof + Ak'^ko {ko + go - 2ego)) 
+2 cos 9 (k'^ko {ko + go) {Ak"^ + 2g2 + (3A;o + go) {ko + 2go)) 
+e (Ak^ + 2g2A;o^ (A;o - go) + 2k^ {2x + A^o' - A;ogo) - A;'A;o (A;o^ 
+4a;go + lO/co^go + ko (-4g^ + 9go 



NL{ko,qo,k,q,9) 



—q^ sin ^ 



2A;2g2 cos^ 9 (ik^^qo'' + C^o' (2g' + go {-ko + go))' + ^k^koqo (g' (2 - 4^ + go (a^o + go 
-2ego))) + 2A;gcos(^)' {^k^^qo^ + 2q^^ko' i^q"" + 2g2A;o' + g'A;ogo - A;o'go + Sg^go' 
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+qo') + '2k% ((2 + ko'qo + ^^Qo {2q' + Qo') + h {q^ (6 - 12^ - (-2 + go')) 
+fc2fco( - ((e - 3) ko'qo') - ^ko'qo {q' (-3 + 9^ + 5 - 1) go') - 2 (-1 + 2^ go(6g^ 
+6g2go^ + go') + h (sq'C - 22q' (-1 + go^ - 9 - 1) go') ) ) + g' (3^ - 1) 
+k% ( (-1 + 50 ko^ - 2x (-1 + 30 go + 2 (-1 + 3^ A;o'go + 2 (1 - 3^ A;ogo') - 
(e - 1) A;o'(2xgo + A;o'go + ko (g' + 3go2))' + k'(2^ko'' + ^ko'qo + 2x'^qo' + (3 - 5^) 
•fco'go' - 2A;o=^go (g' (4^ - 2) + (9^ - 5) go^) + ^o' ( - (g' - 1)) - 2g^ (7^ - 5) go^ 
+ (9 - 110 qo") + 2koqo (- (g' (-1 + 0) + 2g'go' + (1 + go') )) - 2kqcofie 
■ (k' (3e - 1) • A;ogo + A;' ( (-1 + 5^ A;o'go - 2x^qo'' + koqo (f (7^ - 3) + 2 - 1) go') 
+fco' (g' (2 - 60 + 4 (-1 + 20 go') ) - ^0 ( - 2q^x^ (-1 + 20 go + 2fco'go (g' (1 + 
-3 (-1 + go') + ^o' (2g'e - (-1 + ?o') + ^o' (g' (1 + + (13 - 110 ?o' 
-13 (e - 1) go') - xko (g' - 1) + g' (17^ - 13) go' + 4 - 1) go') - 2A;o'go (4g' 

• (e - 1) + g' (14^ - 13) go^ + 6 - 1) go') ) + A;^ (2eA;o'go + A;o'go (g' - 3) + 2 

• m - 5) go') + ^o' (g' (1 - 70 + (-3 + 70 go') - 2ego' (2g' + ^qW + go') - ^o' 

• (g^ (1 + - 4g2 (-4 + 50 go' + (9 - 70 go') - 2A;ogo (g' (3 - 50 + g' (4 - 30 go' 

+ (1 + ?o') ) ) + cos^ e (2k^^qo^ + 2g^e^o' (x + ko'' + koqo) ' + 2k'' qo ( (1 + ^o'go 

+2ego (2g2 + go^) + ko (g^ (6 - 140 + (1 + go') ) + ( - ((-3 + A:o'go') 

-2A;o'go (g' (17^ - 5) + (3^ - 5) go') + 2ego' (8g' + Sg'go' + go') + A;o' (4g^ (5^ - 1) 

+2g2 (19 - 230 ?o' - 3 (-3 + go') + 2A;ogo (g' (12 - 220 + 2g' • (5 - 40 go' 

+ (1 + ?o')) + k'ko[- ((-1 + ^o'go') + 2A;o'go(g' (1 - 50 - 3 (-1 + go') 
-8g2 (2e - 1) go (2g' + 3g2go' + go') + ko' (l6gt + 2g2 (11 - 150 ?o' - 13 (C - 1) go') 
-4A;o'go (g' ( - 5 + 4^) + q' (-14 + 150 go' + 3 (-1 + go') + 4fco (g^ (1 + - 2g' 

• (-7 + 80 go' + g' (11 - 130 5o' - (-1 + 50^)))) • 

The vertex dressing in the soft-soft Mr-kernel is performed as in the 3d-case. 
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Infrared Expressions 



D.l Massive Higgs 

As already argued in section \4.2.2\ the only solution without further assumption is that of 
ghost dominance. This then requires the calculation of Igt and Igg in ()4.18|) and ()4.19|) 
only. The latter can be obtained straightforwardly when using the ansatz ()4.16|) and the 
general formula 

{2nY^ (4vr)f T{d + a + P)T{-a)T{-P) ^ ' ^ ^ ^ 

valid for finite integrals. It yields 

qICat^ '2'^-'\d - Ag){2 + d{C - 2) - 4^?(C - 1) - QT^d - 2g)V{2g - f ) 



(47r)f {d-l)g^T{l±^yT{gY 
Note that the integral is convergent if and only if 

d-l d-2 , 

^>9> (D.2) 

where the equality on the upper boundary requires the result to exist only in the sense of 
a distribution. 

Performing the same calculation for the ghost self-energy is more complicated. Since it 
is, in general, a divergent quantity, equation (jP.ljl cannot be applied directly. It is neces- 
sary to regularize and then renormalize the expression. This can be done in a momentum 
subtraction scheme 03] or via dimensional regularization. Here, the latter was performed 
by applying the standard rules of dimensional regularization [2H]- This immediately gives 
a finite result in odd dimensions, as dimensional regularization is in this case already a 
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renormalization prescription. However, by doing so, a divergent quantity has been re- 
moved which is formally eliminated by setting — Z3 equal to this quantity. Regularization 
and renormalization in even dimensions can be performed using the MS-prescription [7T] . 
This procedure yields the same result as the momentum subtraction scheme. The range 
allowed for g in ()D.2|1 permits Igt not only to have a divergence of logarithmic or linear 
order, but also quadratic or cubic divergences in even or uneven dimensions, respectively. 
Using a subtraction scheme, it would be necessary to include the next term in the Taylor 
expansion. On the other hand, dimensional renormalization directly yields 

J _ qICa 2^~^n4^(rf - 3)rf + 21+^^(1 + g - dg))V{^ - g)V{-g)V{2g) 
(47r)l {2-d + 2g){d + 2g)V{^^-2g)V{g)V{^^ + g) 

Note that this expression becomes negative already for values allowed by ()D.2|) . e.g. for 
(yf > 3/4 in three dimensions, thus reducing the allowed range and leading to the plots in 
figure 



D.2 Massless Higgs 

If the Higgs is massless, then the Higgs-loop in the gluon equation ()4.1ip is no longer 
suppressed in the infrared and can in principle be as leading as the ghost-loop or even 
dominate. Thus three solutions can be found in the infrared. As the arguments from 
section 14.2.21 are still valid, all must satisfy t < —1. 

In the first case, the Higgs and the ghost exchange their roles. The ghost behaves 
then like a tree- level particle in the infrared with = 0, while the Higgs diverges in the 
infrared and drives the gluon. In the second case, / = g, and the gluon is driven by the 
combination of both loops. In the third case, the ghost and the gluon behave as in the 
massive case and / = 0: The Higgs behaves as a massless tree- level particle in the infrared. 
In all cases gluons are confined. The first solution may correspond to a phase structure 
similar to a Higgs- Anderson phase, but different to the conventional 'perturbative' Higgs 
phases. The physical interpretation of the second scenario is lacking and the properties of 
this phase would be quite pecuhar. The last case has similar physics as the massive case. 
Note that in the first two cases gauge symmetry is not necessarily intact any more, and 
may be broken by the interaction with the Higgs field. This may restrict these phases to 
theories which do not emerge from a pure Yang-Mills theory, as this would be difficult to 
reconcile with Elitzur's theorem |7SI. Hence the Gribov condition p.lUp will be dropped 
for the first two cases for now. As these cases offer a different structure than the massive 
case, they will be investigated a little more closely. 
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D.2. Massless Higgs 



After employing a renormalization prescription in the Higgs equation as in the case of 
the ghost equation, the infrared system becomes 

^ = Z,{l-5g{)+y-"^-^-'AgAJGT{y) 

^9 



y'-^-3-'AMHT{y). 



Keeping c? = 3 for simphcity, and noting that in both cases 



and 



i-HT 



l-GH 



^ = -^(i + ^) (D.3) 



- 11/ - 15/2 + 14/3)r(-2/) 



(/ - 3 + 32/2 + 44/3 + i6/4)r(i - 21) 
2-2+4Z^; - 1 + /C)r(-2/) sec(2/7r) sin2(/7r) 



(/-l)(4/-l)/r(|-2/) 
the conditional equation for the first case is 

(1 + /)(1 + 2/)(3 + 2/)(/ - 1 + /C)(sec(2/7r) - 1) 



1 



16(/-l)/(8 + /(l + 2/)(7/-ll)) 



This has no solution atC = linl/4</< 1/2, but has e.g. the solution / = 0.362568 at 
^ = 3. This already indicates that in these situation the allowed projections differ from 
the normal case. Since for a broken gauge symmetry the STI ()2.18|) may not be valid 
anymore and the gluon may acquire longitudinal components, this equally well indicates 
that the current truncation is inadequate for this case. 
In the second case, the condition for a solution becomes 

. _ Igg A\ Igh _ A^ Iqg Igh 
~ ^ Al I„T ~ AlU ^ ' 



'GT i-HT ^HT J-HT 



yielding 



Igh Igt 



Igt Iht — IggIht 

(3(C - 2) + 2^(59 - Ag\C - 1) + C - 6^(10 + C)) + 64(^ - l)gcsc\g'K)) 

%{l + g){l + 2g){?, + 2g){g-l + gC) 
■{S + g{l + 2g){Jg-ll)). 
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Only one infrared coefficient is left undetermined, thus it is likely that here the infrared 
behavior will be linked stronger to finite momentum solutions to allow for compensation 
both in the Higgs equation and the ghost equation in the infrared. Note that again a 
solution does not exist for ( = 1 for 1/4 < g = I < 1/2, but it exists for C = 3 with value 
I = g = 0.340408. Thus also here the validity of the truncations made is debatable. 

Nonetheless, these investigations show that additional massless degrees of freedom 
may change the infrared properties significantly. Especially in QCD it turns out that no 
solutions have been found yet with Nf > 4 massless quarks ,51j in this approach, despite 
the possibility of chiral symmetry breaking. 



Appendix E 

Perturbative Expressions 



This appendix contains the perturbative calculations to leading order in the dimensionally 
reduced theory of chapter El Replacing all full quantities in the truncated DSEs ()4.9M.lljl 
by their tree-level values, i.e. setting all dressing functions equal to one and using the tree- 
level vertices (16.14116.191) instead of the full vertices, the standard perturbation theory to 
one-loop order is obtained. It is then possible to calculate the leading-order perturbative 
dressing functions. For the ghost, this leads to 

G(p)-^ = 1-^ (E.l) 

to leading order, independent of the presence of a Higgs. Note that a Landau pole at 
p = 0(5'!) is present. 

The calculation of the gluon self-energy is a little more complicated, since each single 
contribution is linearly divergent, although the sum is finite in three dimensions due to 
the STI ()2.18j) . Corresponding problems are most easily circumvented by contracting the 
gluon equation with the Brown- Pennington projector ()4.3|) . Performing the calculation in 
pure Yang-Mills theory yields 

= 1 - IM^. (E.2) 



64p 



Including the Higgs yields 



Z{pV = 1 - + ^ ( - esc-^ I Jl + ^ 1 1 • (E.3) 




llgCU ^ qICa / ^rrih _ p" + Ami 
QAp IQixp \ p p^ 

Since non-perturbative effects already arise at order g^/p^ in the present truncation scheme, 
the only interesting part is the one for p ^ g^, rrih, leading to 
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Finally, the Higgs self-energy is 

= .^^^sJC^rn^^oJC.! i / rni+^ _ ^\ ^^^^^ ( _ ,„A , 
p'^ p Airp p Zirp \\p J \ y vrih + p j I 

(E.5) 

where the third term is the tadpole contribution. The leading contribution is 

H{pr = 1 - (E.6) 

The coupling constant h is determined by requiring that the sum of the tadpole kernels 
already generates a finite integral, i.e. independent of the regularization scheme. This is 
obtained by requiring ()4.4|) . 

Note that, in three dimensions, resummation produces effects only from order on. 
By dimensional arguments alone, only tree-level expressions in the loops can contribute 
at order g^. Therefore (jE.lj) . ()E.4j) and ()E.6|1 already constitute the resummed solutions. 
This is confirmed by the numerical calculations presented in section 

Note further that this also ensures that gauge symmetry is intact to one-loop order 
in the regime of applicability of leading-order resummed perturbation theory. The viola- 
tion of gauge symmetry in this approach is not stronger than in ordinary leading-order 
perturbation theory. 



Appendix F 
Numerical Method 



In order to numerically solve the DSEs, the method described in [HHl HOI lEHl HUtij was 

extended and improved in detail. The dressing functions were split into three parts. For 
sufficiently small momenta of the order of 10~'^g^, they are replaced by their infrared 
behavior ()4.14M.l'Bj) . For sufficiently large momenta, of the order 10^ g^, they are replaced 
by their leading order perturbative form due to ()E.1|) . ()E.4|) and ()E.6|) . In the remaining 
part, the dressing functions were separated into an analytic factor, which interpolates 
between the analytic infrared and ultraviolet behavior, and a modification function. The 
logarithm of this modification function is then expanded in Chebychef polynomials. This 
factorization is not necessary for the g = 1/2-solution, albeit it lowers the needed CPU 
time significantly. However, it was not possible to find the other solution branch prior to 
this improvement. The analytic interpolation functions are 

gl + q 



The coefficients of the Chebychef polynomials are determined using a global rather than 
a local Newton method |lU7j. This made it possible to solve the equations without any 
further prior knowledge about the solution. For a solution, it is required that relations 
fl4.Hl|l and ()4.21|1 are fulffiled at least up to 10~^ besides fulffiling the system of equations 
(14.9114. 11]) to a much higher precision. 

It turns out that simultaneous ffis of the Chebychef coefficients and Ag, A^, and Ah 
are not sufficiently stable. Hence an iterated procedure has been used where ffist a global 
Newton method was employed and then new infrared coefficients were extracted from the 
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requirement that the functions have to be continuous in the infrared. This was iterated 
until convergence of the infrared coefficients was achieved. A detailed account of the 
method will be given in |108j . 

To add the hard modes posed another problem. Due to the recursive nature of the 
subtraction prescription ()5.in|l of the spurious divergences, the Matsubara sum easily 
switches from logarithmic to quadratic divergence. The only way found up to now to solve 
this problem is to recalculate the renormalization constants in each Newton iteration step. 
By this in principle each Newton step again starts as a first step, but already with quite 
good starting values. This is sufficient to stabilize the algorithm and provides the solutions 
found in chapter El 
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